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ABSTRACT 

FACULTY OF PHYSICAL AND APPLIED SCIENCES 
SCHOOL OF PHYSICS AND ASTRONOMY 

Doctor of Philosophy 
PHENOMENOLOGICAL ASPECTS OF THE E(6)S.S.M. 
by Jonathan Peter Hall 

The work in this thesis explores various phenomenological aspects of the i^gSSM 
with a particular focus on the inert neutralino sector of the model and on the dark 
matter implications. The E'eSSM is a string theory inspired supersymmetric 
extension to the Standard Model with an Eq grand unification group. The model 
provides a solution to the hierarchy problem of the Standard Model, provides an 
explanation for neutrino mass, and has automatic gauge anomaly cancellation. 

The inert neutralino sector of the E'eSSM and the dark matter that naturally 
arises from this sector is studied for the first time. Limits on the parameter space 
from experimental and cosmological observations relating to the inert neutralino 
dark matter are explored and the consequences for Higgs boson phenomenology are 
investigated. In plausible scenarios it is found that the couplings of the lightest 
inert neutralinos to the SM-like Higgs boson are always rather large. This has 
major implications for Higgs boson collider phenomenology and leads to large 
spin-independent LSP-nucleon cross-sections. Because of the latter, scenarios in 
which E'gSSM inert neutralinos account for all of the observed dark matter are now 
severely challenged by recent dark matter direct detection experiment analyses. In 
plausible scenarios consistent with observations from both cosmology and LEP the 
lightest inert neutralino is required to have a mass around half of the Z boson mass 
if it contributes to cold dark matter and this means that tan(/3) cannot be too large. 

A new variant of the E'gSSM called the i^eZf SSM is also presented in which the 
dark matter scenario is very different to the inert neutralino cold dark matter 
scenario and in which the presence of supersymmetric massless states in the early 
universe modifies the expansion rate of the universe prior to Big Bang 
Nucleosynthesis. The new dark matter scenario is consistent with current 
observations and the modified expansion rate provides a better explanation of 
various data than the SM prediction. The prospects for a warm dark matter 
scenario in the EeSSM are also briefly discussed. 
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Overview 



In chapter [T] we present an introduction to the SM with a focus on EWSB and other 
aspects relevant for subsequent chapters such as gauge anomaly cancellation and the 
invisible decay width of the Z boson. Motivations for extensions of the SM such as 
the hierarchy problem and neutrino mass are explored and various notation is fixed. 

In chapter [2] we motivate TeV scale softly broken supersymmetry as a possible 
extension to the SM. A summary of supersymmetric Lagrangians is presented and 
various notation is fixed. Further concepts such as grand unification and 
universality of soft mass parameters are introduced. 

In chapter [3] the E'gSSM is motivated and introduced. This chapter contains 
previous work that has been carried out on the subject of the -EgSSM and provides 
background information relevant for chapters O [U and [71 

In chapter m the subject of dark matter is introduced. Information about the 
thermal dark matter relic density calculation relevant for chapters O [H and [7] is 
provided. We also provide an introduction to thermal relic dark matter in 
supersymmetric models. 

Chapter [S] contains work that was first published in paper I. This work 
represents a first study of the inert neutralino sector of the -EeSSM and the dark 
matter that naturally arises from this sector. 

Chapter [6] contains work that was first published in paper II. This work 
represents a more in-depth study of the inert neutralino and chargino sectors with a 
particular focus on physics relating to the Higgs boson. In plausible scenarios it is 
found that the couplings of the lightest inert neutralinos to the SM-like Higgs boson 
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are always rather large. This has major implications for Higgs boson collider 
phenomenology and leads to large spin-independent LSP-nucleon cross-sections. 
Because of the latter, scenarios in which E'gSSM inert neutralinos account for all of 
the observed dark matter are now severely challenged by recent dark matter direct 
detection experiment analyses. In plausible scenarios consistent with observations 
from both cosmology and LEP the lightest inert neutralino is required to have a 
mass around half of the Z boson mass if it contributes to cold dark matter and this 
means that tan(/3) cannot be too large. 

Chapter [7] contains work that was first published in paper III with the 
exception of section 1731 which contains work that is original to this thesis. In this 
chapter a new variant of the E'gSSM called the £^6^2 presented in which the 

dark matter scenario is very different to the inert neutralino CDM scenario and in 
which the presence of super symmetric massless states in the early universe modifies 
the expansion rate of the universe prior to BBN. The dark matter scenario is 
consistent with current observations and the modified expansion rate provides a 
better explanation of various data than the SM prediction. In section 17.51 the 
prospects for a warm dark matter scenario in the i^gSSM are briefly discussed. 

Summary and conclusions are found in sections 15.51 16.41 and 17.61 and in chapter [8l 

Notation relating to Weyl, Majorana, and Dirac Spinors and to the doublet 
representation of S\J (2) is fixed in appendices |A] and |Bj 
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Chapter 1 

The Standard Model 



The SM is an effective QFT describing the known particles and their interactions 
with the known forces of Nature excluding gravity. It is not currently known how to 
construct a consistent theory that unifies quantum field theory with our current 
best understanding of gravity which is the classical theory of general relativity. One 
candidate for the fully quantum description of gravity describing Nature is string 
theory, but whatever the correct description corrections due to the effects of 
quantum gravity are not expected to become relevant unless the energies involved 
in a process approach the Planck scale Mp ~ 10^® GeV, or alternatively unless one 
wishes to consider length or time intervals as small as Mp^ . We therefore expect to 
be able to use QFT, neglecting the effects of quantum gravity, at energies far below 
the Planck scale. Whilst the general framework of QFT is not expected to be valid 
above the Planck scale, the SM is itself only an effective QFT and is expected only 
to be valid below roughly the TeV scale — the energy scale currently being probed 
at the LHC. The reasons for this are outlined in section [L6l 

The SM contains our current best understanding of the observed particles and 
forces excluding gravity. The observed mesons and baryons that we observe are 
bound states of SM quarks, which are charged under the strong nuclear force 
described by QCD and there are also SM leptons which are free fundamental 
particles such as the electron. In terms of forces, the SM comprises the strong force 
of QCD as well the GWS theory of EWSB describing both QED and the weak force 
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responsible for nuclear decay. The description of the SM given in this chapter is 
largely based on the one given in ref. [6]. 

1.1 Gauge Symmetry and Matter Content 

The SM is a Yang-Mills QFT with a gauge symmetry group 

GsM = SUi3),(g)SU{2)L0U{l)Y. (1.1) 

Is it is a direct product of the SU{3) gauge symmetry describing QCD and the 
SU{2) X U{1) gauge symmetry of the electroweak theory — the unified theory 
describing both electromagnetism (QED) and the weak nuclear force. The observed 
fermionic matter of the SM can be described by LH Weyl spinors in 3 + 1 
dimensions forming chiral representations of Gsm as shown in table [TTTJ A RH Weyl 
spinor may be expressed as a LH one using the CP conjugation (charge conjugation 
and parity) operation. The notation for spinors used throughout is explained in 
appendix |S1 

The SM also includes a fundamental complex scalar doublet field known as the 
Higgs doublet whose VEV is responsible for the spontaneous breaking of 
SU{2)l X U{1)y down to the C/(1)em of QED, as per the GWS theory of 
EWSB [tlElE], and for the generation of fermion masses. Although the evidence 
for ESWB is overwhelming (for a review see ref. |lOj), the mechanism for this 
symmetry breaking is currently unknown, although it must have the correct 
custodial symmetry leading to the observed mass relation between the heavy 
electroweak W and Z gauge bosons. The SM assumes the GWS theory in which 
electroweak symmetry is spontaneously broken by the VEV of a fundamental 
complex scalar field — the Higgs doublet dD [121 [IS [H] . The VEV of this field 
is also able to generate masses for all of the SM fermions. 

The three components of the fundamental (3) and antifundamental (3) 
antitriplets of SU (3)c are known as colours (red, green, and blue) and anticolours 
(antired, antigreen, and antiblue) respectively. Since the EWSB vacuum respects 



4 





SU{3)c 


SU{2)l 


C/(l)y 


LH quark doublet 


Ql 


3 


2 


+ 1/6 


RH down- type quark 




3 


1 


+ 1/3 


RH up-type quark 




3 


1 


-2/3 


LH lepton doublet 


Ll 


1 


2 


-1/2 


RH charged lepton 


^R 


1 


1 


+ 1 


Higgs doublet 


H 


1 


2 


+ 1/2 



Table 1.1: The SUi'S)c and 311(2)1 representations and the 
U{1)y charges (hypercharges) of the SM matter fields, as LH Weyl 
spinors, and of the SM Higgs doublet H. 



SU{3)c, redefinitions of the three colours and three anticolours by SU{3)c 
transformations does not change the description of the physics. The SU{2) l gauge 
symmetry, however, is spontaneously broken by the vacuum and it makes sense to 
label the components separately. We define the third direction of weak isospin 
such that the Higgs VEV is an eigenstate of with eigenvalue — V^- Here we use 
T"- for generators of a general SU (2) l representation and for the generators of 
the 2 representation specifically. Since the direction in SU (2) l space of the Higgs 
VEV defines which direction will be uncharged under the unbroken f/(l)EM; electric 
charge will then commute with T^. This choice of direction for the Higgs VEV is 
arbitrary and has no effect on the physics, since any other equivalent choice would 
be related by a SU{2)l gauge transformation that leaves the Lagrangian invariant. 
Choosing the Higgs VEV to be an eigenstate of r^, the upper and lower 
components of the doublet in the standard basis r" = a""/2 are then eigenstates of 
electric charge. We write the quark doublet Ql and lepton doublet L^, as 



Ql 




and 



Ll 



(L2) 



eL 



The upper component is an eigenstate of with eigenvalue +1/2 and the lower 
component an eigenstate of with eigenvalue — 1/2. The charge under the 
unbroken C/(1)em of a field can be written 



Q = T^ + Y, (L3) 
where is understood to stand for the relevant eigenvalue. This is because it is the 
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gauge transformation with this combination of generators — )• (1 + ida{T^ + Y))H 
that leaves the Higgs VEV {H) invariant, i.e. {H) is uncharged under this 
combination of generators which must then represent the unbroken U{1). 

With one copy of each of the fields listed in table 11.11 we can describe what is 
known as the first generation of SM matter. This comprises the strongly interacting 
up and down quarks — two colour triplet Dirac fermions that are formed from the 
four Weyl spinor colour triplets of one copy of Ql, ur, and da — and also the 
electron — a Dirac fermion formed from bl and en — and a LH neutrino u, or 
equivalently the RH antineutrino v'^. 

The Lagrangian of the SM for this first generation, including all possible 

J3: 



renormalisable, gauge invariant, and Lorentz invariant terms 



IS 



C = -J^"^"^V + t^^^''^M^*+'^Yukawa + /:Higgs, (1-4) 



where 



/:Yukawa = - h"" d^j^H^ Q ^ - u^j^H.Q ^ - h'^ e^j^H^ L ^ + c.c. (1.5) 

and Cmggs contains the gauge invariant kinetic term and scalar potential of the 
Higgs scalar field shown in section [331 ^jia^P^'i/'i is the gauge invariant kinetic 
term for all LH Weyl spinors ipi, with the relevant gauge covariant derivative for 
each field tpi. For the gauge kinetic term 

= d^At-d^A^ + g^'^^r'^AlAl (1.6) 

and the adjoint index a runs over all of the generators of Gsm- The gauge coupling 
can have a different value for each of the three simple subgroups of Gsm and 



^We do not address the strong CP problem and will consistently neglect terms of the form 
A'^^ Afiv The contribution to this term from electroweak gauge bosons is always a total derivative 
and has no effect on the observable physics. However, the contribution to this term from QCD 
eventually, after chiral matter phase rotations removing 5 of the 6 arbitrary complex phases of CKM 
matrix, has an independent and arbitrary coefficient. This coefflcient should be very close to zero in 
order for the theory to be consistent with the non-observation of CP-violating effects from the QCD 
sector, such as an electric dipole moment for the neutron, but theoretically the smallness of this 
coefficient is not explained in the SM in a natural way. This is known as the strong CP problem — 
an unsolved problem in particle physics. 
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the gauge group structure function /"^^^ vanishes when a, b, and c do not belong to 
the same simple subgroup. The dot stands for the SU{2) invariant contraction of 
two SU{2) doublets given in ([Ril) 



V 



tl 
il 



V 



t2 
i2 



ilt2 - 142 . (1.7) 



The vacuum state of the Higgs potential is supposed to spontaneously break 
SU{2)l X U{1)y so that classically expanding around the true electroweak vacuum, 
rather than H = 0, we can write 

H = {H) + (p, (1.8) 



where the Higgs VEV 



I 







(1.9) 



and is an eigenstate of with eigenvalue —1/2 ((5 = 0). 

It is important to note that the unbroken gauge symmetry of the SM does not 
allow for any fermion mass terms — either Dirac or Majorana. However, in the 
EWSB vacuum the VEV of the Higgs doublet will generate Dirac fermion mass 
terms proportional to the Yukawa couplings h in (jl.Sp for all fermions other than 
the LH neutrino 



^Yukawa = ( -'-!=: <^^B^ L T^^W T^^iJ^L + C"^' ) • (l-^O) 



In the SM the Higgs VEV generates Dirac masses for all of the observed Dirac 
fermions, but does not induce any neutrino masses. 

Below the EWSB scale, one can integrate out the and Z electroweak gauge 
bosons that acquire masses from EWSB and write an effective theory with the 
gauge symmetry SU{2>)c x C/(1)em- The Lagrangian for one generation contains the 
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EWSB-induced mass terms of (|1.1U|) . Each mass term couples a LH and RH spinor 
together into a massive Dirac state. With both the LH and RH component of each 
Dirac spinor taken together, each Dirac spinor forms a real representation of 
SU{3)c X U{1)em- The LH neutrino forms a real representation on its own since it 
is a singlet — completely uncharged under the effective gauge group. Since this 
effective theory is non-chiral, containing pairs of LH and RH spinors that are 
equally charged under the gauge group, it is invariant under parity P and charge 
conjugation C separately. At low energy left- and right-handedness are only 
distinguished fundamentally in weak nuclear decay processes which violate 
separately both C and P maximally since the massive bosons only couple to 
LH states. 



1.2 Gauge Anomaly Cancellation, Generations, and the 
Invisible Z Boson Decay Width as Measured at 
LEP 

In Nature we observe that there are in fact at least three complete copies, known as 
generations, of the Weyl fields listed in table 11.11 The particle content of three 
complete generations has now been directly observed and there exists evidence, 
outlined in this section, that indicates that there are no more generations of any of 
the SM representations beyond these three. We shall label these three generations 
with the Roman indices i, j, etc., with the notation for the particles in these 
generations as complied in table [L2l 



1.2.1 Gauge anomalies 

Anomalies are quantum mechanical effects that violate one or more symmetries of 
the classical Lagrangian. A gauge anomaly is a quantum mechanical effect that 
violates some gauge symmetry. QFTs with gauge anomalies are inconsistent since 
gauge symmetry is required to cancel the unphysical degrees of freedom of the 
massless gauge bosons — the longitudinal space-like and time-like polarisations. In 
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Basis 


Flavour 


iviass 


j-iXi Qowii qucLrK 


"Ll - "L 


"Ll - "L 






"L2 ~ *L 


LH bottom qucirk 


W — h' 


"L3 - "l 


RH down quark 


rl 

"Rl - 


— ^ 

- «R 


l-< 1 — 1 OT T* o n f'i' t \ ^ 1 n o 1^ 

lYXl Gtictllgc CJUclIK 




= ^i? 


rLXT UUlbUlil C^Udliv 




— h 


LH up quark 


«L1 = 


= Ul 


LH charm quark 


^L2 


= Cl 


LH top quark 


«L3 


= h 


RH up quark 


^iJl = 


= «R 


RH charm quark 


^_R2 


= ^R 


RH top quark 




= tR 


LH electron 


^Ll = 


= ^L 


LH muon 


^L2 " 


= Ml 


LH tau lepton 


^L3 ' 


= ^L 


RH electron 


^Rl '- 


= eR 


RH niTion 




= 


RII tau lepU)U 


'■m 


= 'r 



LH electron neutrino 
LH muon neutrino 
LH tau neutrino 
Light neutrino 1 
Light neutrino 2 
Light neutrino 3 



'^1 
, ,/ 



Table 1.2: The notation for the three generations of fermionic 
matter of the SM. Where the flavour and mass eigenstate eohimns 
are combined the flavour and mass eigenstates are equal by def- 
inition. In the down quark sector the mass eigenstates are then 
rotated with respect to the flavour eigenstates by the CKM ma- 
trix. In the neutrino sector the mass eigenstates are rotated with 
respect to the flavour eigenstates by the PMNS matrix, which is 
analogous to the CKM matrix of the quark sector. The CKM ma- 
trix is relatively close to the identity, whereas the PMNS is close to 
tribimaximal form, meaning that the mass and flavour eigenstate 
bases are very different from each other. 
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4-dimensional QFTs such as the SM, gauge anomahes arise at one-loop level via 
triangle diagrams of the form 



oc ^"''^ = tr[T"{^^^^}], (1.11) 



where is the group generator corresponding to the adjoint index a of the gauge 
boson labelled a. The trace is a sum over all LH Weyl spinors running around the 
loop, each in some representation, and also a trace over the generator indices of the 
relevant representation. The anticommutator comes from considering each Weyl 
fermion running around the loop in both directions. Such diagrams must sum to 
zero for all combinations of different gauge boson external legs in order for the QFT 
to be consistent. 

Instead of any given RH spinor in the representation r, one may consider the CP 
conjugate state which is a LH spinor in the representation f (see appendix This 
gives a contribution to A""^"^ equal to tr [T^{T^, T^}]. Since = —Tf^ this 
contribution is in fact equal to -tT[T'^{T'^. ,T^}], which is minus the contribution 
from a LH spinor in the representation r. So, if one has an equal number of LH 
spinors in each of the representations r and r of the entire gauge group, or 
equivalently an equal number of LH and RH spinors in the representation r, then 
these states, which taken together from the real representation r © f, do not 
contribute to the gauge anomaly. If the representation r is explicitly real then this 
condition is automatically satisfied. This is the case in the SU{3)c x C/(1)em 
effective theory. 

However, since the SM is a chiral theory, meaning that the fermionic matter 
cannot be written as above in terms of their representations under Gsm, it takes 
some more work to compute the gauge anomalies in (jl.lip associated with the 
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various combinations of SM gauge bosons. Eventually one computes that if one has 
a complete generation of the chiral fermions listed in table 11.11 then all of the gauge 
anomalies do indeed cancel, but this is not the case for just the leptons or for just 
the quarks separately [B|. We therefore conclude that the SM is gauge-anomaly-free 
as long as it contains only complete generations of matter, i.e. it contains the same 
number of generations of quarks and leptons. 

1.2.2 The effective number of neutrinos contributing to the 
invisible Z boson decay width 

The best evidence for the number of generations comes from the number of 
neutrinos as inferred from the invisible decay width of the Z boson measured at 
LEP, i.e. the partial decay width of the Z boson into particles that do not show up 
in the detector. The effective number of neutrinos at LEP A'^^ is defined by 

r(Z ^ invisible) = N^i^T{Z vv) (1.12) 

where the decay width on the left is measured and the decay width on the right is 
calculated assuming that is a massless LH neutrino, v is the corresponding RH 
antineutrino v^. The result from LEP [15j is 

iV^^g^P = 2.984 ±0.008 (1-sigma), (1.13) 

leading to the conclusion that the number of neutrinos is 3, this being the closest 
integer to the central measured value and 2-sigma away. This means that there is 
no fourth generation neutrino with a mass lower than about half of the Z boson 
mass. This is taken as evidence for there being only three generations of leptons 
and, since the numbers of generations of quarks and leptons should be equal in 
order to have gauge anomaly cancellation, only three generations of quarks also. 
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1.3 The Higgs Potential and GWS EWSB 



In the SM EWSB is caused by the non-zero VEV of a single Higgs scalar doublet 
H. The most general gauge invariant and renormalisable form of ^Cniggs appearing 
in (fOI) is 

-Cffiggs = {V^'H)^{V^H)-V{H), (1.14) 

with 



V{H) = m'^H^H + X[H^Hf. (1.15) 

The parameter A must be positive in order for the Higgs potential to be bounded 
from below. If the mass-squared parameter m? is also positive, or zero, then 
classically V{H) has a minimum at = 0. In this case H = Q \s the true vacuum 
and GsM remains unbroken. If, however, is negative, then the degenerate 
minima of V{H) occur on the surface given by 



2 



H^H = (1.16) 



points on which are related by arbitrary SU{2)l gauge transformations. Using our 
earlier definition of the Higgs VEV ()1.9p we can then identify 



2 

-m 



X 



(1.17) 



In this case in the EWSB vacuum Gsm is spontaneously broken to 
SU{3)c X C/(1)em- SU{3)c remains unbroken since H is a singlet under this group. 
Since it is a scalar particle, it is also a singlet under the Lorentz group and its VEV 
therefore does not spontaneously break Lorentz symmetry. 

The covariant derivative acting on H is given by 

V^H = {d^-ig2W;T^ -y2g'B^)H (1.18) 
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where W"" are the three SU{2)l gauge bosons, r° are the generators of SU(2)l 
(T") for the fundamental doublet representation (r" = cr"/2), and is the gauge 
boson for U{1)y under which H has charge +1/2. (72 and g' are the SU{2)l and 
U{1)y gauge coupling constants respectively. Expanding the kinetic term around 
the Higgs VEV as in (jl.Sp we find that in the EWSB breaking vacuum we generate 
the gauge boson mass terms 



r - i!^ 

-L-Higgs ~ 2 4 ■ 



2 ^ 2- 



gliWlf + gliWlf + {g,Wl - g'BX\ + • • • . (1.19) 
It is useful to define the positively and negatively charged SU (2) l gauge bosons 



^iK^'K) (1-20) 



as well as the usual = it iT"^ . In the EWSB vacuum we read off that these 
particles have a mass 



V 



The (correctly normalised) mass eigenstate 



\^{g^Wl-g'B^), (1.22) 



where 



-9 = \^9i + 9", (1.23) 



has a mass 



mz = 9^, (1.24) 



leaving the orthogonal combination 



A, = k9'W'^+92B^) (1.25) 
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massless. This is the photon — the gauge boson of the unbroken L'"(1)em which 
corresponds to the the unbroken combination of generators Q = T^ + Y. 

The general covariant derivative, neglecting the SU (3)c gluon terms, 



-D^ = d^-ig2W^T''-Yg'B^ (1.26) 
can then be written in terms of the gauge boson mass eigenstates as 



. 52_ _ g2^g^ z _ ieQB^, (1.27) 



where 



e = 



929 



cw = cos{'dw) = ^) and (1-28) 



9 '' 



implying that mz = mw / cw ■ 

In the unbroken theory the Higgs complex scalar doublet has four real degrees of 
freedom. After EWSB the three massive gauge bosons each acquire one extra 
degree of freedom from the complex scalar doublet, corresponding to the 
longitudinal polarisation that exists for a massive vector boson, but not for a 
massless one. The remaining one degree of freedom belongs to a real scalar, known 
as the SM Higgs boson. 

We can work in the unitarity gauge in which the three Goldstone modes of the 
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Higgs doublet are set to zero and we expand around the EWSB vacuum 








(1.29) 



V 



v + h 



where h is the (canonically normahsed) real scalar known as the Higgs boson. Since 
in the basis that we have chosen v appears in the real part of the lower component 
of H, this is also the direction corresponding to the massive boson state h. The 
other directions are flat and correspond to the massless Goldstone modes whose 
degrees of freedom contribute to those of the massive gauge bosons. 

We can expand V{H) in the unitarity gauge in order to find the mass of the 
Higgs boson h. We find 



Like the induced fermion masses, the mass of the Higgs boson itself is proportional 
to the Higgs VEV v, but also to an unknown coupling constant A. The value of 
V = 246 GeV is determined from the masses of the and Z bosons, but although 
in the GWS theory this combination of < and A > are determined, the 
individual values of these parameters are not determined unless the Higgs boson 
mass is known. At the time of writing the Higgs boson is currently the only particle 
of the SM yet to be discovered. By looking for the process e~^e~ — t- Zh at LEP, a 
lower limit on the SM Higgs boson mass of 114.4 GeV is obtained [10]. Recent LHC 
analyses from CMS [l6] and ATLAS [13 US IHl EQl EI] between them exclude the 
existence of a SM Higgs boson with a mass between 145 and 288 GeV or between 
296 and 466 GeV at a 95% confidence level. 



V{H) = —h^ + —-h^ + ... 




(1.30) 



from which we read off a mass-squared for the real scalar h 



ml = -2n? = 2Au^. 



(1.31) 
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1.4 Induced Dirac Fermion Masses, the CKM Matrix, 
and Neutrino Mass 



Including all three generations of SM matter, (jl.lOp becomes 



■Yukawa 



^/2 



V 




FA^Lj + ^^<^L, + ^^4%- + c-c.) + • • • . (1.32) 



The Yukawa coupling matrices h^, hF , and may be made diagonal if one 
performs unitary transformations on the fermion fields in flavour space, i.e. 



for each of the fields dm, dn, URi, un, ejn, and CLi- Specifically for the LH quarks 
we write 



This is biunitary diagonalisation of each of the three Yukawa coupling matrices and 
the basis where these matrices are diagonal is the mass eigenstate basis. In the 
gauge invariant fermion kinetic term in the Lagrangian these transformations leave 
everything invariant apart from the couplings of the fermions to the heavy 
bosons coming from the covariant derivative (|1.27p . If one begins with non-diagonal 
Yukawa coupling matrices and then transforms to the mass eigenstate basis that 
diagonalises them, the Lagrangian term coupling quarks to bosons transforms 



(1.33) 



and u^i^UijULj. 



(1.34) 



<^^^'^W+d^^ + c.c. ^ ulia^^W+V^^d^^ +C.C. 




(1.35) 



producing a non-diagonal, unitary flavour mixing matrix 



V = u^^u- 



(1.36) 



known as the CKM [22l [23] matrix. 

As well as the mass eigenstate basis we also define a flavour eigenstate basis in 
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which the coupHngs to the heavy are diagonal. By convention we choose the 
up-quark flavour basis to be equal to the mass eigenstate basis. The down-type 
quark flavour basis is then 

du = Vi.dy (1.37) 

This convention is summarised in table 11.21 The most general form of V can 
contain three angles and six complex phases. However, complex phases in (|1.33p 
cancel out of (jl.32p and therefore complex phases in (jl.34p can be defined to 
remove five of these six phases. (One of the six phases in ()1.34p can be 
parametrised as an overall phase for all six transformations which cancels out of the 
right hand side of (jl.35p .) The SM CKM matrix can therefore be parametrised by 
three angles and one complex phase. This complex phase, which is responsible for 
CP- violating effects, is quite small. Although the CKM matrix is close to being the 
identity, the flavour eigenstates of the quark sector are not quite equal to the mass 
eigenstates. This means that there is a non-zero probability amplitude for a 
boson to couple together quark mass eigenstates of different generations. The 
bosons therefore contribute to FCCCs. The Z boson and the photon (and also the 
gluons) do not contribute to FCNCs since the transformations ()1.33p leave terms 
coupling neutral bosons to fermions, from (jl.27p . invariant. 

1.4.1 Neutrino mass and the type-I see-saw mechanism 

In the SM, which does not include RH neutrinos, LH neutrinos are exactly massless 
since both explicit mass terms and renormalisable terms coupling them to the Higgs 
VEV are forbidden by the gauge symmetry. However, in Nature we now know that 
neutrinos oscillate |10l [M] — a mechanism that requires them to have different 
masses, with the mass eigenstates being rotated with respect to the flavour 
eigenstates. There should be a mixing matrix for the lepton sector, analogous to 
the CKM matrix of the quark sector, known as the PMNS [25^ l26l [27] matrix. 

If we define the charged lepton mass eigenstates to be the eigenstates of flavour 
then the neutrino mass eigenstates will be rotated with respect to the flavour 
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eigenstates by the PMNS matrix. Neutrinos that are produced in some flavour 
eigenstate will then be in a superposition of mass eigenstates. If the mass 
differences involved are small enough, as they must be since the masses themselves 
are small, then the neutrino will propagate coherently as this superposition, but 
with each mass eigenstate component evolving at a different rate, causing 
interference. Therefore a neutrino that is produced as one flavour and propagates 
for some distance may, when it eventually participates in another charged weak 
current interaction, be measured as a different flavour with some probability. 

Whereas the CKM matrix is relatively close to the identity, neutrino oscillation 
data indicates that the PMNS matrix is close to tribimaximal form [28\ I24j. meaning 
that the mass and flavour eigenstate bases are very different from each other. 

One may in principle add to the SM model matter content some number of RH 
neutrinos ^ that couple to the Higgs field and lepton doublet, inducing Dirac mass 
terms after EWSB. Such RH neutrinos would have to be uncharged under Gsm in 
order for the term '^I'^H.Li to be gauge invariant. This in turn means that 
Majorana mass terms for the RH neutrino could also be added to the SM 
Lagrangian. These RH neutrino masses would be unrelated to EWSB — their scale 
associated with some new physics. Let us assume that these Majorana masses are 
much larger than the Dirac neutrino masses induced by EWSB. For one generation 
we may write a neutrino mass term 



where 971 is the RH neutrino Majorana mass, and m is the Dirac mass equal to 
some Yukawa coupling times v. For S!Jt ^ m there is one eigenvalue approximately 
equal to 9Jt and another approximately equal to m^/9Jt. In this case there is then a 
light mass eigenstate that is almost, but not quite, v and that has a Majorana mass 
that is suppressed relative to the EWSB scale. This is the type-I see-saw 
mechanism |29t [5U] . The principle still holds for three generations of RH and LH 
neutrinos, with three of the states arising from the 6x6 mass matrix having 




(1.38) 
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non-zero, but suppressed, masses. 



1.5 Baryon and Lepton Number Conservation 

The renormalisable SM Lagrangian, without RH neutrinos, is invariant under two 
extra U{1) global symmetries known as U{1)b and U{1)l, corresponding to baryon 
and lepton number conservation respectively. Quark fields Ql, dji, and uji have 
baryon number B = +1/3 and lepton number L = 0, with the CP conjugate 
antiquark fields having B = —^/'i and L = 0. The lepton fields Li and have 
-6 = and L = +1, with antileptons having L = — 1. Both B and L are conserved 
by the classical renormalisable Lagrangian, but are anomalous if gauged and are 
violated non-perturbatively. However, the global symmetry U{1)b-l^ corresponding 
to the conservation of the combination B — L, happens to be anomaly free in the 
SM if gauged and globally is conserved even non-perturbatively. 

U{1)b~l iSi however, broken explicitly by Majorana neutrino mass terms. 
Majorana mass terms for the light neutrino mass eigenstates imply (and are implied 
by) the existence of neutrinoless double beta decay [311 [32] (see figure II. ip — a 
process in which baryon number remains unchanged but lepton number is changed 
by 2. Experiments searching for neutrinoless double beta decay (see for example 
CUORE [33], EXO [M], GERDA [Sg, MAJORANA [36], NEXT [37], and 
SNO-|- [38j) will eventually determine the nature of neutrino mass — Dirac or 
Majorana. Because of the smallness of physical neutrinos masses, this 
J7(l)B_i-breaking effect if the neutrino is Majorana in nature would be 
corresponding rather small, with the physical Majorana mass appearing in the 
matrix element of any such process. 

To date no processes violating either baryon or lepton number have ever been 
directly observed. In addition to neutrinoless double beta decay, another example of 
such a process would be proton decay — the decay of a proton with B = +1 into a 
final state with B = 0. This is a process which, unlike neutrinoless double beta 
decay, needs not necessarily violate U{1)b-l- 
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Figure 1.1: A diagram for a neutrinoless double beta decay pro- 
cess induced by the existence of a Majorana neutrino mass. 



1.6 The Hierarchy Problem of the SM 



As previously stated, the Higgs boson mass in the SM is not determined by the 
other known parameters of the model. There do, however, exist various theoretical 
bounds p9]. Prom p.3ip we see that the Higgs boson mass-squared is proportional 
to the self-coupling A. The running coupling A in the loop-corrected potential is 
required to remain positive in order for the EWSB vacuum to be stable. For low 
values of the running coupling A the coupling decreases with increasing energy 
scale. Depending on the cut-off energy scale A that one requires the model to be 
valid up to the vacuum stability requirement puts a lower bound on the SM Higgs 
boson mass — a bound that increases with A [40J. At the same time, for larger 
values of A the coupling increases with energy. Too large values of A below some 
cut-off energy scale A therefore render the perturbation theory invalid. Here there 
exist uncertainties associated with the using of perturbation theory to try to assess 
where perturbation remains valid, but nonetheless the requirement that there is no 
Landau pole in A below A puts an upper bound on the Higgs mass — a bound that 
decreases with A. This large Higgs mass effect can also be seen non-perturbatively 
in Lattice calculations \4:1\ I42|. Furthermore, too large Higgs boson masses lead to a 
non-unitarity of the S-matrix for certain processes where unitarity is preserved via 
cancellations between divergent diagrams involving virtual Higgs bosons and 
divergent diagrams involving virtual longitudinal polarisations of massive weak 
gauge bosons such as WW scattering. Unitarity bounds should not be violated in 
renormalisable theories |43| . 
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Because of these various upper bounds the SM Higgs boson should have a mass 
below around the TeV scale in order for the theory to be valid. This however 
introduces a naturalness problem into the theory — the unnatural hierarchy 
between the EWSB scale and the Planck scale which is about 16 orders of 
magnitude greater. The reason why this is considered unnatural is because in the 
SM the Higgs doublet is a doublet of fundamental complex scalars. (They are 
fundamental scalars as opposed to composite scalars which would be expected only 
to appear in some effective theory of the constituent particles that one would only 
expect to be valid up to some energy scale associated with the confinement.) 
Fundamental scalars are a problem in non-supersymmetric theories, because their 
masses receive radiative corrections proportional to the masses of any particles that 
they couple to [441 SS] ■ Their self-energy Feynman diagrams are quadratically 
sensitive to the highest mass scales in the theory. 

For example, let us consider the one-loop contribution to the self-energy diagram 
of a fundamental scalar that couples to a fermion of mass mp with a Yukawa 
coupling Ap using dimensional regularisation with a mass scale parameter fj, 




(1.39) 



da 



(4vr) 



A2 



7 + 1 

e 



In 



47r/i" 
^F 



, (1.40) 



where Ap = rrip — a{l — oi)(f' and the number of dimensions d = 4 — e. This 
contribution contains a part that has a pole at the 4-dimensional limit e — )• and 
additional finite parts including a part proportional to mp. The coefficient of Ap in 
is proportional to the logarithm which is order one. 



If we define the renormalised scalar propagator to have a pole where the energy 
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equals the renormalised mass = Mp then we have 



= 



+ {Zm -1) + {Z- 



(1.41) 



where T? is the total correction to the scalar mass-squared in the propagator from 
loop corrections to and counter term insertions in the scalar propagator at one-loop 
order. Zm and Z are the scalar mass term and wavefunction renormalisations 
respectively and —lA^ is the total one-loop correction to the scalar self-energy 
diagram. 

If the SM were valid up to arbitrarily large energy scales, with no new physics 
existing at higher energy scales, then there would be no problem. The 
renormalisation constants may be defined to cancel the poles in e of the constant 
and coefficients and, since the theory is renormalisable, such poles would then be 
cancelled by counter terms at all orders in perturbation theory. The additional 
finite corrections would be at most of order the top quark mass. 

However, this is not the case. Even if no physics comes in earlier the SM cannot 
be valid above the Planck scale where contributions from quantum gravity should 
become important. This being the case, it is not clear that using dimensional 
regularisation and integrating momenta up to infinity is a valid prescription, but if 
one alternatively uses a cut-off regulator, with a momentum cut-off at some energy 
scale where the theory ceases to be valid, then one still obtains the generic result 

- lAl^ = -iXl [pole + Cpmj, + •••], (1.42) 

where Ci? is some order one coefficient and the 'pole' is now of order the cut-off 
scale squared. If the scalar also couples to another scalar of mass ms with a 
coupling constant As then this also gives a contribution of the generic form 

- lAl = iXs [pole -I- Cgms H ] . (1.43) 

Regardless of the physical interpretation, it is still the case that if these pole parts 
are cancelled by the renomalisation constants at this level then, since the theory is 
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renormalisable, the poles will still be cancelled by counter terms at all orders. This, 
however, cannot be said of any large finite corrections proportional to boson or 
fermion masses associated with some new high energy physics |45j . 

If one requires the renormalised mass to be much smaller than such large 
additional finite contributions to then one one may also define the 
renomalisation constants to almost completely cancel these additional finite terms, 
leaving the small desired mass, at fixed order. The unnaturalness arises when one 
then goes to higher order in perturbation theory. While the new poles that arise at 
this order will be exactly cancelled if one also includes all diagrams containing 
counter terms up the relevant order, new large finite corrections will also arise that 
will not in general be cancelled by the finite parts of the counter terms. The 
renormalised mass is therefore expected to be of order the largest mass scale in the 
theory unless the finite parts of the counter terms are carefully retuned at every 
order in perturbation theory. 

The SM Higgs mass is thus sensitive to any new physics that might exist at or 
below the Planck scale. Since we know that QFT itself is not expected valid at the 
Planck scale, it is unreasonable to assume that there is not some new physics at 
some scale far higher than the EWSB scale. If the Higgs boson couples to this new 
physics at all then the Higgs boson mass should be at least of order this scale unless 
one is willing to accept large tunings at every order in perturbation theory to make 
it such that the large contributions cancel, leaving a Higgs mass of order the EWSB 
scale. This is the hierarchy problem of the SM. 

If, for example, we include radiative contributions from RH neutrinos with 
masses of order 10^^ GeV, then the additional finite corrections to the Higgs boson 
self-energy of order 10^^ GeV^ should be tuned to almost cancel leaving a physical 
Higgs boson mass-squared 24 orders of magnitude smaller. 

This hierarchy problem leads us to conclude that some new physics must exist at 
or around the TeV scale to stabilise the Higgs mass. The most common theories 
motivated to solve the hierarchy problem have solved it by assuming that the Higgs 
boson is a composite, as in the case of technicolour theories; by assuming that the 
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Planck scale is in fact around the TeV scale, as in the case of large extra 
dimensions; or by assuming that the Higgs mass is stabilised by super symmetry, as 
in the theories that we shall introduce in chapters [2] and [3l All of these theories 
involve the existence of new physics at the TeV scale. 

1.7 Unsolved Problems in Particle Physics 

Although the main motivation for new physics, particularly at the TeV scale 
currently being probed by the LHC, is the hierarchy problem, there are many other 
questions left unanswered by the SM. We shall briefly mention some of them in this 
section. Neutrino mass has already been discussed, but many other questions about 
fundamental fermion mass also remain unanswered. Although the induced fermion 
masses are allowed in the SM, the Yukawa couplings are measured and not 
predicted. Theories that attempt to explain the sizes and values of these Yukawa 
couplings as well as the striking difference between the CKM and PMNS matrices 
are known as theories of flavour. These typically invoke some new symmetry known 
as flavour symmetry with the spontaneous breaking of flavour symmetry producing 
the observed patterns of Yukawa couplings (see for example refs. [36 l ITTl HH] ). 

In the SM there are also three independent and unexplained gauge couplings. 
Grand unification (introduced in section 12. 4p proposes that the SM gauge group is 
in fact the remnant of some larger spontaneously broken gauge group with a single 
gauge coupling. In such a scenario the SM gauge couplings, running up in energy, 
should unify to the same value at some energy scale associated the breaking of the 
larger GUT group. In the SM the couplings do not in fact unify, but they do in the 
supersymmetric models introduced in the next chapter. 

The baryon asymmetry of the universe is another problem. It is not known why 
the universe appears to be made almost entirely of matter and not antimatter. 
Although the SM technically satisfies the Sakharov conditions 09] — conditions 
required for the existence of baryogenesis processes that could have created this 
asymmetry — of baryon number non-conservation (non-perturbatively) and CP 
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violation (via the CKM matrix), the small CP-violating phase of the CKM matrix 
is not thought to be large enough to have been the origin of the observed baryon 
asymmetry. 

Although most of the baryonic mass in the universe is to some extent 
understood (being mostly due to QCD colour confinement rather than the Higgs 
mechanism), dark matter, discussed in chapter [H and dark energy are completely 
unaccounted for in the SM. 

The strong CP problem (see footnote 1 in section [TTT]) is another unsolved 
problem. 

The origin of the gauge group, matter representation, and number of space-time 
dimensions is also not understood, although the theory must be consistent with 
respect to gauge anomaly cancellation and the existence of stable atoms and stable 
gravitational orbits is obviously necessary for our existence. 
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Chapter 2 



Supersymmetry and Grand 
Unification 

The description of supersymmetry given in this chapter is largely based on the 
descriptions in refs. |5m 145). The notation for fermion spinors used is given in 
appendix lAl 

Supersymmetry is a symmetry relating particles of different spin. In a theory 
with some amount of supersymmetry each particle, possessing a given spin and 
other internal quantum numbers, necessarily comes as part of what is known as a 
supermultiplet — an association of particles that have different spins but all other 
quantum numbers the same. The particles in these supermultiplets are then 
transformed into each other by supersymmetry transformations that leave the 
supersymmetric Lagrangian invariant. The size of the supermultiplets describing 
the theory depends on the number of conserved supercharges J\f of the 
supersymmetry algebra. If the supersymmetry is preserved by the vacuum then the 
particles of different spin that make up a supermultiplet are all degenerate in mass 
as well as having all other quantum numbers the same. 

Plausible low energy models for physics beyond the SM can be constructed using 
M = 1 supersymmetry. This theory contains the following types of supermultiplet: 
A chiral supermultiplet containing a complex scalar and a LH (for a left chiral 
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supermuliplet, RH for a right chiral supermultiplet) Weyl spinor; a vector 
supermultiplet containing a spin-1 real vector and a Weyl spinor; and a graviton 
supermultiplet containing the spin-2 graviton and a spin-3/2 gravitino. Each of 
these supermultiplets separately contains the same number of physical bosonic and 
fermionic degrees of freedom. The (as we will see is necessary, multiple) Higgs 
scalars whose VEVs are responsible for EWSB must then be part of chiral 
supermultiplets containing the same number of fermionic degrees of freedom in the 
form of spin-i/2 fields. These fermions are known as Higgsinos. Each Weyl fermion 
matter field of the SM must be part of either a chiral or vector supermultiplet and 
in plausible models they are all contained in chiral supermultiplets. The scalar 
superpartners of the quarks and leptons are known as squarks and sleptons 
respectively. In super symmetric gauge theories the massless gauge bosons form 
vector supermultiplets along with Weyl fermions known as gauginos. 

Although not the original motivation for supersymmetry itself [51], the main 
motivation for what is known as TeV scale softly broken supersymmetry (see for 
example ref. [52j) is that it provides a solution to the hierarchy problem of the SM. 
In supersymmetric theories scalar self-energies do not have the quadratic sensitivity 
to high energy scales that are the origin of the SM hierarchy problem. The 
quadratic terms due to fermions in loops such as ()1.39p are cancelled by quadratic 
terms due to the bosons from the same supermultiplet. These boson terms have the 
opposite sign since they do not have the extra minus sign associated with a fermion 
loop. In another sense, the non-existence of this quadratic sensitivity comes about 
because the Higgs scalar itself is part of a supermultiplet and must remain 
degenerate in mass with its non-scalar superpartners. For A/" = 1 chiral 
supermultiplets the complex scalar must remain degenerate with the Weyl fermion. 
Since Weyl fermions do not have the quadratic sensitivity to high energy scales, the 
quadratic contributions to the scalar self-energy must necessarily cancel. 

Clearly the scenario described contradicts observation if supersymmetry is 
preserved by the vacuum since it invokes the existence of many new unobserved 
particles that are degenerate in mass with observed particles and have similar 
interactions. In realistic models supersymmetry, and this mass relation, must be 
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broken. We will see that the TeV scale soft breaking scenario provides a solution to 
the hierarchy problem, but predicts that there should be observable superpartners 
with masses not too far above the TeV scale. 



2.1 Superpotentials 

The renormalisable Lagrangian of an = 1 supersymmetric gauge theory is 
specified by specifying the gauge group, the gauge group representations of the 
chiral supermultiplets, and what is known as the superpotential W. The 
superpotential is a chiral object, being a dimension-3 holomorphic function of 
complex scalars from either purely left or purely right chiral supermultiplets. Here 
we will work purely with left chiral supermultiplets as is canonical. 

Let a supersymmetric gauge theory contain left chiral supermultiplets, labelled 
with i, each containing a complex scalar (pi and a LH Weyl spinor -0^. Furthermore, 
let the superpotential 



1 1 

VV = -mij(j)i(j3j + -Xijk(l)i4'j4'k, (2.1) 



with rriij = niji and Xijk similarly symmetric in all of its indices. Terms in the 
superpotential with mass dimension greater than 3 are non-renormalisable. The 
renormalisable, supersymmetric, gauge invariant Lagrangiarj^ is then 



1 



2 



where 



dW 1 

Fi = -g— = -rriijCpj - -\ijk(t)j<t'k (2.3) 



^We do not address non-renormalisable operators in supersymmetric theories. Although the 
effects of dimension-5 operators are interesting and potentially important, they have not been sys- 
tematically studied in the i^eSSM. 
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and 



D' 



(2.4) 



The gaugino A"" with adjoint index a is a LH Weyl spinor, so A"''^ is a RH Weyl 
spinor. If the gauge group is a direct product of simple subgroups then the gauge 



The MSSM is minimal in the sense that it introduces as few new particles as 
possible to the particles of the SM. To this end one begins by simply assigning all of 
the fields in table fTTT] to left chiral supermultiplets. This immediately creates a 
number of problems and in the MSSM these are solved in a way that introduces as 
few new fields as possible. 

Firstly there are two problems related to the assigning of the Higgs doublet H to 
a chiral supermultiplet, but both have the same solution. The first of these 
problems is that the Weyl fermion superpartner of the Higgs scalar doublet 
contributes to the gauge anomaly (jl.lip . The inclusion of this field gives extra 
non-zero contributions to gauge anomalies and therefore makes the gauge theory 
anomalous. The second of these two problems is that, since the superpotential must 
be a holomorphic function of complex scalars from purely left (by convention, 
alternatively right) chiral supermultiplets, superpotential terms coupling H to 
down-like squarks and charged sleptons are forbidden by the U{1)y gauge 
symmetry. This means that the mass inducing couplings to down-like quarks and 
charged leptons that appear in (jl.Sp cannot be present in the supersymmetric 
Lagrangian. The minimal solution to both of these problems is the same and it is to 
have two Higgs scalar doublets as in table [2Tl 

The extra contributions to gauge anomalies from the Higgsinos then cancel since 
their charges are opposite and together they from a real representation of Gsm- 
(We choose to write all SU{2) antidoublets as doublets since they are equivalent, as 



coupling constant g^^^ 



can have a different value for each of these subgroups. 



2.2 The Matter Content of the MSSM 
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Supermultiplet 


Boson 


Fermion 


SU{3)c 


SU{2)l 


U{1)y 


LH quark doublet chiral 


Ql 


Ql 


3 


2 


+ 1/6 


LH down-type antiquark chiral 


dl 


d^R 


3 


1 


+ 1/3 


LH up-type antiquark chiral 




Ur 


3 


1 


-2/3 


LH lepton doublet chiral 




Ll 


1 


2 


-1/2 


LH charged antilepton chiral 




4 


1 


1 


+ 1 


Down-type Higgs doublet chiral 


Hd 




1 


2 


-1/2 


Up-type Higgs doublet chiral 


Hu 


Hu 


1 


2 


+ 1/2 


Gluon vector 


G^' 


Gluino G 


8 


1 





SU (2) L gauge vector 




Wino W 


1 


3 





U{1)y gauge vector 




Bino B 


1 


1 






Table 2.1: The SU{3)c and SU{2)l representations and the 
U{1)y charges of the supermultiplets of the MSSM. 



shown in appendix [Bl) The contribution to the gauge anomaly due to gauginos is 
automatically zero since gauginos are necessarily in the adjoint representation 
which is real. Since Ll and Hd have the same quantum numbers one might think 
that a more minimal solution would be to declare these fields to be part of the same 
supermultiplet, but in practice such models prove unrealistic. 

The most general renormalisable and gauge invariant superpotential containing 
the fields in table 12.11 is 

W = fjiH,.H^ + h%ul,H^.QL^ 

- hf^d%,H,.QL, - hf^e%,H,.LLk + AW, (2.5) 

where 

+ Ct''LL^■Hu + l^^Jf^^'R^d'RA,. (2.6) 



If both of the Higgs scalars acquire VEVs 



(Hd) 




and {Hu) 



V2 



(2.7) 
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then this superpotential yields Dirac mass terms equivalent to those in (|1.3'2p 



>C Yukawa = \= ( hf'vddLd r- + hY.VuuLu 



V2 



+ hrjVdel,eLj + c.c.] +■■■ . (2 



In order for the and Z boson masses to be the same as their SM values we 
require 

= vj + vl (2.9) 
We therefore define an angle /3 such that 

tan(/3) = — (2.10) 

Vd 

=^ Vd = V cos(/3) and 
Vu = ■usin(/3). 

The Yukawa coupling matrices /i^, hF , and must then be multiplied with 
respect to those of the SM by factors of 1/ cos(/3), 1/ sin(/3), and 1/ cos(/3) 
respectively. With increasing tan(/3) the hierarchy between the top and bottom 
Yukawa couplings is lessened, but all of these Yukawa couplings are greater than 
their SM values for all angles (3. 



2.2.1 i?-parity 

The terms in AW, however, are dangerous since they all violate either lepton or 
baryon number conservation. Most importantly they lead to Lagrangian terms that 
allow protons to decay into final states with zero baryon number. These decays are 
mediated by squarks and require both baryon and lepton violating terms from AW 
(see figure [2T]) . The most common solution in the MSSM is to impose an 
additional discrete 1^2 symmetry on the fields in the superpotential. This is known 
as i?-parity, which we denote ^2^, and is defined such that the Higgs scalars are 
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even and all of the matter scalars (squarks and sleptons) are odd. When interpreted 
as a symmetry of the Lagrangian, angular momentum conservation implies that the 
fermionic superpartners have opposite i?-parity. The gauge bosons must be 
Zg^-even and the gauginos are then Z2^-odd. Oddness under 1^2^ is the meaning of 
the tilde over the squark, slepton, Higgsino, and gaugino fields. 

The imposition of forbids all of the terms in AW and makes the 
renormalisable MSSM Lagrangian invariant under global U{1)b and U{1)l. In the 
SM only gauge invariance is required in order for these to be global symmetries 
since the squarks the sleptons do not exist. The 1^2^ -odd particles, denoted with a 
tilde thoughout, are known as the supersymmetric particles and the imposition of 
1^2^ means that the lightest supersymmetric particle is absolutely stable. It is 
therefore the case that in the MSSM a discrete symmetry imposed in order to 
prevent rapid proton decay also leads to the existence of a new stable particle that 
may be a plausible candidate for dark matter. 



2.2.2 The /i problem of the MSSM 

The other problem with W in ()2.5p is that it contains the bilinear mass term //. 
This is a supersymmetry respecting parameter that a priori has no relation to 
either the EWSB or supersymmetry breaking scales. The problem is that in order 
to achieve EWSB, with v of the correct magnitude, the parameter ^ should be of 
order the EWSB scale. In the MSSM as written, however, it is not clear why it 
should not be either of order the Planck scale or zero. The /i problem refers to a 
fine-tuning that has to be imposed on the ^ parameter once. Supersymmetry does 
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at least mean that the parameter is stable at the EWSB scale under radiative 
corrections, even as it is not explained. 



2.3 Soft Supersymmetry Breaking 



Although unbroken supersymmetry is easily ruled out, even in the case of 



spontaneously broken supersymmetry the relationship 



(2.11) 



where A^^ is the mass-squared matrix for all real scalars and Ai'^ is the 
mass-squared matrix for all Weyl spinors in chiral supermultiplets, still holds at 
tree level in the absence of gauge anomalies j45j . It is trivially satisfied in the case 
of unbroken supersymmetry since the two real scalars and Weyl spinor from each 
chiral supermultiplet are degenerate. Because of this relation it has not been 
possible to create a realistic model for supersymmetry breaking in the MSSM 
without introducing extra physics. 

In order to create realistic models one usually invokes the existence of some 
other sector, known as the hidden sector, in which supersymmetry is spontaneously 
broken. In this hidden sector scenario the visible sector (containing SM matter, 
here the MSSM sector) does not itself cause spontaneous supersymmetry breaking, 
but supersymmetry breaking effects are communicated to it somehow from the 
hidden sector. 

In this hidden sector scenario it is useful to parametrise the kinds of 
modifications to the visible sector Lagrangian that spontaneous supersymmetry 
breaking in the hidden sector can cause in ignorance of the exact mechanism of 
supersymmetry breaking. It is therefore useful to list the gauge invariant mass 
terms that may be induced in the visible sector. Firstly one can have SSB masses 
for all scalars (f) of the form 
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If there is more than one copy of a scalar with the same quantum numbers, 
i.e. there is more than one generation, then there can be more structure to the mass 

doublets one can have 



-^L"^L/l, (2-13) 
and for the RH charged sleptons one can have 

-^Ri^lij^Rj^ (2-14) 

whereas the only possible soft scalar mass-squared term involving the down-type 
Higgs doublet is 

-HWhH,. (2.15) 



Note that and Hd do not have the same quantum numbers since Ll is Z^-odd. 
Secondly one can have SSB masses for gauginos 



M('') A^^^ A'' + c.c. 



(2.16) 



If the gauge group is a direct product of subgroups, then in general the gauginos 
associated with each subgroup can have a different gaugino mass. For example, for 
the SM gauge group one can have soft gaugino mass terms 



^ MgG'^'tG'a + AfsIV^'^^T^" + Mi^'^'^t^a + c.c. . (2.17) 



These mass terms always exist for gauginos since they are necessarily in real 
representations of the gauge group. They do not exist for chiral fermions. Thirdly 
one can have SSB trilinear terms. For each trilinear term that is allowed to appear 
in the superpotential 

>^ijk4>i4'j4'k (no sum on z, j, /c) (2.18) 



^Since sleptons are scalars there is no concept of slepton handedness. The handedness refers to 
the handedness of the fermionic superpartner. 
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one can have the soft supersymmetry breaking Lagrangian term 



\jkA\ijk(t^i(t>j<i>k (no sum on k), 



(2.19) 



where A has mass dimension 1. In fact it is tlie case that for any term tliat can 
appear in the superpotential one can have a corresponding SSB Lagrangian term 
that is equal to the superpotential term multiplied by some new supersymmetry 
breaking parameter with mass dimension 1. In the MSSM there then exists the SSB 
breaking term corresponding to the fj, term 



There terms are known as soft since they only involve new parameters that have 
positive-definite mass dimension. The supersymmetric relationships between the 
dimensionless couplings involving bosons and fermions that lead to the cancellation 
of the quadratic sensitivity of fundamental scalars to arbitrarily high scales are 
preserved. If these new parameters with the dimensions of mass are roughly of 
order some scale associated with SSB then the consequences are the following: 
Firstly, while the observed quarks and leptons only acquire masses proportional to 
the EWSB scale, the unobserved gauginos, squarks, and sleptons acquire masses 
proportional to the SSB scale, allowing their current non-observation to be 
naturally explained. Secondly, contributions to the finite radiative corrections to 
the Higgs boson self-energy will be at most of order the SSB scale since this is the 
scale of differences between bosonie and fermionic masses within supermultiplets. 
Therefore if the SSB scale is not too far above the EWSB scale then the hierarchy 
problem is still solved. It is therefore believed that if supersymmetry is the solution 
to the hierarchy problem then the squarks, sleptons and gauginos, while currently 
unobserved, should not have masses too far above the TeV scale and should 
therefore be discovered at the LHC. This is what is known as TeV scale softly 
broken supersymmetry. 

If the squarks and sleptons are present at or not too far above the TeV scale 



(2.20) 
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then these soft supersymmetry breaking terms can in general lead to FCNCs in 
contradiction with observation. This problem can be avoided if each of the 3x3 
soft scalar mass-squared matrices, such as those in ()2.13p and ()2.14p — namely 
"^Qij' "^dij' "^uij ' "^iij' ""^eij — proportional to the identity and if for each 
of the 3x3 Yukawa matrices — hY, , h£ , and hfi — all nine associated soft trilinear 
couplings are equal. Explicitly this means 

rrip^j = mjpS^j y F £ {Q,d,u, L,e} and (2.21) 
= Ag yi,j y G e {U,D,E}. (2.22) 

The SSB gaugino masses and trilinear couplings may also in general have large 
phases that lead to large CP- violating effects, again in contradiction with 
observation. The actual restrictions on and relationships between these soft mass 
parameters will depend on the nature of the SSB mechanism. 



2.4 Grand Unification 



Grand unification is the idea that just like U{1)em is a remnant of the 
spontaneously broken electroweak gauge symmetry group, so the SM gauge 
symmetry group Gsm is a remnant of some still larger group that is spontaneously 
broken by some mechanism at some GUT scale. The further assumption is that this 
GUT group should not be a direct product of simple groups, as Gsm is, but should 
itself be simple, with a single gauge coupling. Grand unified theories can offer an 
explanation for the charges of the observed SM particles and for the observed values 
of the three gauge couplings at low energies. The three gauge couplings would be 
equal to some single GUT group gauge coupling at the GUT scale and below the 
GUT scale, where the GUT symmetry is spontaneously broken, the three couplings 
would then run, with different beta functions, to the low energy values that we 
observe. Grand unification therefore makes the prediction that, running to higher 
energy, the three gauge couplings should unify at some scale. Since we know the 
values of the couplings at low energy, if we know the beta functions in some theory 
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then this prediction can be tested — the beta functions depending on the whatever 
new physics exists between the EWSB scale and the GUT scale. 

A U{1) gauge theory is of course invariant under a rescaling of the gauge 
coupling as long as the charges are also appropriately rescaled, but if the U{1) is a 
remnant from some spontaneously broken larger group then the U{1) charges of the 
particles, forming some representation under the larger group, will then be fixed. 

In the SM, with no new physics coming in above the EWSB scale, the couplings 
do not unify. They do, however, unify if one assumes supersymmetry spontaneously 
broken at the TeV scale. In the MSSM the scale of this unification is around 
10^® GeV (see for example ref. |45j ) . 

The smallest possible GUT group is SU{5) [53]. SM matter can arise from 10 
and 5 representations which decompose under 



SU{5) SU{3)c0 SU{2)l0U{1)y (2.23) 



as 



3 A, J ^ 1 1,2,- J l'-). (2.25) 



This is one generation of SM quarks and leptons as long as one uses the correctly 
GUT normalised U{1)y gauge coupling, which we can read off as 



51 = Jp- (2.26) 



The two Higgs doublets of the MSSM, however, do not form a complete 
representation of SU{5). If they are from 5 and 5 representations then one must 
explain why the colour triplets from these representations are not present at low 
energy while the Higgs doublets are. This is known as the doublet-triplet splitting 
problem. 
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2.4.1 Unification of SSB masses 



At the GUT scale in order for the soft gaugino mass terms to be gauge invariant 
under the GUT group all soft gaugino masses, like all gauge couplings, must be 
equal. This GUT scale soft gaugino mass is known as Mi/2- Grand unified gauge 
symmetry implies that these soft gaugino masses should be unified at the GUT 
scale, but the further assumption is sometimes made that not only should the 
unifications (|2.2ip and (|2.22p be imposed at the GUT scale for phenomenological 
reasons, but that all of the soft scalar masses in (|2.21|) should be equal to a unified 
soft scalar mass mo and that all of the soft trilinear couplings in ()2.22p should be 
equal to a unified trilinear coupling Aq. This this known as the constrained scenario 
or sometimes gravity mediated supersymmetry breaking. In the MSSM this 
constrained scenario is known as the cMSSM. 

Gravity mediated supersymmetry breaking assumes that supersymmetry 
breaking is communicated to the visible sector only by non-renormalisable 
operators that are suppressed by the Planck mass. If the supersymmetry breaking 
scale in the hidden sector is Mg then the visible sector SSB masses will be of order 
the gravitino mass ^.3^2 ~ M^/Mp [50j. A further assumption is that the 
non-renormaliasable operators should be completely flavour blind and that the SSB 
parameters should be unified at the Planck scale. In the constrained scenario, 
however, the unification relations are all applied at the nearby GUT scale. 

One success of this scenario is that if the soft Higgs masses start off equal to mo 
at the high scale then they will typically be driven negative by radiative corrections 
on the way down to the EWSB scale, allowing for EWSB if one also has an 
appropriate // parameter. Whether in a constrained scenario or not, this is known 
as radiative EWSB. The soft Higgs masses will be of the correct order of magnitude 
for EWSB since they will be of order the SSB scale. 
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Chapter 3 



The EqSSM 

The main theoretical shortcoming of the MSSM as a model describing TeV scale 
softly broken supersymmetry is the problem. The model also predicts the tree 
level result that the lightest Higgs boson must have a mass smaller than 
mz\ cos(2/3)|. Large loop corrections must then push the Higgs mass above the LEP 
limit in order for the model not to be ruled out. In practice when this is done the 
model is quite fine-tuned [54j. In light of the shortcomings of the minimal model, it 
is worth considering supersymmetric models that have a non-minimal structure at 
the TeV scale. 

The E'eSSM [55\ \5U\ [57] is a string theory inspired supersymmetric model based 
on an Eq GUT group. The low energy gauge group contains an extra U{1), called 
U{1)n, under which the RH neutrinos that arise in the model are not charged. This 
means that the RH neutrinos may acquire large intermediate scale Major ana 
masses. This choice, that the low energy gauge group is Gsm (8 U{1)n, defines the 
model. The U{l)]sf gauge symmetry is spontaneously broken at low energy by a 
SM-singlet field — charged under the extra U{1)n, but a singlet under Gsm- This 
field radiatively acquires a VEV which is naturally of order the SSB scale, meaning 
that there is Z' boson with an induced mass of order the TeV scale. This SM-singlet 
VEV also induces an effective /i parameter, also naturally of order the SSB scale, 
with the /i term of the MSSM being forbidden by the enlarged gauge symmetry. 

Although Eq is not a group without complex representations, complete 
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representations of Eq are nonetheless free of gauge anomalies. In the E^eSSM 
automatic gauge anomaly cancellation is thereby ensured by allowing three 
complete 27 representations of Eq to survive down to the low energy scale. These 
three 27s contain the three generations of known matter, however they also contain 
the VEV acquiring Higgs doublets and SM-singlet. This means that there are two 
extra copies of the Higgs doublets and SM-singlet in the low energy particle 
spectrum. Whereas in the MSSM the Higgs doublets do not form a complete 
representation of the potential SU{5) GUT group, in the i^gSSM supcrmultiplets 
with the quantum numbers of Higgs doublets are contained within each of the 
fundamental 27 representations of the GUT group that also each contain one 
generation of SM matter. 

In the E^eSSM only one generation of Higgs doublets and SM-singlets, defined to 
be the third, acquires the required VEVs and is known as 'active'. The other two 
generations, the first and second, of Higgs doublets and SM-singlets do not acquire 
VEVs and these are known as 'inert'. Furthermore, in the ii'eSSM it is assumed 
that the inert generations have suppressed Yukawa couplings to SM matter, 
suppressed due to some flavour symmetry. This means that new FCNCs from the 
enlarged Higgs sector are suppressed and also explains why the inert generations do 
not radiatively acquire VEVs. 

3.1 Gauge Symmetry and Matter Content 

The subgroups of the Eq GUT group may be written 

Eq D 50(10) ® C/(1)^ 

D SU{5)^U{1)^^U{1)^ 

D SU{3)c^ SU{2)l^U{1)y ®U{1)^®U{1)^. (3.1) 
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In the E^gSSM is spontaneously broken at the GUT scale directly to 
SU{3)c 5C/(2)l U{1)y U{1)n, where 



U(1)n = cos(i?)C/(l);^ +sin(i?)C/(l)^ (3.2) 

and tan(??) = \/T5. This is such that the RH neutrinos that appear in the theory 
are completely uncharged. Three complete 27 representations of Eq then survive 
down to low energy in order to ensure gauge anomaly cancellation. They 
decompose under the SU{5) U{l)]y subgroup as [58] 

27 ^ (10, ^ ) e (5, ^ 



5,-^]e(5,-^]e{l,-^](B{l,0]. (3.3) 



\/4oy V ' \/4o 

The first two terms contain normal matter, whereas the final term, which is a singlet 
under the entire low energy gauge group, contains the RH neutrino, or technically 
the LH antineutrino The second-to-last term, which is charged only under 
U{1)n, contains the SM-singlet S. The third generation SM-singlet acquires a VEV 

(53) = ^ (3.4) 

which, as we shall see, generates the effective term and spontaneously breaks 
U{1)n leading to a mass for the Z' boson. The remaining two terms contain the 
down- and up-type Higgs doublets Hd and Hu, but also contain SU{3)c triplets. 
These exotic coloured states are known as D and D — the antitriplet from 5 and 
the triplet from 5 respectively. Only the third generation of Higgs doublets acquires 
VEVs 



{H%) = ^ = ^cos(/3) and 



(H^^,) = ^ = -|sin(/3). (3.5) 
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Supermultiplet 


Boson 


Fermion 


r" 




^/5/3Q^ 




LH quark doublet chiral 


Ql 


Ql 


3 


2 


+ 1/6 


+1 


LH down-type antiquark chiral 




d^R 


3 


1 


+ 1/3 


+2 


LH up-type antiquark chiral 




Ur 


3 


1 


-2/3 


+1 


LH lepton doublet chiral 


Ll 


Ll 


1 


2 


-1/2 


+2 


LH charged antilepton chiral 




4 


1 


1 


+ 1 


+1 


LH antineutrino chiral 


N" 


Arc 


1 


1 








Down-type Higgs doublet chiral 


Hd 


Hd 


1 


2 


-1/2 


-3 


Up-type Higgs doublet chiral 


Hu 


Hu 


1 


2 


+ 1/2 


-2 


SM-singlet chiral 


S 


Singlino S 


1 


1 





+5 


Exotic colour antitriplet chiral 


D 


D 


3 


1 


+ 1/3 


-3 


Exotic colour triplet chiral 


D 


b 


3 


1 


-1/3 


-2 


Gluon vector 


CM 


Gluino G 


8 


1 








SU{2)l gauge vector 




Wino W 


1 


3 








U{1)y gauge vector 




Bino B 


1 


1 








U{l)]\f gauge vector 




Bino' B' 


1 


1 









Table 3.1: The SU{3)c and SU{2)l representations and the E'e 
GUT normalised U{\)y and U{1)n charges of the supermultiplets 
of the SgSSM. 



The charge assignments of the matter of the supermultiplets of the E^gSSM are 
summarised in table 13.11 



The low energy gauge invariant superpotential 



>Vo + Wi +W2, 



(3.6) 



where 



+ Wjk^RiHuj-QLk + h?jkd%iHdj-QLk + hfik^'ki^dj-^Lki (3-7) 



m = 9TjkDiQLj-QLk + 9ljkD,d'iiju'iif„ and 



(3.8) 



(3.9) 



It is now clear that the effective /i parameter is given by 



A333'S 

V2 ' 



(3.10) 
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generating the term fiHds-Hus in the superpotential. The fj, problem is solved since 
s is of order the SSB scale and A333 is perturbative. 

3.1.1 Discrete symmetries of the superpotential 

It should be noted that simply due to gauge invariance the superpotential of the 
-EgSSM is already invariant under the imposed on the MSSM provided that the 
exotic D and D bosons and the SM-singlet bosons are interpreted as being 1^2^ -even 
along with the Higgs doublets. The squarks and sleptons, including the RH 
sneutrinos, are Zg^-odd. The [/(l)B_L-violating terms of the MSSM superpotential 
that matter parity is invoked to forbid are never present in the renormalisable 
-EgSSM superpotential since they would violate the extra surviving U{1)n' gauge 
symmetry. Importantly, all of the C/(l)B_2,-preserving MSSM terms are gauge 
invariant with the exception of the /i term. Again the 1^2^-0(1(1 states are known as 
the supersymmetric particles and in the -EgSSM the LSP is automatically stable. 

In order for non-diagonal flavour transitions arising from the Higgs sector to be 
suppressed, the superpotential is assumed to obey an approximate symmetry 
known as Under this symmetry all of the fields in the superpotential other 
than ^3, Hfi^, and Hu3 are odd. It is this approximate symmetry that distinguishes 
between the active and inert generations of Higgs doublets and SM-singlets, with 
the inert generations having suppressed couplings to matter and not radiatively 
acquiring VEVs. This approximate symmetry suppresses Ajj^ couplings of the 
forms Aq,33, A3q,3, A33Q,, and Xa/S'y, where a,/3,^ £ indexing the inert 

generations only. Such an approximate symmetry, with a stable hierarchy of 
couplings, can be realised in E'gSSM flavour theories such as the one proposed in 
ref . [59] . The symmetry cannot be exact or else the lightest of the exotic coloured 
states would be absolutely stable. The existence of such stable coloured exotics 
contradicts observation [5U] . 

Although the C/(l)B-L-violating terms of the MSSM are forbidden by gauge 
symmetry, since the I2 cannot be exact another exact discrete symmetry must be 
imposed on the superpotential in order to avoid rapid proton decay caused by the 
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1 

+ 


1 

+ 


1 

+ 




'-'3' -"dS' -"u3 


1 


1 


1 


1 






+ 


+ 
















+ 


+ 







Table 3.2: The charges of the fields of the SgSSM supcrpoten- 
tial under various exact and approximate Z2 symmetries that the 
superpotential may or may not obey. 1i2^ is aheady a symmetry 
due to gauge invariance. Either or is imposed in order to 
avoid rapid proton decay. X2 is an approximate flavour symmetry. 
i e {1,2,3} and a € {1,2}. 

terms in the Wi and W2 that involve the exotic coloured states. There are two 
ways to impose an appropriate Z2 symmetry on W that lead to baryon and lepton 
number conservation. The first option is to impose a symmetry called under 
which only the sleptons, including the RH sneutrinos, are odd. In this case the 
superpotential is equal to Wo + Wi and the model is known as the ii^eSSM-I. U{1)b 
and U{1)l are symmetries of the renormalisable superpotential if the exotic 
coloured states D and D are, respectively, diquarks and antidiquarks, with B = ±2/3 
and L = 0. The second option is to impose a symmetry called under which both 
the sleptons and the exotic D and D bosons are odd. In this case the superpotential 
is equal to Wo + W2 and the model is is known as the E'gSSM-II. U{1)b and U{1)l 
are symmetries of the superpotential if the exotic coloured states D and D are, 
respectively, antileptoquarks and leptoquarks, with B = ^1 and L = =pl. 

All of these potential exact and approximate discrete symmetries of the 
superpotential ()3.6p are summarised in table 13.21 

It should be noted that, although the matter that survives down to low energy 
form three complete 27 representations of the broken Eq, with the exception of the 
uncharged RH neutrinos, to ensure anomaly cancellation, the imposed exact 
discrete symmetries and approximate flavour symmetries do not commute with Eq. 



46 



3.1.2 Non-Higgs supermultiplets and RH neutrinos 

It is known that in the model as presented thus far the gauge couphngs, though on 
course to unify, do not unify below the Plank scale. The beta functions above the 
SSB scale are modified compared to those of the MSSM by the existence of the 
extra matter. For example, above the SSB scale the QCD beta function is in fact 
zero at one-loop order. This issue can be solved by having the Eq GUT group be 
broken to an intermediate group before being broken finally to Gsm '^U{1)n as 
shown in ref. |61j . 

The canonical solution [55 \ \56\ \57\. however, is to introduce into the 
superpotential a bilinear term involving extra fields, known as non-Higgs fields, 
from extra incomplete 27 and 27 representations known as 27' and 27' 

W' = fi'H'H', (3.11) 

where H' is the field from 27' and H' is the corresponding field from 27'. These 
supermultiplets taken together do not spoil gauge anomaly cancellation. To some 
extent this solution reintroduces the /i problem, but /u' is not required to be related 
to the EWSB scale and in order to observe satisfactory gauge coupling unification it 
is only required that /i' < 100 TeV. The unification of the gauge couplings in the 
E^eSSM can then be achieved for any phenomenologically acceptable value of 03 at 
the EWSB scale consistent with the measured low energy central value. This is 
unlike in the MSSM where significantly higher values of 03 are required at the 
EWSB scale, well above the central measured value |57j . 

Since RH neutrinos are completely uncharged they can acquire very heavy 
Majorana masses, allowing for a type-I see-saw mechanism. Furthermore, in the 
early universe the heavy RH neutrinos, which can each decay into final states with 
lepton number either -|-1 or —1, can create a lepton asymmetry, leading to 
successful leptogenesis [62] . 
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3.2 U{1) Gauge Boson and Gaugino Mixing, EWSB 
Scale Gaugino Mass Relations, and Z-Z' Mixing 



In the low energy Lagrangian of the E'gSSM as well as the U{V) gauge boson kinetic 
terms contained in (12.21) 



\b^^^B^, - ^B'^'B'^,, (3.12) 



where the first term is the kinetic term for the U{1)y gauge boson B, with 

B^u = d^By — duBfj_, and the second term is the kinetic term for the U^l)^ gauge 

boson B', the term 



'-^B'^-^B'^^ (3.13) 



is also gauge invariant. At the GUT scale the coefficient sin(x) must be equal to 
zero since this kinetic mixing term violates the Eq gauge symmetry. Furthermore, 
this E'g-breaking mixing term is not induced by radiative corrections as long as only 
complete representations of Eq survive down to low energy. If the non-Higgs 
supermultiplets are present, however, a non-zero sin(x) can be induced [55] . 

Making the change of variables [6^ 



Bf, B^- B'^tanix)^ 

K 

B' ^ (3.14) 

^ cos(x) 

the mixing term (j3.13p is eliminated from the Lagrangian, but in the covariant 
derivative one must make the substitution 



signs' ^ (g[-^^-g,Q^ian{x)]B' = gfQ''^B', (3.15) 
V cos(x) / 
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using an effective g'^ coupling and effective cliarges 



,cfr 



and 



cos(x) 



iCff 




(3.16) 



However, even in tlie presence of non-Higgs doublets the EWSB scale relations 
= 9i = 9i ^-iid Q'^^ = are expected to be satisfied at one-loop level to within 
one- loop accuracy [55] . 

3.2.1 Soft Gaugino Masses 

In the SSB breaking part of the Lagrangian the Sg-violating soft mass term 



can also be induced at low energy, even though it is forbidden at the GUT scale. 
Along with the gauge kinetic mixing, however, this soft gaugino mass maxing is also 
expected to be small [4J. 

If at the GUT scale the soft gaugino masses Mg = Afg = = M[ = M^j^ and 
M^-^ = 0, as required by gauge invariance, then, due to the RGEs, at the EWSB 
scale one expects M{ w w V2M2 » M;^^ [55]. 

3.2.2 Z-Z' mixing 

The three VEVs Vd-, Vu, and s do not just induce diagonal masses for the Z and Z' 
bosons, but also induce a mixing term. The induced Z-Z' mass-squared matrix is 



MiiB^^^B"" + c.c. 



(3.17) 




(3.18) 
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where 



2 5^ 2 



mlz' = ~^v^Q^cos\(3)-Q^sm'ifi)], and 



= 9?vHQPcos\p) + Qrsin\l3))+gfQrs^ (3.19) 



with Q^^g the C/(l)Ar charges a down-type Higgsinos, up-type Higgsinos and 
singhnos respectively, given in table 13. 1[ The mass eigenstates are then 

Zi = Z cos{azz') + Z' sm{azz') and 

Z2 = -Z sin{azz') + Z' cos{azz'), (3.20) 



where 

1 / 2m2 \ 

Oizz' = 7: arctan — k ^ . (3.21 

2 \mi.,-mi.J 

Experimental limits on the Z2 boson mass and on the mixing angle azz' are 
model dependant since in different models that involve a Z' boson the couplings of 
that Z' will depend on the model. In the -EgSSM the most recent limit on the Z2 
boson, set by the ATLAS collaboration ^64j, searching for dilepton resonances, is 
mz2 > 1520 GeV at a confidence level of 95%. This analysis is for a Z' boson 
associated with the extra 11(1)^ of the E'eSSM, but neglects any other matter 
beyond that of the SM. When decays of the Z2 boson into inert neutralinos (inert 
Higgsino and singlino dominated mass eigenstates) are considered the Z2 width 
tends to increase by a factor of about 2 [65]. This then means that the branching 
ratio into leptons is decreased by a factor of about 2. Estimating the eff'ect of 
halving this expected branching ratio on the analysis in ref. [6l] one can read off a 
95% confidence level lower bound of around 1350 GeV. At the times of the 
publications of papers I, II, and III the most recent available limits were 
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mz2> 861 GeV [BS], m^a > 865 GeV [B7], and > 892 GeV [B5] respectively, all 
at confidence levels of 95%. Limits on the angle azz' typically require it be less 
than order 10^^ [69j. This means that neglecting azz' and setting mz^ = mz and 
mz2 = TTT-Z' ^ 9iQs^ is in most cases an excellent approximation. 

3.3 EWSB and the Active Higgs Boson Mass 
Eigenstates 

The EWSB active Higgs potential of the two active Higgs doublets Hd = Hds and 
Hu = Hu3 and the active SM-singlet S" = 5*3 is 



V{Ha,H^,S) = X'\S\'{\H/ + \HJ'\+X'\H^.H^ 



2 

+ ml\S\' + ml\Hf + ml\Hf 



+ 



XA^SH,.H+ c.c. 



+ A, (3.22) 



where mg, m^, and are the soft scalar masses for S, H^, and respectively 
and A represents the contributions from loop corrections. Once again = (t"/2 in 
SU{2)l doublet space and H^.H^ = H+ - H^H^. We define A = A333 and Ax is 
then the corresponding SSB parameter. 

Initially this EWSB sector involves ten degrees of freedom. Four of these, 
however, are massless Goldstone modes which provide the longitudinal polarisations 
of the massive W^, Zi, and Z2 bosons. When CP invariance is preserved the other 
six degrees of freedom form one charged complex scalar, one CP-odd pseudoscalar, 
and three CP-even real Higgs states. The masses of the charged and pseudoscalar 
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Higgs bosons are 



—^s--v +m^ + A± and 



(3.23) 



2 V2XAx 



sin(293) 

where A± and A^i are loop corrections and 



(3.24) 



tan((/p) = — sin(2/3). 



(3.25) 



The CP-even active Higgs sector comprises D\tH^, yicH^ and yicS. In the 
field-space basis 



h H N 



rotated by f3 with respect to the standard interaction basis such that 



DltH^a = -j=[hcos{^)-Hsm{f])+Vd], 







-j= ( h sin(/3) + H cos(/3) + w„ 



and 



V2 



(3.26) 



the mass matrix for the CP-even Higgs sector is [5] 







1 9V 




\ 






V dvdp 


dvds 






1 d^V 


1 d'^V 


1 d'^V 






V dvdf3 




V dsd(3 








1 9V 


d^v 




v 


dvds 


V dsd/3 


92 s 


/ 



Ml 



MI2 



Ml, 



Mfg Ml 



(3.27) 
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where 




11) 




(^-|).^sin(4/5) 



+ {Qu - Qd) {Qd cos\f3) + sm\f3)) sm(2/3) + A12, 





MI3 = --^^,cos(2/3) + ^(Q^-Q^)Qf z;.sm(2/3) + A23, 
Mfg = -^7;sin(2/3)+A27;s + 5f fccos2(/3) + Q^sin2(/3))Q^z;s + Ai3, and 



In p.28p Aij are the contributions from loop corrections which in the leading 
one-loop approximation are rather similar to the ones calculated in the NMSSM. 
Explicit expressions for A^ in the leading one-loop approximation are given in 
ref. [5]. Since the smallest eigenvalue of the mass-squared matrix ()3.27p is always 
less than its smallest diagonal element, at least one Higgs scalar in the CP-even 
sector, approximately h, always remains light, i.e. < ^'^11 • the leading 
two-loop approximation the mass of the lightest Higgs boson in the EgSSM does 
not exceed about 150-155 GeV. The field-space state h has couplings to SM matter 
identical to those of the SM Higgs boson for all values of tan(/3). 

When the visible sector SSB mass scale and the active SM-singlet VEV s are 
considerably larger than the EWSB scale, the mass-squared matrix ()3.27p has a 
hierarchical structure and the masses of the heaviest Higgs bosons are closely 
approximated by the diagonal entries M22 and M33 |55) . As a result the mass of 
one of the two heavier CP-even Higgs bosons, predominantly H, is approximately 
niA while the mass of the other, predominantly A^, is approximately niz'- When 
X> g'l vacuum stability requires niA to be considerably larger than mz' and the 
EWSB scale so that the qualitative pattern of the Higgs spectrum is rather similar 
to the one that arises in the PQ-symmetric NMSSM |70] . In this limit the heaviest 
CP-even, the CP-odd, and the charged states are almost degenerate with masses 




V2 
2V2 s 



sin(2/3) + m|, + A33. 



(3.28) 
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around uia [55] . 
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Chapter 4 



Thermal Relic Dark Matter 



Non-baryonic dark matter is an unknown form of matter that is believed to make 
up the majority of the matter energy density of the universe. It interacts either 
very weakly or not at all electromagnetically. The existence of dark matter was first 
proposed when analysis of orbital motion within the Coma galaxy cluster using the 
virial theorem implied that there was more mass present than just that of the 
visible baryonic matter }71j . 

Currently the evidence for the existence of cosmological dark matter is very 
strong. Its existence is inferred from galactic rotation curves [72] and from various 
measurements of galaxy clusters (see for example ref. [73]). There is also evidence 
from observations of mass inferred from gravitational lensing (see for example 
refs. [74:\ I75j). but our best measurements of the amount of cosmological dark matter 
come from fits to CMB data in the context of the standard cosmological model 
ACDM (CDM plus dark energy). Such fits to WMAP data [76] in particular give 

nsh^ = 0.0227 ±0.0006 (1-sigma) (4.1) 
for the present baryon energy density and 

nBuh"^ = 0.110 ±0.006 (1-sigma) (4.2) 
for dark, non-baryonic matter, where h w 0.73 is the reduced Hubble parameter and 
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U is the energy density divided by the critical density pc- It is thought that the 
majority of dark matter must be non-relativistic in order for the observed large 
structure formation to be explained. 

A standard assumption for pre-BBN cosmology is that the DMP was at some 
time prior to BBN in thermal and chemical equilibrium with the photon and other 
species still themselves in equilibrium with the photon. At some time in the past it 
would have then decoupled from equilibrium and under this assumption one can 
predict the relic density today of the DMP in some model if one knows all of the 
model parameters. The chemical decoupling happens roughly when the particle's 
inelastic interaction rate (maintaining chemical equilibrium) becomes less than the 
expansion rate of the universe H = a/ a. When this freeze-out occurs the number 
density of the frozen-out species typically remains much larger than it would have 
been if the species had remained in equilibrium with the photon as the universe 
cooled. If such a thermal relic particle has a freeze-out temperature that is 
much less than the mass of the particle such that the particle was non-relativistic at 
freeze-out then it is known as CDM. 



4.1 The Boltzmann Equation 



Let us assume that in some model some number of particle species, labelled with i 
in order of ascending mass rrii, are odd under some symmetry such that the 
lightest one is stable and the DMP. The evolution of the cosmological number 

density of a Z^^-odd particle species in the early universe can be expressed as 



rii 



-SHrii - ^{(7ijV-j) {n-rij - n^'^n'^) 
j 

- {(^XjiVjx){^jnx - ^r^x) 



(4.3) 
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The first term accounts for Hubble expansion and the second term accounts for 
annihilations with other Z^-odd particles, including self-annihilations. The third 
term represents the decays of Zg'-odd particles i into other Zg'-odd species j as 
well as decays of other Z^-odd species into species i. The final term represents the 
inelastic scattering of supersymmetric particles i off of -even particles X into 
other Zg'-odd species j and vice versa [771 178] . 

Summing up these equations yields the somewhat simpler expression 

-3Hn - - (4.4) 



It should be noted that, assuming that all heavier Z^^-odd particles decay into the 
DMP with not too long a lifetime, after thermal freeze-out the relic DMP number 
density will subsequently becomes equal to n. 

During thermal freeze-out the annihilation rates of the -odd particles become 
small compared to the expansion rate of the universe and their number densities 
become larger than their (non-relativistic) equilibrium values. The universe 
expands too fast for the number densities to track their equilibrium values. Let us 
assume, however, that these states inelastically scatter off of SM states X 
frequently enough that the ratios of the number densities of the -odd particles 
do maintain their equilibrium values during the time of thermal freeze-out. We 
shall refer to this as condition A and assuming that it is satisfied we have 



eq eq 

eq ^ ) 

n- n-^ n ngq 



(4.5) 



which allows us to rewrite ()4.4p as 

n = -3Hn- {av){n^ -nl^), (4.6) 
where ngq = n^'^ and the effective cross-section 

eq eq 
Tl ■ Tl ■ 



^ 3 """^ 
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Here we see that much heavier Z^-odd states, with correspondingly smaller 
non-relativistic equilibrium number densities, would be present in smaller numbers 
during the DMP's thermal freeze-out and annihilation cross-sections involving them 
would be less important. 

4.2 The Freeze-Out Temperature 

The energy density of one relativistic species of {boson, fermion} is 



Pi = 9i{hVs}-T^, (4.8) 



where Tj is the temperature of that species and gi is the number of degrees of 
freedom. We define an effective number of relativistic degrees of freedom ^efif for the 
whole system by writing 

P = 9eS^T\ (4.9) 

where p is the total density of relativistic matter and T = T-y is the photon 
temperature. The effective number of degrees of freedom ges takes into account the 
factor of 7/8 for fermions and also takes into account the fact that some species no 
longer in equilibrium with the photon may have a different temperature. The 
entropy density of a single species of {boson, fermion} with temperature Tj is 

Si = gAhVsj^Ti (4.10) 

and similarly an effective number of relativistic degrees of freedom ^eflf for the 
whole system is defined by 

s = heS^T\ (4.11) 

The numbers ^efif as ^eff will differ when any species i has a different temperature to 
the photon, with hes containing factors of {Ti/T)^ and ^eff containing factors of 
{Ti/Tf. The number density of a non-relativistic species, which the cold DMP 
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should be at freeze-out, is 



\ 3/2 / 

mu \ ( -mi 



rii = gi{^] exp(^) (4.12) 



and the energy density is simply 

Pi = mini. (4.13) 
In a radiation dominated universe the expansion rate is then given by 



H = -i^P = Jp9eS^T = hgesT , (4.14) 



where we define the constant ki for future convenience. We can approximate the 
freeze-out temperature by equating an effective DMP interaction rate with the 
radiation dominated expansion rate 



It is useful to scale the temperature by the DMP mass and define 

X = —. (4.16) 
mi 

One can then use the expression for the non-relativistic DMP number density to 
derive the transcendental equation 

" ln{^Mmi{av)P) - y2ln{x^)' ^^'"^^^ 

where 



To see how n evolves after freeze-out we first note that for isentropic expansion 
the total entropy density of the system s oc a~^, where a is the scale factor of the 
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universe. This means that s/s = 



—3H. Defining 



n 
s 



(4.19) 



we find 



^g,Mmi{av) (y2 _ y2j (4,20) 



where 



9* 



1 + ^ d/leff 
3/leff dT 



(4.21) 



This equation can be used to find the rehc density today numerically. 

By integrating from x = today to x = x^ at freeze-out one can determine the 
value of y today and the current DMP relic density is 

O = (4.22) 

Pc 

The entropy density today is dominated by the cosmic microwave and neutrino 
backgrounds. The CMB temperature is measured and the neutrino temperature 
can then be calculated as in section 17. 3[ 

The freeze-out temperature x^ depends only logarithmically on the effective 
cross-section, as in (j4.18p . but {(tv)^ is critical to determining how small y is driven 
during the time around freeze-out, before interactions become negligible. 
Subsequently, after the period of thermal freeze-out, y approximately remains 
constant as the universe expands. 
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4.3 Super symmetric Dark Matter 



In supersymmetric theories with i?-parity this iJf plays the role of and the 
DMP is the LSP |79j . In such theories the LSP is typically either the lightest 
neutralino or the gravitino, depending on the nature of the SSB mechanism which 
determines the typical scale of the gravitino mass relative to the visible sector SSB 
masses. In gravity mediated supersymmetry breaking the LSP is typically the 
lightest neutralino. 

A sub-weak-strength interacting neutralino is generally considered a good 
candidate for LSP dark matter [801 IHT] . Neutralinos do not typically form Dirac 
states and as such a neutralino DMP's relic abundance in standard cosmology is 
determined by thermal freeze-out and not by matter-antimatter asymmetry as in 
the case of baryons. 

Thermal relic neutralino dark matter has been widely studied in the 
MSSM [Ha [831 El [85] and cMSSM [861 [871 [HHl [SSI [90] . A successful dark matter 
scenario may be realised if the LSP is the lightest neutralino and there are various 
successful regions of parameter space that have different dominant annihilation 
mechanisms. For example there is the bulk region, which involves annihilation via 
t-channel slepton exchange; the focus point region, which involves annihilation via 
t-channel chargino exchange; and the funnel region, which involves annihilation via 
s-channel Higgs boson exchange. There are also regions corresponding to 
coannihilation with staus or stops. 

Typically in these scenarios ~ 1/20, meaning that the dark matter is indeed 
cold. Since, again, the freeze-out temperature is only logarithmically dependent on 
the effective cross-section, as in (j4.18p . this approximate value does not vary 
significantly for a wide variety of sub-weak-strength interacting neutralino dark 
matter scenarios, including the E'eSSM scenarios described in chapters [5l [6l and [3 
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Chapter 5 



Dark Matter in the E'qSSM 



In this chapter, which contains work that was first pubhshed in paper I, we 
present a study of neutralino dark matter in the presence of inert Higgsinos and 
singhnos, using the extended neutralino sector of the -EgSSM as an example. The 
study here should be compared to the study of dark matter in the USSM in ref. [1]. 
The particle content of the USSM, in addition to the states of the MSSM, also 
contains a SM-singlet S and a Z' boson together with their fermionic superpartners 
the singlino S and the gaugino B' . The existence of these interaction states can 
modify the nature of the neutralino LSP. In this study we include the above states 
of the USSM and also the extra inert doublet Higgsinos and singlinos predicted by 
the E'gSSM, but not included in the USSM — Hd2, Hdi, Hu2, Hui, S2, and Si. We 
do not, however, include the corresponding inert scalars which do not play a role in 
the heavy inert scalar limit. We also do not include any of the exotic coloured D 
and D states since in general we would not expect them to play a significant role in 
the calculation of the dark matter relic abundance. 

We study neutralino dark matter in the SgSSM, as defined above, both 
analytically and numerically, using micrOMEGAs [91j. We find that results for the 
relic abundance in the -EgSSM are radically different from those of both the MSSM 
and the USSM. This is because the two inert generations of doublet Higgsinos and 
singlinos predicted by the SgSSM provide an almost decoupled neutralino sector 
with a naturally light LSP that can account for the CDM relic abundance 
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somewhat independently of the rest of the model. In plausible scenarios the LSP 
annihilates predominantly through an s-channel Z boson. 

Imposing the conditions that the LSP has a mass greater than half of the Z 
boson mass, so that the LSP does not contribute at all to the Z boson invisible 
decay width, and accounts for all of the observed dark matter implies that tan(/3) 
must be less than about 2. Apart from this requirement on tan(/3), the very 
stringent constraints on MSSM or USSM parameter space that come from requiring 
that the model explains the observed dark matter relic density become completely 
relaxed since in the E'gSSM the neutralino dark matter depends almost exclusively 
on the parameters of the almost decoupled inert neutralino sector. We expect 
similar results to apply to any singlet extended supersymmetric model with an 
almost decoupled inert neutralino sector comprising extra generations of inert 
Higgsinos and singlinos. 

In section Em we discuss the inert neutralino sector of the EgSSM, introduce the 
effective model that we study, and highlight the most important couplings for our 
analysis of the LSP dark matter relic density. In section 15.21 we display the 
complete neutralino and chargino mass matrices of the E'gSSM. In section we 
present some analytical results that provide useful insights into the new inert sector 
physics. These results are subsequently used to understand and interpret the results 
of the full numerical dark matter relic density calculation using micrOMEGAs which 
are presented in section 15.41 The conclusions are summarised in section 15.51 

5.1 The Trilinear Higgs Yukawa Couplings 

The most important couplings in our analysis are the trilinear couplings between 
the three generations of down- and up-type Higgs doublets and SM-singlets 
contained in the superpotential of the E'eSSM ()3.7p 

Kjk^iHdj-Huk = \jk{SiH^-H^j^ — SiH^jH^i.). (5.1) 
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The trilinear coupling tensor Ajj^ consists of 27 numbers which play various roles. 
The purely third family coupling A = A333 is very important, because it is the 
combination fi = Xs/^/2 that plays the role of an effective // term in this theory. 
Some other neutralino mass terms, such as those involving S, are also proportional 
to A. The couplings of the inert Higgs doublets to the third generation SM-singlet 
Aq^ = Xsa/S directly contribute to neutralino and chargino mass terms for the inert 
Higgsino doublets, fdap = ^a3(3 and /u^/j = A^^a directly contribute to neutralino 
mass terms coupling an inert doublet Higgsino to an inert singlino. 

The 13 Higgs trilinear couplings mentioned thus far are the only couplings that 
obey the proposed symmetry. This approximate flavour symmetry is proposed 
in order to prevent FCNCs in the SM matter sector by eliminating non-diagonal 
flavour transitions originating from the Higgs sector. The cannot be exact as 
discussed in subsection 13.1.11 If Xijk obeyed exactly then, as we will see below, 
the neutralino mass matrix (and also the chargino mass matrix) would be 
decoupled into two independent systems and the lightest from each sector would be 
absolutely stable. We shall refer to the Z^-breaking couplings involving two third 
generation fields as x^a = ^3a3, Xua = Aasa, and Za = \a33- The notation for the 
Xijk couplings used is compiled in table 15. 1[ 



ijk 


333 


3a/3 


a3/3 q/33 


33a 3a3 a33 




A 




fdaj3 fuaf} 





Table 5.1: The abbreviated notation for the Xijk couplings. 



The 8 remaining Z2 -breaking couplings Xap-y are of less importance for our 
study. As long as only the third generation Higgs doublets and SM-singlet acquire 
VEVs then these couplings do not appear in the neutralino or chargino mass 
matrices. Additionally, they only appear in Feynman rules that involve the inert 
Higgs scalars and we assume that these are given SSB masses that are heavy enough 
such that these particles do not contribute to any processes relevant for this study. 
Similarly we neglect the exotic coloured D and D states since we expect them to be 
too heavy to play a significant role in the dark matter relic density calculation. 
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5.2 The Neutralino and Chargino Mass Matrices 



In the MSSM there are four neutralino interaction states — the neutral wino, the 
bino, and the two neutral Higgsinos. In the USSM two extra states are added - 
the singlino and the bino'. In the conventional USSM interaction basis 



Arint 
^^USSM 



B H° S & 



(5.2) 



and neglecting bino-bino' mixing, as justified in ref. [3] (see also subsection 13. 2. ip . 
the USSM neutralino mass matrix 



N 

USSM 



V 



Ml 





-\g'vd 


\g'vu 





M2 


\gvd 


-\gvu 




\gvd 





-M 


\9'vu 


-\9Vu 























Q^g'iVd 


Qug'x^u 








\ 



^ Hd gi^d 

Q^g[s 
>^5is M[ ) 



,(5.3) 



where M^^ M2, and M( are the soft gaugino masses and g are the \J{\)fq 
charges of down-type Higgsinos, up-type Higgsinos, and singlinos respectively, given 
in table [XT] In the E'gSSM this is extended. We take the full basis of neutralino 
interaction states to be 



int 



ATintT 
^^USSM 



Bd2 



Bd\ 



^1 



(5.4) 



The final six states are the extra inert doublet Higgsinos and singlinos that appear 
in the full -EgSSM. Under the assumption that only the third generation Higgs 
doublets and singlet acquire VEVs the full Majorana mass matrix is then 



M 



N 



^^-'USSM 


B2 


Bi 




^22 


^21 


Bi 


^21 





\ 



(5.5) 
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where the submatrices involving the inert interaction states are given by 



■^aB — -^Ba — 



1 

71 



K^s fu^aV sm{P) 

X/SaS fd0aV C0S{P) 

y fuaiBV sm{/3) fda^v cos{(5) 



and the Zl'^-breaking submatrices by 



Br. 



1 

71 







XdaS ZaV sin(/3) 

XuaS 2;Q'ycos(^) 

x„a^^sin(^) XdaV cos{/3) 





Similarly we take our basis of chargino interaction states to be 



a 



int 



c; 



int 



where 



Hi 

Hf 

V J 



and 



a 



(w-\ 



int 



H- 



H 



(12 



The corresponding mass matrix is then 



M 



c 



( pt\ 
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where 



M2 



\^mw sm{/3) 











P = 



V^mw cos(/3) 




(5.11) 









V 









One can already see from (j5.6p from that a typical feature of the ii^gSSM is that 
the LSP is composed mainly of inert singlino and ends up being typically very light. 
One can see this by inspecting the submatrices Aap and assuming a hierarchy of the 



s ^ v in order to satisfy the experimental limits on the Z2 boson mass. At the time 
of the publication of paper I the experimental lower limit was 861 GeV, from 
ref. [66] • The current limit is around 1350 GeV as discussed in subsection I3.2.2I 

For both the neutralinos and the charginos we see that if the Z^-breaking 
couplings are exactly zero then the inert parts of the neutralino and chargino mass 
matrices becomes decoupled from the USSM parts. However, as previously 
discussed, although approximate decoupling is expected, exact decoupling is not 
and will therefore not be considered. 

5.3 Analytical Discussion 

It will be useful to get some analytical understanding of the calculation of the dark 
matter relic abundance coming from the new neutralino/chargino physics of the 
ii^gSSM before looking at the results of the full numerical simulation. To this end, 
in this section, we consider just one inert generation consisting of two inert Higgs 
doublets and one inert SM-singlet. We label this generation as the first generation. 
We shall assume that the Z^^-breaking Yukawa couplings of the first Higgs 
generation to the third conventional Higgs generation are large enough to allow the 
neutralino/chargino states of the USSM to decay into the LSP, formed mostly from 
inert neutralino interaction states, but also small enough such that we can consider 



form Xai3S ^ f{u,d)ai3V- This is a natural assumption since we already require that 
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the inert neutralinos to be approximately decoupled from the rest of the neutralino 
mass matrix for the purposes of obtaining an analytical estimate of the mass 
eigenstates. This all amounts to considering the single block An of the extended 
neutralino mass matrix (j5.5p . 



5.3.1 The neutralino masses and mixing for one inert generation 



Within the first generation we use the basis 



N 



int 



(5.12) 



and the neutralino mass matrix is then, from ()5.6p . 



A = A 



11 



1 

72 



X's fuV sin (3 

X's fdV cos f3 

y fuV sin /3 fdv cos f3 J 



(5.13) 



where A' = An = A311, fd = fdii = A131, and /„ = fuu = A113. As discussed earlier, 
it is natural to assume that X's ^ fv and this will lead to a light, mostly first 
generation singlino lightest neutralino. 

Finding the eigenvalues of the matrix A amounts to solving a reduced cubic 
equation. Expanding in fv/X's the three neutralino masses from the first 
generation are 



1 fdfu V ■ /r,a\ , 

mi = — ^-sm(2/3) + ..., 



(5.14) 



X' s mi 



and 



(5.15) 



X' s mi 



(5.16) 



The lightest state A'^i, with mass mi, is mostly singlino (as we will confirm below) 
and the two heavier states have nearly degenerate masses, split by mi. At /3 = or 
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7r/2 the lightest neutraUno becomes massless. This is when only one of the third 
generation active Higgs doublets has a VEV. 

We shall define the neutralino mixing matrix N by 



^aj^ab^b ^ (no sum on i) (5.17) 

with superscripts indexing the interaction states and subscripts indexing the mass 
eigenstates. The lightest state is then made up of the following superposition of 
interaction states: 



(5.18) 



Again expanding in fv/X's 



Ni 



fdV la-, , 

-— cos(/3) + . 



fuV . „^ 



+ 



... J 



(5.19) 



confirming that the LSP is mostly singlino in this limit. The other eigenvectors, 
which determine the composition of neutralinos 2 and 3, are 



af + hj + 



(no sum on z), 



(5.20) 



where 



62 = 02 = — [/dCOs(/3)-/„sin(/3)]-i + --- and (5.21) 



h'i, = 03 



X's 



[/dCos(/3) + /„sin(/3)]-V 



(5.22) 
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Figure 5.1: S-channcl LSP annihilation diagrams. 



Note that a, 6 1 and that 02 and 62 Aip sign at fd cos(/3) = /„ sm{/3) whereas 03 
and 63 are always positive. Very approximately these eigenvectors are then 



No 



1 

V2 



1 

V y 



sign(/nS/3 - fdCfs) H and 



(5.23) 



1 

V2 



(A 

1 



+ 



(5.24) 



Under the assumptions of this section the chargino from the first generation is 
the first generation charged Higgsino with mass X's/^/2. 



5.3.2 Annihilation Channels 

From (j5.14p it can be seen that the LSP mass mi is proportional to v'^/s and so is 
naturally small since v <^ s. To understand this, recall that Z-Z' mixing leads to 
two mass eigenstates — Z^^Z' and Zi ~ Z — and limits on Z-Z' mixing and on 
the Zi mass place lower limits on s that imply that f <C s must be satisfied. For 
example, when s = 3000 GeV the Z2 mass is about 1100 GeV and v'^/s 20 GeV. 
The LSP mass further decreases as s becomes larger in the considered limit. In 
practice it is quite difficult to arrange the parameters such that the LSP mass 
exceeds about 100 GeV, although this depends on the sizes of Yukawa couplings 
that one is willing to accept (an issue explored more thoroughly in chapter [6]) . 
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In view of the above discussion the LSP is expected to be relatively light. When 
determining the important early universe annihilation channels we therefore begin 
by looking at s-channel annihilation, which can result in lighter mass final states. 
The most important diagrams are shown in figure 15.11 and it will turn out that the 
most important of these annihilations are those with a Z boson in the s-channel. 
The Z-N\-N\ coupling in this diagram is suppressed by a factor 

relative to the Z-neutrino-neutrino coupling under the assumptions of this section 
since the LSP only couples through its small Higgsino components. This coupling 
vanishes completely at fdCos{/3) = /uSin(/3), which is when the LSP contains a 
completely symmetric combination of H^^i and H^i- While in the MSSM a Higgsino 
dominated LSP would be expected to be such a symmetric combination of 
down-type and up-type (active) Higgsino, with mass around /i, an inert neutralino 
LSP in the SeSSM a priori has no reason to be close to such a symmetric 
combination. 

Full gauge coupling strength s-channel Z boson annihilations tend to leave a 
relic density that is too low to account for the observed amount of dark matter, but 
in this model the coupling of the mostly singlino LSP to the Z boson is typically 
suppressed, as it only couples through its doublet Higgsino admixture, leading to an 
increased relic density. As X's decreases the proportion of the LSP that is made up 
of inert doublet Higgsino, rather than inert singlino, increases. This can be seen in 
()5.19p . This then increases the strength of the overall Z-Ni-Ni coupling. The 
inclusive cross-section for s-channel annihilation through a Z-boson is therefore 
highly dependent on X's, which affects both the coupling and the LSP mass mi. 
The effect of independently increasing the coupling is always to increase the 
cross-section, but the effect of independently increasing the LSP mass can be to 
either increase or decrease the cross-section, depending on which side of the Z 
boson resonance the centre-of-mass energy is on in typical collisions during the 
period of thermal freeze-out. S-channel annihilation through the lightest Higgs 
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Z,W 





Z,W 



Figure 5.2: T-channel LSP annihilation diagrams. 



boson will also become important if typical LSP collisions are on resonance. 

The most important of the potential t-channel processes are shown in figure E21 
In practice these channels will not play a significant role compared to the s-channel 
annihilations considered previously, but we discuss them for completeness. The 
t-channel particle for these processes is one of the neutralinos or the chargino of the 
first generation. In the first diagram — t-channel annihilation to active third 
generation Higgs scalars — the couplings are approximately just / couplings of the 
first generation and appropriate mixing matrix elements. With the inert chargino 
or with inert neutralino 2 or 3 in the t-channel the diagram is approximately inert 
singlinos annihilating with an inert doublet Higgsino in the t-channel and the 
couplings are approximately just fd and /„ for producing Hd and Hu interaction 
states respectively. The LSP mass is smaller than the other masses by a factor of 
order j . With the LSP itself in the t-channel the first diagram therefore receives 
an enhancement of order jv^ for the t-channel propagator at low momentum, but 
has a suppression of order /s^ in the couplings due to the LSP only containing 
doublet type first generation Higgsinos with amplitudes of order v/s. 

The second diagram in figure [52] represents annihilation to massive gauge 
bosons. To very good approximation these bosons only couple to weak isospin 
doublets and not to SM-singlets (since Z-Z' mixing must be very small). These 
diagrams therefore have a suppression of order j relative to the full gauge 
interaction strength due to the couplings even with an inert chargino or with inert 
neutralino 2 or 3 in the t-channel. On top of this suppression these diagrams also 
receive an additional suppression of order j in the couplings, but an 
enhancement of order jv^ in the propagator when the LSP is in the t-channel. 
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This second type of diagram has a greater chance of being kinematically ahowed 
than the first. 

As previously stated, inert Higgs scalars are assumed heavy and annihilation to 
and/or through these particles is not considered. It should be noted though that 
these particles have suppressed couplings to SM matter due to the approximate 
symmetry and diagrams for the annihilation of LSPs into SM matter that involve 
these particles would be suppressed by these couplings. 



5.4 Numerical Analysis 

We now turn to the full model, in which the LSP is determined from the neutralino 
mass matrix in (j5.5p . There are two copies of the inert generation considered in the 
previous section as well as six unknown mixing parameters between the two 
generations. In general, after rotation to the mass eigenstate basis, we expect that 
two states are much lighter than the rest — both inert-singlino-like in the A's ^ fv 
limij]. 

In this section we use numerical methods to predict the relic density. We first 
diagonalise the neutralino, chargino, and Higgs scalar mass matrices numerically. 
Having done this micrOMEGAs 2.2 [9T| is then used to numerically compute the 
present day relic density, including the relevant annihilation and coannihilation 
channel cross-sections and the LSP freeze-out temperature . micrOMEGAs 
achieves this by calculating all of the relevant tree level Feynman diagrams using 
CalcHEP. The CalcHEP model files for the considered model are generated using 
LanHEP |92j . The micrOMEGAs relic density calculation assumes standard cosmology 
in which the LSP dark matter was in equilibrium with the photon at some time in 
the past, numerically solving (|4.20|) . 

^An exception to this would be the large M[ limit in which the LSP could originate from the 
lower-right block of the USSM neutralino mass matrix H5.3[l due to a mini see-saw mechanism as 
discussed in ref. [3]. 
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5.4.1 The parameter space of the model 



As justified in section \J7I\ we assume that differences between the GUT normahsed 
couphngs of the two U{1) gauge groups U{1)y and U{1)n, as weh as the mixing 
between the two groups, is neghgible, giving g[ ~ 0.46. The free parameters are 
then the trilinear Higgs couplings Xijk, the singlet VEV s, tan(/3), the soft A333 
coupling Ax, and the soft gaugino masses. It will turn out that the soft gaugino 
masses usually have little effect on the dark matter physics. One can see this by 
observing the neutralino mass matrix (|5.5p where the USSM terms coming from the 
soft gaugino masses do not directly mix with terms from the new ii^gSSM inert 
sector. The active scalar Higgs doublet and SM-singlet SSB masses are determined 
from the minimisation conditions of the scalar potential (|3.22|) given s, v, tan(/3), 
and Ax- 

In the following analysis we shall choose s = 3000 GeV and fi = 400 GeV which 
gives A = 2V2 /15 « 0.19 and makes the Z2 mass about 1100 GeV. Although much 
of the physics is highly dependent on s, this specific choice of s does not limit the 
generality of the results obtained since s always appears multiplied by a Yukawa 
coupling. This is explained in more detail below. We also choose 
M]^ = M[ = M2/2 = 250 GeV. These relations between the SSB gaugino masses are 
motivated by their RG running from the GUT scale (see subsection l3.2.ip . but the 
value is not. In this analysis the squarks and sleptons will not play a significant role 
in the calculation of dark matter relic abundance since the LSP will always be much 
lighter. We choose equal SSB sfermion masses Mg = 800 GeV and set the stop 
mixing parameter Xt, defined by 

^' = - 

where At is the SSB parameter associated with the top Yukawa coupling, to be 
equal to ^/GMg, resulting in large loop corrections to the lightest CP-even (SM-like) 
Higgs mass as in ref. [55]. This is known as the maximal mixing scenario and 
results in a lightest CP-even Higgs mass in excess of 114 GeV for all parameter 
space considered. The SSB A coupling Ax is set by choosing the pseudo-scalar 
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Higgs mass tua, from (|3.24|) . We choose itia = 500 GeV. 

We initially assume the Z|^-breaking Ajjfc couplings to be small (0.01) for the 
following analysis. The main properties of the physics can then be seen by varying 
the three parameters A' = A22 = An, / = fd22 = /«22 = fdu = fuu, and tan(/3). 
The first and second generation mixing couplings are set such that A2i,i2 = eA' and 
/(cZ,u){2i,i2) = f/- Assuming this parameter choice the sub-matrices of the neutralino 
mass matrix (15.61) become 



A 



22 



Au 



1 

V2 



^ X's /wsin(/3) 

X's /wcos(/3) 

^ fvsm{(3) fvcos{(3) 



\ 



A 



21 



eA 



22- 



and (5.26) 



(5.27) 



This simple parametrisation is sufficient for illustrating the generic properties of the 
physics. Deviations from this parametrisation are discussed afterwards. 

With the above parametrisation, the two generations are approximately 
degenerate when the mixing terms are not very large. In this case the LSP and the 
NLSP will each contain approximately equal contributions from each interaction 
basis generation. 

Finally, it is worth remarking that, assuming the above parametrisation, the 
effect on the neutralino and chargino inert sectors of changing s is simply equivalent 
to that of changing A' (although the Z2 mass will depend on s). This means that 
the following results are applicable for any experimentally consistent values of s as 
long as one accordingly scales A'. 



5.4.2 The neutralino and chargino spectra 

Figure [5^ shows how the spectrum of chargino masses varies with A'. Although the 
plot is for tan(/3) = 1.5, as one can see from (|5.1ip the inert sector of the chargino 
mass matrix has no dependence on tan(/3), with the mass terms just being 
proportional to the SM-singlet VEV. The almost constant masses are those mass 
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0.05 0.1 0.2 

A' 

Figure 5.3: Inert chargino masses (magnitude only) against A' 
with / = 1, e = 0.1, tan(/3) = 1.5, s = 3000 GeV, and -breaking 
Xijk couplings set to 0.01. 

eigenvalues coming mostly from the USSM sector — the third generation charged 
Higgsino and the wino. The charginos coming mostly from the inert sector vary 
with A' as expected and drop below the LEP lower limit around 100 GeV [93] at 
some value of A', depending on the value of s. The effect of the e = 0.1 mixing 
between generations can be seen in the splitting between the two inert sector 
charginos. Where lines cross in figure 15.31 the chargino masses are of opposite sign. 
When chargino mass lines of the same sign approach each other, they veer away 
from each other at the would-be crossing point due to interference. 

Figure [53] shows how the spectrum of neutralino masses varies with A'. The inert 
neutralino spectrum is dependent on tan(/3), but each of the qualitative features 
can be understood. We see the two light neutralino states that become heavier as A' 
decreases from unity until the approximation X's S> fv breaks down. At this point 
/f sin(/3) begins to dominate and the LSP mass decreases with decreasing A' as the 
dominance of fvsm{f3) becomes greater. In this low A' region the LSP is no longer 
mostly inert singlino, but is mostly inert up-type Higgsino. The six almost 
unvarying neutralino masses are those mostly from the USSM sector, which is not 
mixing very much with the inert sector. We have already seen that the inert sector 
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0.05 0.1 0.2 

A' 

Figure 5.4: Inert neutralino masses (magnitude only) against A' 
with / 1, e = 0.1, tan(/3) = 1.5, s = 3000 GeV, and -breaking 
Xijk couplings set to 0.01. 

chargino masses continue to be set by A' as we go down into the low A' region, 
resulting in light charginos in this region. By contrast, the four heavier inert sector 
neutralinos begin to be governed by the fv terms rather than the \'s terms in the 
low A' region and therefore approach a constant value in this region. 

As in the case of the charginos, the effect of the e = 0.1 mixing can be seen in 
the splitting between the two light neutralinos and the four heavier inert 
neutralinos which are both split by this mixing and further split by the light 
neutralino mass as predicted in the previous section. 

Figure 15.51 shows how the composition of the LSP in terms of the inert 
interaction states varies with A'. The behaviour in the X's ^ fv limit is as 
predicted in ()5.19p . We also see how the dominant component of the LSP changes 
from inert singlino to inert up-type Higgsino in the low A' region. 

5.4.3 The dark matter relic density 

Using the parametrisation in ()5.26p and ()5.27p we use micrOMEGAs 2.2 to 
numerically compute the present day relic density. Figure \5M shows a contour plot 
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0.05 0.1 0.2 
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Figure 5.5: The component structure of the LSP in terms of the 
mert interaction states against A' with f — 1, e — 0.1, tan(/3) = 1.5, 
s = 3000 GeV, and Z|^-breaking A^-fc couphngs set to 0.01. 
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Figure 5.6: Contour plot of the LSP mass and reUc density ri^/i^ 
regions in the (A', tan(/3))-plane with s = 3000 GeV, e = 0.1, and 
f — 1. The red region is where the prediction for fl^h^ is consistent 
with the measured 1-sigma range of fioivi/i^- In the region to the 
right of the hatched hue the LSP mass is less than half of the Z 
boson mass. 
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Figure 5.7: Contour plot of the LSP mass and relic density ft^h"^ 
regions in the (A',/)-plane with s — 3000 GeV, e — 0.1, and 
tan(/3) = 1.5. The red region is where the prediction for ft^h^ 
is consistent with the measured 1-sigma range of fioM^^- In the 
region to the right of the hatched line the LSP mass is less than 
half of the Z boson mass. 

of the LSP mass and predicted relic density regions in the (A', tan(/3))-plane, 

with s = 3000 GeV, e = 0.1, and / = 1. We focus on small values of A' < 0.4 since 
for large A' the LSP is a very hght, predominantly inert singlino state which does 
not annihilate very efficiently through any channel, leading to a too high relic 
density > ilDM^i^- (Such regions are shaded dark green.) As A' is decreased 

below 0.3 the LSP mass increases and approaches about half of the Z boson mass 
and there is a region where the prediction for ^^h'^ is consistent with the measured 
1-sigma range of Q^Mh'^- (Such regions are shaded red.) When the LSP mass is 
around 40 GeV it contains enough inert doublet Higgsino such that s-channel 
annihilation via the Z boson becomes strong enough to account for the observed 
relic density. As the LSP mass is increased further from 40 GeV and approaches 
45 GeV, the annihilations before freeze-out become on resonance for annihilation 
with a Z boson in the s-channel and the predicted relic density becomes lower than 
that observed. (Such regions are shaded light green.) 

In the regions where the LSP mass is less than half of the Z boson mass the LSP 
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will contribute to the effective number of neutrinos as inferred from the invisible Z 
boson width at LEP defined in ()1.13|) . The same couplings that lead to a successful 
relic density, via annihilations through an s-channel Z boson, also mean that there 
may be a significant contribution to the invisible Z boson width. For the present 
discussion it is assumed that such points, with an LSP mass lower that half of the 
Z mass, are unsafe from an experimental point of view. A detailed discussion of the 
validity of such points is postponed until chapter O Note that in the MSSM this 
issue does not arise since either the LSP is bino-like, and so does not couple to the 
Z, or is Higgsino- or Wino-like in which case it would have accompanying almost 
degenerate charginos and therefore must have a mass greater than about 100 GeV 
in any case. Here we can have an inert Higgsino/singlino LSP with a mass lower 
than half of the Z boson mass while still having experimentally consistent 
inert-doublet-Higgsino-like charginos . 

We note at this point that the requirement that the LSP mass exceeds 45 GeV 
implies low tan(/3) and this is the reason for the restricted range of tan(/3) in 
figure [5^ This can be seen from ()5.14p where we found that the LSP mass should 
be approximately proportional to sin(2/3), i.e. to the product of the two doublet 
Higgs VEVs, which is maximized at sin(2/3) = 1, corresponding to tan(/3) = 1. In 
the E'gSSM an experimentally acceptable lightest Higgs mass can be achieved even 
with tan(/3) as low as about 1.2 [55j, so having low tan(/3) is not a problem in such 
models. 

Decreasing A' further results in LSP masses above 45 GeV and to the left of the 
hatched line in figure other successful relic density regions (shaded in red) 
appear. These regions are punctuated by the light Higgs resonance, leading to the 
interesting double loop shape of the successful red regions to the left of the hatched 
line in figure ESI In these regions the LSP can have a mass significantly larger than 
half of the Z boson mass, moving far enough off the Higgs and Z boson resonances 
that annihilation is weakened just enough to give the observed relic density. 

However, another effect is observed as A' decreases. The composition of the LSP 
changes from being singlino dominated to being Higgsino dominated. For low A' the 
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cross-section begins to increase with decreasing A', even where this stih corresponds 
to increasing LSP mass, leading to a lower relic density. This is because the inert 
doublet Higgsino components in the LSP rapidly grow, as can be seen in figure [531 
At low A', when the LSP is largely inert doublet Higgsino, annihilation is too 
strong, leading to the predicted relic density being lower than that observed (as 
indicated by the light green shading in figure 15. 6p . The effects of the t-channel W 
and then Z pair production channels can also be seen as they each become relevant. 

Furthermore, for the entire successful region to the left of the hatched line in 
figure ESI the lightest chargino is heavy enough to be consistent with experiment, as 
can be seen on figure [5^ This result will be recreated for all high enough values of 
s. For larger values of s the successful regions and corresponding inert chargino 
masses are shifted down by the corresponding amount in A'. 

When X's ^ fv lowering / results in a lower LSP mass, as in ()5.14p . It also 
extends the range of A' in which this approximation is valid, i.e. it moves the 
boundary of the previously discussed low A' region to be further down in A'. 
Figure \57l\ shows the LSP mass and predicted present day relic density for different 
values of A' and / with e = 0.1 and tan(/3) = 1.5. The shifting of the successful 
region, where the LSP mass is above mz/2, down in A' at lower values of / is 
apparent. At lower values of tan(/3) this successful region extends further down in 
/. It should be noted that in order to predict the correct dark matter relic density, 
A' should be much smaller than / and that this disparity becomes greater if s is 
increased. Increasing s effectively just shifts all of the features on figures [531 and [5771 
to the left. 

5.4.4 Deviations from the considered parametrisation 

Breaking the relation fu{22,ii) = /d(22,ii) can have similar effects to those of 
changing tan(/3). However, because these parameters cannot be too high (in order 
for the theory to be perturbative up to the GUT scale) and because lowering them 
to much less than unity makes the LSP too light, tan(/3) can be varied much more 
freely than the /„/ fd ratio. 
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The effect of increasing the inert generation mixing parameter e is to increase 
the various mass splittings between similar inert mass eigenstates. Increasing the 
mixing between the first and second generations thus results in a lighter LSP, 
shrinking the successful region, and a lighter lightest chargino, potentially 
inconsistent with current chargino non-observation. 

The physics of the inert sector when deviating from the currently considered 
parametrisation is studied much more carefully in chapter [6l 

As long as the LSP is still mostly from the inert sector, as considered here, other 
parameters do not greatly affect the dark matter physics and are effectively free. 
Squark and slepton parameters do not affect the dark matter physics of the 
considered model. Top and stop loops can have a significant effect on the lightest 
Higgs mass, but as long as this mass is experimentally allowed then these 
parameters are also not constrained by the requirement that the model produce 
successful dark matter. 

5.5 Summary and Conclusions 

In this work we studied inert neutralino dark matter arising in supersymmetric 
models with extra inert Higgsinos and singlinos. As an example we considered the 
extended neutralino sector of the SgSSM. This work represents a first study of the 
inert neutralino sector of the SeSSM and it is found that in this model the LSP 
does typically arise from this sector. We studied this novel dark matter scenario 
both analytically and numerically, using micrOMEGAs. 

The dark matter scenario differs greatly from those of the MSSM and USSM 
since the two inert neutralino generations provide an almost decoupled neutralino 
sector with a naturally light LSP that can account for the CDM relic abundance 
somewhat independently of the rest of the model. Although the SgSSM has two 
inert generations, the presence of the second inert generation is not crucial to the 
dark matter scenario. 

In the successful regions where the observed dark matter relic density is 
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reproduced the neutralino mass spectrum is well described by the analytical results 
of section (531 In this region the LSP is mostly inert singlino and has a mass 
approximately proportional to f^/s, as in (j5.14p . and as X's is decreased the LSP 
becomes heavier and also less inert singlino dominated, picking up significant inert 
doublet Higgsino contributions. 

To avoid potential confiict with high precision LEP data we considered the case 
where the LSP mass is above half of the Z boson mass. Since the LSP mass in 
(j5.14p is proportional to sin(2/3), we found that such regions of parameter 
space in which the dark matter relic density prediction is consistent with 
observation require low values of tan(/3) — less than about 2. Depending on the 
value of the singlet VEV s, the f[u,d)ai3 trilinear Higgs coupling parameters should 
also be reasonably large compared to the Aq,^ ones. 

One of the main messages arising from this work is that neutralino dark matter 
could arise from an almost decoupled sector of inert Higgsinos and singlinos and 
that if it does then the parameter space of the rest of the model is completely 
opened up. For example, if such a model is regarded as an extension of the MSSM 
then the lightest MSSM-like supersymmetric particle is not even required to be a 
neutralino and could even be a sfermion. In the E'gSSM the lightest MSSM-like 
supersymmetric particle can decay into the inert LSP via Z^^-breaking Ajj^ 
couplings that need not be extremely small. 

Similar results should apply to any singlet extended SSM with one or more 
extra, inert generations of Higgsinos and singlinos with a trilinear Higgs coupling 
tensor equivalent to that in ()5.ip . 
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Chapter 6 



Novel Higgs Decays in the 

The discovery of the Higgs boson, the last remaining undiscovered particle of the 
SM, is one of the main goals of the LHC. The strategy for Higgs searches depends 
on the decay branching fractions of the Higgs boson into different channels. Physics 
beyond the SM may affect the Higgs decay rates to SM particles and give rise to 
new Higgs decay channels necessitating a drastic change in the strategy for Higgs 
boson searches. (For recent reviews of non-standard Higgs boson decays see 
refs. \94: \ 195 1 [96].) I^i particular there exist several extensions to the SM in which 
the Higgs boson can decay with a substantial branching fraction into particles that 
cannot be directly detected. 

These invisible Higgs boson decays can occur in supersymetry, with the lighest 
Higgs boson decaying into neutralino LSPs. In some regions of MSSM parameter 
space the lightest Higgs boson decays into the lightest neutralino with a relatively 
large branching ratio, giving rise to invisible final states if i?-parity is 
conserved [97]. LEP and Tevatron data allow the neutralino LSP to be sufficiently 
light such that the decays of the lightest Higgs boson into these neutralinos is 
kinematically allowed and such light neutralinos can annihilate efficiently through a 
Z boson pole resulting in a reasonable density of dark matter. 
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The presence of invisible decays considerably modifies Higgs boson searches and 
makes discovery much more difficult. If the Higgs boson is mainly invisible then the 
usual visible branching ratios will be dramatically reduced, preventing detection in 
the much studied channels at the LHC and Tevatron. In the case where invisible 
Higgs boson decays dominate it is impossible to fully reconstruct a resonance and it 
is very challenging to identify the Higgs boson at collider experiments, i.e. the 
quantum numbers remain unknown. At e^e^ colliders, the problems relating to the 
observation of an invisible Higgs boson are less severe \98\ [991 llOOj since it can be 
tagged through the recoiling Z. The LEP exclusion of Higgs boson masses up to 
114.4 GeV applies even in the case of invisibly decaying Higgs bosons [lOlj and 
similar limits could apply to Higgs bosons decaying into soft lepton pairs some 
fraction of the time, as happens for some of the novel Higgs decay scenarios 
discussed in this chapter. 

Higgs boson searches at hadron colliders, however, are more difficult in the 
presence of such invisible decays. Previous studies have analysed Zh and Wh 
associated production |102[ I103[ I104j as well as tth production jl05[ I106j and tiVV 
and bbVV production [107j as promising channels, where h is the Higgs boson and 
the Vs stand for vector bosons. The possibility of observing an invisible Higgs 
boson in central exclusive diffractive production at the LHC was studied in |108j . 
Another proposal is to observe such an invisible Higgs in inelastic events with large 
missing transverse energy and two high Et jets. In this case the Higgs boson is 
produced by VV fusion and has a large transverse momentum resulting in a signal 
with two quark jets with distinctive kinematic distributions compared to Zjj and 
Wjj backgrounds [TIM fTM fTTn] . 

In this chapter, which contains work that was first published in paper II, we 
consider novel decays of the lightest Higgs boson and associated collider signatures 
within the E'eSSM. If the Yukawa couplings of the inert neutralino sector are 
required to be small enough such that perturbation theory remains valid up to the 
GUT scale then the masses of the two lightest inert neutralino states are expected 
to be smaller than about 60-65 GeV. As a result the lightest inert neutralino tends 
to be the LSP. As noted in the previous chapter such an inert neutralino can give 
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an appropriate contribution to the dark matter density, consistent with recent 
observations, if it has a mass around 35-50 GeV. In this case the hghtest Higgs 
boson decays predominantly into inert neutrahno states and the usual Higgs boson 
branching ratios to SM particles are less than a few percent. 

In section 16.11 we look in more detail at the inert sector of the model and the 
couplings of the inert neutralinos, inert charginos, and active Higgs bosons are 
specified. Novel decays of the lightest CP-even Higgs state and dark matter 
constraints are discussed in section 16. 2i In section 16.31 we specify some benchmark 
points and discuss the experimental constraints and predictions. The conclusions 
are summarised in section 16.41 

6.1 Inert Charginos and Neutralinos 

In our analysis we will assume that Z^-violating couplings are small and can be 
neglected. This assumption can be justified if one takes into account that the 
Z^-violating operators can give an appreciable contribution to the amplitude of 
K^-K^ oscillations and give rise to new muon decay channels such as — >• e~e^e~ . 
In order to suppress processes with non-diagonal flavour transitions the Yukawa 
couplings of the exotic particles to the quarks and leptons of the first two 
generations should be smaller than 10~'^-10~^. Such small Zg^-violating couplings 
can be neglected in the first approximation. 

In this approximation, and given the assumption that only Hd, and S 
acquire non-zero VEVs, the charged components of the inert Higgsinos do not mix 
with the MSSM-like chargino states. The neutral components of the inert Higgsinos 
and inert singlinos also do not mix with the USSM-like neutralino states. If 
symmetry was exact then both the lightest state in the ordinary neutralino sector 
and the lightest inert neutralino would be absolutely stable. Therefore although 
Z^-violating couplings are expected to be rather small we shall assume that they 
are large enough to allow either the lightest USSM-like neutralino state or the 
lightest inert neutralino state to decay within a reasonable time — the lighter of the 
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two being the stable LSP and the dark matter candidate. 
In the basis 



AT.™* , = { ffO frO d frO ffO d 
-"inert I -n^r, 11. ,o On -Oji 



inert ^ -fJ^g J^u2 '-'2 ^dl ^ul '-'l 

the inert part of the neutrahno mass matrix is given by 

/ 



-^inert 



V 



^22 ^21 



4T A 

^21 ^11 



(6.2) 



with the submatrices given in (|5.6|) . The inert part of the chargino mass matrix P, 
given in (15. lip , may be written 

From (|6.2p and (|6.3p one can see that in the exact symmetry Hmit the 
spectrum of the inert neutrahnos and charginos in the -EgSSM can be parametrised 
in terms of Xa)3, fdap-, fuap-, tan(/3), and s. In other words the masses and couplings 
of the inert neutrahnos are determined by 12 Yukawa couplings, which can in 
principle be complex, tan(/3), and s. Four of the Yukawa couplings mentioned 
above — A^/j — as well as the VEV of the SM singlet field s set the masses and 
couplings of the inert chargino states. Six off-diagonal Yukawa couplings define the 
mixing between the two generations of inert Higgsinos and singlinos. 

In the following analysis the VEV of the active SM-singlet field is chosen to be 
large enough (s > 2400 GeV) so that experimental constraints from ref. [67] on the 
Z2 boson mass {m^^ > 892 GeV) and Z-Z' mixing are satisfied. Since the 
publication of paper II the limit on the Z2 mass in the E'eSSM has increased as 
discussed in subsection 13.2.21 In order to avoid the LEP lower limit on the masses 
of inert charginos [93j the Yukawa couplings Aq,^ are chosen such that all inert 
chargino states are heavier than 100 GeV. In addition, we also require the validity 
of perturbation theory up to the GUT scale and this constrains the allowed range of 
all Yukawa couplings. 
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The theoretical and experimental restrictions specified above set very strong 
limits on the masses and couplings of the lightest inert neutralinos. In particular, 
our numerical analysis indicates that the lightest and second lightest inert 
neutralinos are always light. They typically have masses below 60-65 GeV. These 
neutralinos are predominantly inert singlino in nature. From our numerical analysis 
it follows that the lightest and second lightest inert neutralinos might have rather 
small couplings to the Z boson so that any possible signal that these neutralinos 
could give rise to at LEP would be extremely suppressed. As a consequence such 
inert neutralinos would remain undetected. At the same time four other inert 
neutralinos, which are approximately linear superpositions of neutral components of 
inert doublet Higgsinos, are normally heavier than 100 GeV. 



6.1.1 The diagonal inert Yukawa coupling approximation 

In order to clarify the results of our numerical analysis it is useful to consider a few 
simple cases that give some analytical understanding of our calculations. The 
simplest case is when all of the Yukawa coupling from the off-diagonal blocks of 
(|6.2p are zero such that 

fda^ai3, and 

faSa/B (no sum on a). (6.4) 

This leads to two decoupled generations with the properties studied in section 15.31 
The mass matrix of inert neutralinos (j6.2p reduces to block diagonal form while the 
masses of the inert charginos are given by 

^c^ = (6-5) 

When fa = fda = fua one can prove using the method proposed in ref. [lllj that 
there are theoretical upper bounds on the masses of the lightest and second lightest 



^ajS — 
fdaP = 
fuaB 
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inert neutralino states. The theoretical restrictions are P ^ /^^e' where 



.2 _ 1 



2 



(6.6) 



The value of decreases with increasing and tan(/3), approaching its 

maximum value 

Ha (6.7) 

as Tn^^ — >• and tan(/3) — >• 1. 

The upper bound on the mass of the lightest inert neutralino also depends on 
the values of the Yukawa couplings fda and fua ■ The theoretical restrictions on 
these couplings due to the requirement that the theory should remain perturbative 
up to the GUT scale become weaker with increasing tan(/3). At large values of 
tan(/3) the upper bounds on \ from (|6.6p becomes rather small and as tan(/3) 
tends to unity the upper bounds on \m^,J again become rasther small, because 
theoretical constraints on /^q, and f^a become rather stringent. Taking both of 
these effects in account the upper bounds on \m^ \ achieve their maximum values 
around tan(/3) ~ 1.5. For this value of tan(/9) the requirement of the validity of 
perturbation theory up to the GUT scale implies that for / = f^n = fui = fd2 = fu2 
f must be less than about 0.6. As a consequence the lightest inert neutralinos are 
lighter than around 60-65 GeV for |?tI(^^| > 100 GeV. 

Using the results from section 15.31 for the compositions of the light neutralinos 
from each inert generation one can derive the couplings of these states to the Z 
boson. We define Rza/s couplings such that the Z-Na-Np coupling is equal to Rzaji 
times the Z-v-v coupling 

Rzc.fi = N^N} - N^Nj + N^N^ - N^Nl (6.8) 
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where is the neutrahno mixmg matrix element corresponding to mass eigenstate 
i and inert interaction state a in the basis (|6.1|1 . 

In the case where the off-diagonal inert Yukawa coupling blocks vanish while 
^ f(u,d)aV the relative couplings of the lightest and second lightest inert 
neutralino states to the Z boson are given by 

Rzap = RzaaSa(3 (nosumona), (6.9) 

where 

Rzaa = ^(/icos2(/3)-/2„sin2(/3)). (6.10) 

This demonstrates that the couplings of A^i and to the Z boson can be very 
strongly suppressed. It becomes zero when \ fda\ cos(/3) = \ fua\ sin(/3), which is when 
Na contains a completely symmetric combination of H^^ and H^a- (|6.10p also 
indicates that the couplings of A'^i and to the Z boson are always small if the 
inert charginos are rather heavy or if f^a and fua are small, i.e. the masses of A'^i 
and ai^e small. 

6.1.2 A27 and pseudo-Dirac lightest inert neutralino states 

In order to provide an explanation for the origin of the approximate symmetry 
that singles out the third generation of Higgs doublets and SM-singlets, and to 
account for tribimaximal mixing and other features of the quark and lepton spectra, 
a A27 flavour symmetry may be applied to the E'gSSM [59j. The addition of the 
A27 flavour symmetry implies an inert neutralino mass matrix with An ~ A22 ~ 0, 
leading to approximately degenerate lightest neutralino states with a pseudo-Dirac 
(see appendix |A]) structure. 

When all flavour diagonal Yukawa couplings Xaa, fdaa-, and fuaa exactly vanish, 
i.e. All = ^22 = 0, the inert neutralinos form Dirac states. In this limit the 
Lagrangian of the E'gSSM is invariant under an extra U{1) global symmetry. Under 



91 



this symmetry the fermionic components of the inert supermultiplets transform 



^2 - 


■> e b2, 


Hd2 - 


, „ici fj 

e Hd2, 


Hu2 - 


■> e nu2, 


Si - 




Hdi - 




Hul - 





In the above Hmiting case the hghtest inert neutraUno is a Dirac state formed 
predominantly from Si and 5*2 . In this case the LSP and its antiparticle have 
opposite charges with respect to the extra global U{1) and this could lead to the 
scenario known as asymmetric dark matter |1121 1113t I114^ I115j . The ADM scenario 
supposes that there could be an asymmetry between the density of dark matter 
particles and their antiparticles in the early universe similar to that for baryons. 
This could have a considerable effect on the relic density calculations. In particular, 
if an asymmetry exists between the number densities of dark matter particles and 
their antiparticles in the early universe then one can get an appreciable dark matter 
density even if the dark matter particle-antiparticle annihilation cross section is 
very large, like in the case of baryons. Furthermore, if most of the dark matter 
antiparticles are eliminated by annihilation with their particles then such an ADM 
scenario does not have the usual indirect signatures associated with the presence of 
dark matter. For example, there would be no high energy neutrino signal from 
annihilations in the Sun. At the same time, a relatively high concentration of dark 
matter particles can build up in the Sun, altering heat transport in the solar 
interior and affecting low energy neutrino fluxes [115j . 

In practice the A27 scenario tells us that we are somewhat away from the above 
limiting case, with a broken global U{1) symmetry leading to almost degenerate 
pseudo-Dirac lightest neutralinos, where the relic density of the LSP can be 
calculated by standard methods. It will turn out that the LSP cannot be too light 
(must be of order mz/'2) in order not to have a too high cosmological relic density. 
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At the same time we will see that the two lightest neutralinos cannot be too heavy 
in order for perturbation theory to be valid up to the GUT scale. In practice this 
means that in realistic scenarios the two lightest inert neutralino states are rather 
close in mass. The A27 scenario provides an explanation for this feature of the 
successful neutralino mass pattern. 

It is worth noting that the results from the previous section can be reinterpreted 
in terms of this scenario. Specifically in the case where = = and 
A21 = A12 a block diagonalisation of the inert neutralino mass matrix (j6.2p results 
in 

A22 A'22 = -A21, and All ^ = ^21, (6-12) 

with A21 = A12 — )• A'21 = A!y2 = 0. This only corresponds to a redefinition of the 
generations 1 and 2 and does not mix fields of different hypercharge. This provides 
a dictionary between these two scenarios 



~ -^22 = 


= A'n 


= A21, 


~/d22 = 


/dll 


= fd2l 


~fu22 — 


full 


= fu21 



(6.13) 
(6.14) 



Rewriting the inert neutralino mass matrix in this block diagonal form also makes 
it clear that the Rzvi coupling vanishes in this limit in the same way that it did for 
the diagonal case in subsection 16. 1.1[ 



6.1.3 The couplings of Higgs bosons to inert neutralinos 

The presence of light inert neutralinos in the particle spectrum of the E'gSSM 
makes possible the decays of the Higgs bosons into these final states. Now and in 
the next section we argue that such decays may result in the modification of the 
SM-like Higgs signal at current and future colliders. Since our main concern in this 
work is the decays of the SM-like lightest Higgs boson, we shall ignore the effects of 
the inert Higgs scalars and pseudoscalars which do not mix appreciably with the 
active scalar sector responsible for EWSB. We also assume that all of the inert 
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bosons are heavier than the hghtest CP-even Higgs boson. 

If all other Higgs boson states are much heavier than the lightest CP-even Higgs 
boson then the lightest Higgs state, approximately given by h, as defined in (|3.26|) . 
manifests itself in interactions with SM gauge bosons and fermions as a SM-like 
Higgs boson. Since within the E'gSSM the mass of this state is predicted to be 
relatively low, its production cross section at the LHC should be large enough so 
that it can be observed in the near future. In this context it is particularly 
interesting and important to analyse the decay modes of the lightest CP-even Higgs 
state. Furthermore, we concentrate on the decays of the SM-like Higgs boson into 
the lightest and second lightest inert neutralinos. 

The couplings of the Higgs states to the inert neutralinos originate from the 
interactions of Hu, H^, and S with the inert Higgs fields in the superpotential. 
Using ()3.26p one can express VizH^, d\e H^, and Die S in terms of the field-space 
basis states h, H, and N. The components of the field-space basis are related to the 
physical CP-even mass eigenstates by a unitary transformation 



U 



H 



(6.15) 



Combining all these expressions together one obtains an effective Lagrangian term 
that describes the interactions of the inert neutralinos with the CP-even Higgs 
mass eigenstates 

Xl'fhmNfN^+c.c., (6.16) 

where 



X 



fey il-m tj 



V2 



[/;^^cos(/3)-C/f^sin(/3))P, 



dij 



-1= (C/,^^ sin(/3) + cos(/3)) P,,,, 



.17) 
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with 



A,,. = \^^NfN^ + X^^NfN] + \^^N}N^ + \^^N}Nl (6.18) 
Fdrj = fdiiNfN^ + fdi^NfN] + f,^^NfN^ + f,^^NfNl and (6.19) 
Fuij = fuuN^Nf + f,,,N!Nf + f,,,NfNf + f^,,NfNf. (6.20) 

The expressions for the couphngs of the active CP-even Higgs scalars to the 
inert neutrahnos become much simpler in the case where the Higgs spectrum has 
the usual hierarchical structure. In this case U is almost the identity. As a 
consequence the couplings of the SM-like Higgs boson to the lightest and second 
lightest inert neutralino states are approximately given by 

Xj^l = --L(^i7^^^cos(/3)+F„„^sin(/3)). (6.21) 

In the case of the diagonal inert Yukawa coupling limit defined in subsection 16.1.11 
and if AqS » f(d,u)aV one can use the expression (|5.19p to find values for Nf and 
N2 and derive approximate formulae for X^^. Substituting into (|6.21|) one obtains 

Xll = -^5ai3 + --- (no sum on a). (6.22) 

This simple analytical expression for the couplings of the SM-like Higgs boson to 
the lightest and second lightest inert neutralinos is not as surprising as it may first 
appear. When the Higgs spectrum is hierarchical, with s ^ v, the VEV of the 
lightest CP-even state v is responsible for all light fermion masses in the i^gSSM. 
As a result we expect that their couplings to the SM-like Higgs can be written as 
usual as being proportional to the mass divided by the VEV. We see that this is 
exactly what is found in the limit of being small. 
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6.2 Novel Higgs Decays and Dark Matter 



The interaction Lagrangian (|6.16p gives rise to decays of the hghtest Higgs boson 
into inert neutrahno pairs with partial widths given by 



No, 



m 



m 



mt 



mt 



1 2 

m - m- 

4 



3.23) 



where A^^ = {1,2} for {a = /?, a 7^ /?}. 

The partial widths associated with these inert decays of the SM-like Higgs boson 
()6.23p have to be compared to decay rates into SM particles. When the SM Higgs 
boson is relatively light (less than about 140 GeV) it decays predominantly into h 
quark and r lepton pairs. The partial decay width of the lightest CP-even Higgs 
boson into Dirac fermion pairs is given by |116j 



r(/ii ^ //) 
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3.24) 



For the case of the decays into r leptons the coupling of the lightest CP-even Higgs 
state to the r lepton normalised to the corresponding SM coupling 



1 



cos(/3) 



C/,\cos(/3)-C/,^sin(/3) 



3.25) 



For a final state that involves h quarks one has to include the QCD corrections. 
In particular, the fermion mass in ()6.24p should be associated with the running h 
quark mass fhii{p). The bulk of the QCD corrections are absorbed by using the 
running b quark mass defined at the appropriate renormalisation scale — the scale 
of the lightest Higgs boson mass /i = m/i^ in the considered case. In addition to the 
corrections that are associated with the running b quark mass there are other QCD 
corrections to the Higgs coupling to b quarks that should be taken into 
account |117j . As a consequence, the partial decay width of the lightest CP-even 
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Higgs boson into b quark pairs can be calculated using (|6.24|) if one sets 



Nc = 3, 

nif = fhi,{mh^), and (6.26) 



2 



where 



^ 5.67^^ + (35.94 - 1.36iV,)%^ and 



(6.27) 



Here we neglect radiative corrections that originate from loop diagrams that 
contain non- SM particlej]. 

From (|6.22|) one can see that in the SeSSM the branching ratios of the SM-like 
Higgs state into the lightest and second lightest inert neutralinos depend rather 
strongly on the masses of these particles. When the lightest inert neutralino states 
are heavy relative to the b quark the lightest Higgs boson decays predominantly 
into NaNp while the branching ratios for decays into SM particles are suppressed. 
On the other hand if the lightest inert neutralinos have masses that are 
considerably smaller than the masses of the b quark and r lepton then the 
branching ratios of the SM-like Higgs into inert neutralino final states are small. 

Constraints on the mass of the lightest inert neutralino can be obtained if we 
require that this particle accounts for all or some of the observed dark matter relic 
density. In the limit where all non-SM fields other than the two lightest inert 
neutralinos are heavy the lightest inert neutralino state in the ii^gSSM is responsible 
for too large a thermal relic density of dark matter. The LSP A^i is composed 
mainly of inert singlino and has a mass inversely proportional to the charged inert 
Higgsino mass. In this limit it is typically very light with \ <^ mz- As a result 



^Radiative corrections that are induced by supersymmetric particles can be very important, par- 
ticularly in the case of the bottom quark at high values of tan(/3). For a review see ref. [118| . 
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the couplings of the lightest inert neutralino to gauge bosons, the SM-like Higgs 
boson, quarks, and leptons are quite small, leading to a relatively small dark matter 
annihilation cross-section into SM particles and giving rise to a relic density that is 
typically much larger than the measured value. Thus, in the limit considered, the 
bulk of the -EgSSM parameter space that leads to small inert neutralino masses is 
ruled out. 

The situation changes dramatically when the mass of the lightest inert 
neutralino increases. In this case the Higgsino components of A^i become larger and 
as a consequence the couplings of A'^i to the Z boson grow. A reasonable density of 
dark matter can be obtained for | ~ Tnz/2 when the lightest inert neutralino 
states annihilate mainly through an s-channel Z boson. It is worth noting that if 
A^i were pure inert Higgsino then the s-channel Z boson annihilation would proceed 
with the full gauge coupling strength leaving a relic density too low to account for 
the observed dark matter. In the i^gSSM the LSP is mostly inert singlino so that 
its coupling to the Z boson is typically suppressed, since it only couples through its 
inert Higgsino admixture, leading to an increased relic density. In practice an 
appropriate value of J^dm^^^ can be achieved even if the coupling of Ni to the Z 
boson is relatively small. This happens when A'^i annihilation proceeds through the 
Z boson resonance. Thus, scenarios that result in a reasonable inert neutralino dark 
matter relic density correspond to lightest inert neutralino masses that are much 
larger than ffn,{mh-^) and hence to the SM-like Higgs boson having very small 
branching ratios into SM particles. 



6.3 Benchmarks, Constraints, and Predictions 

In order to illustrate the features of the SgSSM mentioned in the previous section 
we specify the set of benchmark points in tables 16.11 16. 2| and 16.31 For each 
benchmark scenario we calculate the spectra of inert neutralinos, inert charginos, 
and active Higgs bosons as well as their couplings, the decay branching ratios of the 
lightest CP-even Higgs state, and the dark matter relic density. micrOMEGAs 2.2 is 
used to numerically compute the present day density of dark matter. 
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Table 6.1: Benchmark scenarios for 133-135 GeV. The 

branching ratios and decay widths of the hghtest Higgs boson; 
the masses of the active Higgs bosons, inert neutraUnos, and 
charginos; and the couphngs of the inert neutrahnos A^i and N2 
are calculated for s ^ 2400 GeV, A = 0.6, Ax = 1600 GeV, 
niQ — niu = Ms = 700 GeV, and Xt — VQMs, corresponding 
to TOft, w « 890 GeV. AV^fP and D are defined in 
and (|6.30p respectively. 
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Table 6.2: Benchmark scenarios for 771/,,^ « 114-116 GeV. The 
branching ratios and decay widths of the hghtest Higgs boson; the 
masses of the active Higgs bosons, inert neutrahnos, and charginos; 
and the couphngs of the inert neutralinos iVi and sue cal- 
culated for s = 2400 GeV, X ^ g[ ^ 0.468, Ax = 600 GeV, 
TOQ = TOti — Ms = 700 GeV, and Xt = \/6Ms, corresponding 
to mh,_ K mz^ ~ 890 GeV. AV^ep and D are defined in (lOS)) 
and (I6.30p respectively. Continued in table 16.31 
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Table 6.3: Continued from table |621 more benchmark scenarios 
for m/ij « 114-116 GeV. Again, the branching ratios and decay 
widths of the hghtest Higgs boson; the masses of the active Higgs 
bosons, inert neutrahnos, and charginos; and the couphngs of the 
inert neutrahnos A^i and are calculated for s — 2400 GeV, 
\ = g{^ 0.468, Ax = 600 GeV, mg = m„ = M, = 700 GeV, and 
Xt = corresponding to m,,, w mz^ « 890 GeV. 

and D are defined in (|6.28p and (|6.30p respectively. 
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6.3.1 Benchmark scenarios 

In order to construct benchmark scenarios that are consistent with cosmological 
observations and colhder constraints we restrict our considerations to tan(/3) < 2. 
The plots in figure [6TT] show that in principle an appropriate value of the dark 
matter density can be obtained even when tan(/3) > 2. However, larger values of 
tan(/3) lead to the lightest and second lightest inert neutralinos having smaller 
masses as discussed in section 15.31 As a result larger couplings of the lightest inert 
neutralinos to the Z boson are required to reproduce the measured value of Vl^yifx^ 
and such light inert neutralinos with substantial couplings to Z boson give a 
considerable contribution to its invisible width leading to a conflict with LEP 
measurements. This is discussed in more detail in the following subsection. 

Even for tan(/3) < 2 the lightest inert neutralino states can get appreciable 
masses only if at least one of the inert chargino mass eigenstates is light 
tUq^ ~ 100-200 GeV. As clarified in sections 16.11 and 15.31 the masses of the lightest 
inert neutralino states decrease with increasing and it is therefore rather 

" (-■1,2 

difficult to find benchmark scenarios consistent with cosmological observations for 
vn^^ > 200 GeV. At the same time we demonstrate (with benchmark ix in 
table 16. 3p that one light inert chargino mass eigenstate is enough to ensure that the 
lightest inert neutralino state gains a mass of order m^/2. 

To obtain the kind of inert neutralino and chargino spectra discussed above one 
has to assume that the couplings A^^ are rather small. They are expected to be 
much smaller than the largest ^daji and ^ua^ couplings. On the other hand, in order 
to get |?Ti^J ~ l^i-TVal ~ f^z/'^ the Yukawa couplings fdap and fuap need to be 
relatively close to their theoretical upper bounds caused by the requirement of the 
validity of perturbation theory up to the GUT scale. Since gauge coupling 
unification and RG fiow determine the low energy value of g'^ , the mass of the Z2 
gauge boson is approximately set by the SM-singlet VEV s only. In our study we 
choose s = 2400 GeV so that the Z2 mass is about 890 GeV. This value of the Z2 
boson mass is just above the lower bound of 865 GeV found in ref. [67j — the most 
recent limit at the time of the publication of paper II — and allows satisfaction of 
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Figure 6.1: Contour plots of (Xf/)^ and various regions in the 
(/,tan(/3))-plane with s = 2400 GeV, fdaa = fuaa = Ka = Va, 
fd2i = /, fu2i = fd2i/a, fdi2 = 1.02/<i2i, /«i2 = 0.98/„2i, and 
A21 = A12 = 0.06, implying that to^,^ ^ = 101.8 GeV. The upper 
plot is for a — 0.75 + 0.25tan(/3) and in the lower plot is for a = 
0.5 + 0.5tan(/3). The red region is where the prediction for fl^h^ 
is consistent with the measured 1 -sigma range of ^jyM^ given in 
(|4.2p . The dark green region corresponds to I? < 3 while the 
pale green region represents the part of the parameter space in 
which D is between 3 and 4. The grey area indicates that D > 4. 
D is defined in (j6.30D . The blue region corresponds to TOj^ > 
mz/'2 while the dark blue region to the right is ruled out by the 
requirement that perturbation theory remains valid up to the GUT 
scale. 
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stringent limits on the mass and Z-Z' mixing that come from precision 
electroweak tests [69] . 

Since we restrict our analysis to low values of tan(/3) < 2 the mass of the SM-like 
Higgs boson is very sensitive to the choice of the coupling A. Stringent LEP 
constraints require A at the EWSB scale to be larger than the low energy value of 
g'^ ~ 0.47 and if one increases A much further then the theoretical upper bounds on 
fdap and fuoL^ from RG running become substantially stronger. As a consequence, 
it is rather difficult to find solutions with jw-^^l ~ l'^7V2l ^ f^z/'^- Therefore in our 
analysis we concentrate on values of A at the EWSB scale less than about 0.6. In 
addition we set stop SSB masses mq = rriu = Mg = 700 GeV and restrict our 
consideration to the maximal mixing scenario with the stop mixing parameter, 
defined in (|5.25|) , Xt = ^/^Mg . This choice of parameters limits the range of 
variations of the lightest CP-even Higgs mass. In the leading two- loop 
approximation the mass of the SM-like Higgs boson varies from 115 GeV for \ = g'^ 
to 136 GeV for A = 0.6. From tables 16. 2| and 16.31 one can see that the large 
values of A > 5^ that we choose in our analysis result in an extremely hierarchical 
Higgs spectrum, as pointed out in section [3^31 In tables W\\ 16. 2| and 16.31 the masses 
of the heavier Higgs states are computed in the leading one-loop approximation. In 
the case of the lightest Higgs boson mass the leading two-loop corrections are taken 
into account. 

The set of benchmark points that we specify demonstrates that one can get a 
reasonable dark matter density consistent with recent observations if 
|m^^| ~ l^-^jl ~ Triz/"^- Our benchmark scenarios also indicate that in this case the 
SM-like Higgs boson decays predominantly into the lightest inert neutralinos A'^i 
and while the total branching ratio into SM particles varies from 2% to 4%. 

Benchmarks i, ii, iv, v, and viii are motivated by the non-Abelian flavour 
symmetry A27 which describes well the observed hierarchy in the quark and lepton 
sectors. As was discussed in subsection I6.1.2t these scenarios imply that all flavour 
diagonal Yukawa couplings Aq,q, fdaa-, and fuaa are rather small. Due to the 
approximate global U{1) symmetry (|6.1ip that originates from the flavour 
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symmetry A27, the spectrum of inert neutralinos comprises a set of pseudo-Dirac 
states. When the masses of the hghtest and second hghtest inert neutrahnos are 
close, or they form an exact Dirac state, then the decays hi — )• N^Nj^ lead to 
missing transverse energy in the final state. These decay channels give rise to a 
large invisible branching ratio for the SM-like Higgs boson. 

In tables 16. H 16.21 and 16.31 benchmarks i, ii, iv— vi, and ix have almost 
degenerate lightest and second lightest inert neutralinos. In some of these 
benchmark points both lightest inert neutralinos are lighter than mz/2 and as such 
the Z boson can decay into NaNp so that the lightest and second lightest inert 
neutralino states contribute to the invisible Z boson width. In other benchmark 
scenarios both of the lightest inert neutralinos have masses above mz /2 and the 
decays Z — t- N^Np are kinematically forbidden. 

When the LSP and NLSP are close in mass, LSP-NLSP coannihilation might be 
an important factor in determining the dark matter relic density. If this is the case 
then the LSP-NLSP mass splitting should be an important factor. Since 
annihilations of two identical neutralinos are p-wave suppressed, one should 
compare f3Rzii with Rzi2 when trying to determine how important coannihilations 
are, where f3 is the relative speed of the incoming particles, approximately 1/6 
around the time of thermal freeze-out. It is useful to consider the following 
situations: With the LSP and NLSP almost degenerate and with equal 
self-annihilation cross-sections, but a negligible coannihilation cross-section, the 
relic density of dark matter would be twice what it would have have been if the 
NLSP had not been present. If, alternatively, the coannihilation cross-section was 
equal to the self-annihilation cross-sections then the existence of this extra channel 
would lead to a lower present day relic density. In this case it would in fact be equal 
to the relic density calculated in the absence of the NLSP. In this way, in such a 
scenario where coannihilations and self-annihilations are about as important as each 
other, the relic density actually ends up being largely independent of the 
LSP-NLSP mass splitting. 

For the benchmarks i and ii this latter situation is approximately the case and 
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the LSP-NLSP mass splitting turns out not to be an important factor. The mass 
sphtting is in fact small — about half a GeV — but if it were larger and the NLSPs 
were made to have frozcn-out much earlier, the relic density would only be 
decreased slightly (in this case by about a tenth). In benchmark iv, even though 
the LSP and NLSP are close in mass, coannihilations are unimportant due to the 
small value of Rzu- In this case increasing the NLSP mass substantially while 
keeping everything else fixed would lead to an approximate halving of the predicted 
relic density, since the NLSP would have decoupled much earlier than, rather than 
at the same time as, the LSP. The only other benchmark scenario where the LSP 
and NLSP are close enough in mass for coannihilations to be potentially important 
is benchmark ix. Here coannihilation is in fact the dominant process and 
changing the LSP-NLSP mass splitting would have a large effect on the predicted 
relic density. In fact, in this scenario if the NLSP were not present the predicted 
relic density would be within the measured range. 

If the mass difference between the second lightest and the lightest inert 
neutralino is around 10 GeV or more then some of the decay products of a N2 that 
originates from a SM-like Higgs boson decay might be observed at the LHC. In our 
analysis we assume that all scalar particles, except for the lightest Higgs boson, are 
heavy and that the couplings of the inert neutralino states to quarks, leptons, and 
their supcrpartncrs arc relatively small. As a result the second lightest inert 
neutralino decays into the lightest one and a fermion-antifermion pair mainly via a 
virtual Z. In our numerical analysis we did not manage to find any scenario with 
|m^^| — |m^J > 20 GeV leading to reasonable values of fi^/i^. Hence we do not 
expect any observable jets at the LHC associated with the decay of a N2 produced 
through a SM-like Higgs decay. However, it might be possible to detect some 
lepton-antilepton pairs coming from decays of the form hi — > N2Na- In particular 
we hope that pairs coming from such decays of the lightest CP-even Higgs 

state could be observed at the LHC. 

Benchmarks iii, vii, and viii can lead to such relatively energetic muon pairs 
in the final states of SM-like Higgs decays. Since the Higgs branching ratios into 
SM particles are rather suppressed, such decays of the lightest CP-even Higgs state 
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might play an essential role in Higgs searches. 

In addition to the inert Higgs decays, the scenarios considered here imply that at 
least two of the inert neutralino states that are predominantly formed from the 
fermionic components of the inert Higgs doublet supermultiplets, as well as one of 
the inert chargino states, should have masses below 200 GeV. Since these states are 
almost inert Higgsinos they couple rather strongly to W and Z bosons. Thus at 
hadron colliders the corresponding inert neutralino and chargino states could be 
produced in pairs via off-shell W and Z bosons. Since they are light their 
production cross-sections at the LHC would not be negligibly small. After being 
produced, inert neutralino and chargino states would sequentially decay into the 
LSP and pairs of quarks and leptons resulting in distinct signatures that could be 
discovered at the LHC. 

6.3.2 Neutralino and chargino collider limits 

The remarkable signatures discussed above raise serious concerns that they could 
have already been observed at the Tevatron and/or even earlier at LEP. For 
example, the light inert neutralino and chargino states could have been produced at 
the Tevatron. The CDF and DO collaborations have set a stringent lower bound on 
chargino masses using super symmetry searches with a trilepton signal jll9l 1120] . 
These searches ruled out chargino masses below 164 GeV. However, this lower 
bound on the chargino mass was obtained by assuming that the corresponding 
chargino and neutralino states decay predominantly into the LSP and a pair of 
leptons. In our case, however, the inert neutralino and chargino states are expected 
to decay via virtual Z and W exchange, decaying predominantly into the LSP and 
a pair of quarks. As a consequence the lower limit on the mass of charginos that is 
set by the Tevatron is not directly applicable to the benchmark scenarios that we 
consider here. Instead, in our study, we use the 95% confidence level lower limit on 
chargino masses of about 100 GeV that was set by LEP [93] . 

In principle the LEP experiments also set constraints on the masses and 
couplings of neutral particles that interact with the Z boson. As mentioned above, 
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when the masses of A'^i and are below mz /2 they are almost degenerate and 
thus the decays of Z — )• N^Np all contribute to the invisible width of the Z boson, 
changing the effective number of neutrino species A'^^^. The contribution of Ni 
and N2 to N^^^^ is given by 

AiVi#^ = «ii + 2ai2 + 022, (6.28) 



where 
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(6.29) 



All three terms in ()6.28p contribute to AA'^^ only if 2|m^^| < mz- In the case 
where only the Z boson decays into NiNi are kinematically allowed 012 and 022 
should be set to zero. If {mfr \ + Imy^ | < mz whilst 2|myy | > mz then only an and 



2ai2 contribute. 



In order to compare the measured value of A'^^'^^p with the effective number of 
neutrino species in the ii^gSSM N^^^ = 3 + AA^^^^ it is convenient to define the 
variable 



D 



ArLEP /\rLEP exp 
exp 

Clt 



.30) 



The value of D represents the deviation between the predicted and measured 
effective number of neutrinos contributing to the Z boson invisible width. It is 
worth pointing out that in the SM, from ()1.13p . D = 2. In the benchmark scenarios 
presented in tables [UTTl 16.21 and 16.31 the value of D is always less than 3. The plots 
in figure [6T] also demonstrate that there is a substantial region of E'eSSM 
parameter space with ^ < mz/2 and D < 3. This indicates that relatively light 
inert neutralinos with masses below mz/2 are not necessarily ruled out by 
constraints on the effective number of neutrinos set by LEP experiments. Indeed, as 
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argued in section 16.11 the Yukawa couplings fdai3 and fuap can be chosen such that 
Rzai3 are very small. The couplings of the lightest and second lightest inert 
neutralinos to the Z boson are relatively small anyway because of the inert singlino 
admixture in these states. Nevertheless, figure ETT] shows that scenarios with light 
inert neutralinos having masses below mz/2 and relatively small couplings to the Z 
boson can lead to appropriate dark matter densities consistent with observation. 

LEP has set limits on the cross-sections of e'^e~ — )• and e~^e~' — )■ C^C^, 

where predominantly — ^ QQ^^i and Ci — )• qq'Ni respectively [121] . Unfortunately 
the bounds are not directly applicable for our study because OPAL limits were set 
for a relatively heavy or Ci only — greater than about 60 GeV. Nevertheless, 
these bounds demonstrate that it was difficult to observe light neutralinos with 
masses less than about 100 GeV if their production cross-sections 
(T(e^e~ — )■ NaNij) < 0.1-0.3 pb. Since at LEP energies the cross-sections of 
colourless particle production through s-channel 'j/Z exchange are typically a few 
picobarns, the lightest and second lightest inert neutralino states in the E'gSSM 
could have escaped detection at LEP if their couplings Rza/s ^ 0.1-0.3. 

6.3.3 Dark matter direct detection 

Another constraint on the couplings of the lightest inert neutralino comes from 
experiments for the direct detection of dark matter. At the time of the publication 
of paper II the most stringent upper limits on the DMP-nucleon elastic scattering 
spin-independent cross-section came from the CDMS collaboration [122j and from 
the first analysis of 11.7 days of data from the XENONIOO experiment |123) . In the 
low DMP mass region relevant for our study, the most stringent of these was the 
XENONIOO limit. In particular the XENONIOO 11.7 day analysis produced a limit 
on the cross-section of 3.4 x 10~^^ cm^ for a 55 GeV DMP at a confidence level of 
90%. This limit remains fairly constant for lower DMP masses and does not 
increase above about 4 x 10~^^ cm^ for even the lowest LSP masses that are 
consistent with our thermal freeze-out scenario. Currently the most stringent limits 
on the spin-independent DMP-nucleon cross-section come from the more recent 
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analysis of 100.9 days of data from XENONIOO [IM]. The best limit is 
7.0 X lO^^'* cm^, which is for a 50 GeV DMP, again at a confidence level of 

Since in the i^gSSM the couplings of the lightest inert neutralino to quarks, 
leptons, and their superpartners are suppressed, the spin-independent part of the 
A'^i-nucleon elastic scattering cross-section is mediated mainly by t-channel SM-like 
Higgs boson exchange. Thus, in the leading approximation the spin-independent 
part of A'^i-nucleon cross-section in the E'eSSM takes the form [1251 [3] 

4m^m^ I yhi FAr|2 fat-i-s 
f^Sl = — 5 — ^1^11^ I ' {b.61) 



where is the nucleon. 



rur = , and 



27 

q=u,d,s Q=c,b,t 



with 



i^Nfrq = {N\mqqq\N) and 



fro ~ ^ ~ JTq- 

q=u,d,s 



Here, for simplicity, we assume that the lightest Higgs state has the same couplings 
as a SM Higgs boson and ignore all contributions induced by heavy Higgs boson 
and squark exchange). Due to the hierarchical structure of the active Higgs boson 
spectrum and the approximate symmetry this approximation works very well. 
Using the experimental limits set on cxsi and ()6.3ip one can obtain upper bounds on 



^The near degeneracy of the lightest and second to lightest inert neutralinos could result in the 
inelastic scattering collisions in which iVi is upscattered off of a nucleus into N2 and this could affect 
the direct detection of A''i in experiments. However, such processes may take place only if the LSP- 
NLSP mass splitting is less than about 100 keV [126) . In the -BeSSM mass splittings of this order 
are not expected to be typical. In the benchmark scenarios considered in tables 16.11 16.21 and 16.31 
the mass splitting is substantially larger and such inelastic nuclear scattering of A'^i does not play a 
significant role. 
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In tables 16. 1| 16.21 and 16.31 we specify the interval of variations of a^i for each 
benchmark scenario. As one can see from (|6.3ip the value of cisi depends rather 
strongly on the hadronic matrix elements — the coefficients that are related to 
the TT-nucleon a term and the spin content of the nucleon. The hadronic 
uncertainties in the elastic scattering cross-section of DMPs and nucleons were 
considered in ref. [125j. In particular, it was pointed out that could vary over a 
wide range. In tables 16.11 16.21 and 16.31 the lower limit on asi corresponds to = 
while the upper limit corresponds to f,^^ = 0.36 (see ref. From tables [STTl 16. 2[ 
and 16.31 and ()6.3ip it also becomes clear that asi decreases substantially when mh^ 
grows. 

Since in all of the benchmark scenarios presented in tables 16.11 16.21 and 16.31 the 
lightest inert neutralino is relatively heavy, with I^tt-jv^I ~ mz/2, allowing for a small 
enough dark matter relic density, the coupling of A'^i to the lightest CP-even Higgs 
state is always large, giving rise to a A'^i-nucleon spin-independent cross-section that 
is of the order of, or larger than, the 90% confidence level bound of ref. |123j . 

The dark matter scenario detailed in the present and previous chapters is now 
severely challenged by the most recent XENONIOO results ll24j. Although these 
results appear to rule out the E'eSSM as a model of dark matter, it should be noted 
that they do not rule out the model per se. Scenarios similar to those in tables 16.11 
and 16.21 but in which the predicted relic density is somewhat less than the 
measured relic density can be consistent with direct detection experiments, 
although since such scenarios would not completely explain the observed dark 
matter relic density they may be considered less well motivated. In such scenarios it 
needs to be assumed that the majority of the observed dark matter is not composed 
of ii^gSSM inert neutralino LSPs, but is composed of some extra matter beyond that 
of the ^eSSM. 
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6.4 Summary and Conclusions 



In this work we considered novel decays of the SM-like Higgs boson in the -EgSSM. 
Particular attention was given to the dark matter that the model predicts and this 
work also represents a more in-depth study of the inert neutralino and chargino 
sectors of the E'eSSM than the previous study presented in chapter [5] and paper I. 

To satisfy LEP constraints we restricted our consideration to scenarios with 
relatively heavy inert chargino states ^ > 100 GeV. In our analysis we also 
required the validity of perturbation theory up to the GUT scale which sets 
stringent constraints on the values of the Yukawa couplings at low energies. Using 
these restrictions we argued that the lightest and the second lightest inert 
neutralinos are always light — they typically have masses below 60-65 GeV. These 
neutralinos are mixtures of inert Higgsinos and singlinos. In the considered model 
the lightest inert neutralino A'^i tends to be the LSP and play the role of dark 
matter while N2 tends to be the NLSP. The masses of A'^i and N2 decrease with 
increasing tan(/3) > 1 and inert chargino masses. 

Because the lightest inert neutralino states are predomiantly inert singlino in 
nature, their couplings to the gauge bosons, active Higgs bosons, quarks, and 
leptons are rather small, resulting in relatively small LSP annihilation cross-sections 
and the possibility of an unacceptably large dark matter density. In the limit where 
all non-SM states except for the inert neutralinos and charginos are heavy a 
reasonable density of dark matter can be obtained if \'m^_^ ^ \ ~ mz/^, where the 
inert LSPs annihilate mainly through an s-channel Z boson. On resonance an 
appropriate value of ^i-^/i^ can be achieved even for a relatively small coupling of 
the LSP to the Z boson. In order to achieve plausible scenarios consistent with 
both LEP and cosmological observations, requiring jw-^^l ~ fnz/2 if A'^i contributes 
to CDM, tan(/3) cannot be too large. 

The main message arising from this work is that within the dark matter 
motivated scenario although the lightest and the second lightest inert neutralinos 
can have small couplings to the Z boson their couplings to the SM-like Higgs state 
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hi are always large. Indeed, we argued that in the first approximation the couplings 
of A^i and N2 to the lightest CP-even Higgs boson are proportional to Iw-^J/v and 
Irrij^^l/v respectively. Since |m^^_^| must be of order mz/2 in order for the theory 
not to predict too much dark matter, these couplings are much larger than the 
corresponding coupling of b quarks to the SM-like Higgs boson. Thus the SM-like 
Higgs boson decays predominantly into the lightest inert neutralino states and has 
very small branching ratios (2%-4%) for decays into SM particles. 

The most recent XENONIOO dark matter direct detection limits |124j now place 
rather stringent constraints on the E'eSSM inert neutralino dark matter scenario. 
As an explanation for all of the observed dark matter relic density the model now 
looks to be ruled out. There do, however, exist scenarios in which the E'gSSM LSP 
accounts for only some fraction of the observed dark matter that are consistent 
with constraints from colliders and cosmology. 
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Chapter 7 



Dark Matter and Big Bang 

in the EqZ^SSM 

In this chapter, which contains work that was first pubhshed in paper III, with 
the exception of section [731 which contains work that is original to this thesis, we 
introduce a new scenario for dark matter in the -EgSSM in which the dark matter 
candidate is just the usual bino. At first sight having a bino dark matter candidate 
seems impossible since, as already discussed, the lightest inert neutralino mass 
eigenstates, predominantly the inert singlinos Sa, naturally have suppressed masses 
and it is very difficult to make them even as heavy as half the Z mass. To overcome 
this we propose that the inert singlinos are exactly massless and decoupled from the 
bino, which is achieved in practice by setting the Yukawa couplings f(d,u)ai3 to zero. 
This is easy to do by introducing a discrete symmetry under which the inert 
singlet scalars Sa are odd and all other bosonic states are even — a scenario we 
refer to as the ^e^fSSM. 

In the Eq'Z2SSM the inert singlinos Sa will be denoted as a in order to 
emphasise their different (massless and decoupled) nature. The stable DMP is then 
generally mostly bino and the observed dark matter relic density can be achieved 
via a novel scenario in which the bino inelastically scatters off of SM matter into 
heavier inert Higgsinos during the time of thermal freeze-out, keeping the bino in 
equilibrium long enough to give the desired relic abundance. As long as the inert 
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Higgsinos are close in mass to the bino this is always possible to arrange — the only 
constraint being that the inert Higgsinos satisfy the LEP constraint of being heavier 
than 100 GeV. This in turn implies that the bino must also be heavier than or close 
to 100 GeV. These constraints are easy to satisfy and, unlike in the inert neutralino 
LSP dark matter scenario, we find that successful relic abundance can be achieved 
within a GUT scale constrained version of the model — the c^e^f SSM — 
assuming a unified soft scalar mass mo, soft gaugino mass Mi and soft trilinear 
mass ^0 fit the GUT scale. 

It is worth noting that studies of the cE'gSSM [HSIEIES] have hitherto 
neglected to study the full 12 x 12 neutralino mass matrix and only considered the 
6x6 mass matrix of the USSM [4J. Although the question of dark matter was 
addressed in the USSM, the requirement of successful relic abundance was not 
imposed on the cE'gSSM in refs. |128t O [65] even though both analyses considered 
the same 6x6 neutralino mass matrix. This is because it was expected that dark 
matter would arise from the inert sector of the ci^eSSM which was not studied. 
When cosmological constraints on inert neutralino dark matter are included in the 
i^gSSM certain trilinear Higgs Yukawa couplings relevant to the inert sector are 
required to be large as we saw in chapters [S] and El In the c-EgSSM these large 
couplings strongly affect the RG running from the GUT scale and we have not been 
able to show that having inert neutralino LSPs consistent with CDM constraints 
can also be consistent with having universal (GUT scale constrained) soft mass 
parameters. Here we shall consider the cii^gSSM with the full 12 x 12 neutralino 
mass matrix, including both the USSM and inert neutralinos, under the assumption 
that the fermionic components of the inert SM-singlet supermultiplets, the two 
inert singlinos, are forbidden to acquire mass by an extra symmetry of the 
superpotential. In practice there is then a 10 x 10 neutralino mass matrix once the 
two massless inert singlinos are decoupled. 

In summary, the main result of this study is that bino dark matter, with nearby 
inert Higgsinos and massless inert singlinos, provides a simple and consistent 
picture of dark matter in the i^gSSM and is consistent with GUT scale unified soft 
parameters. We also consider the effect of the presence of the two massless inert 
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singlinos in the £^6^2 ^^M on the effective number of neutrinos contributing to the 
expansion rate of the universe prior to BBN, affecting ^He production. Current fits 
to WMAP data |129j favour values greater than three, so the presence of additional 
contributions to the effective number of neutrinos is another interesting aspect of 
the ii'e^f SSM. In practice we find that the additional number of effective neutrino 
species is less than two, due to entropy dilution, depending on the mass of the Z' 
boson which keeps the inert singlinos in equilibrium. 

The ii^e^fSSM is introduced and its neutralino sector is explored in section [TTTl 
The details of the dark matter calculation are presented in section 17. 2i A'cg is 
defined and calculations of its value in the -Eg^f^S^ presented in section [7^ 
Some benchmark points are presented in section [731 The possibility of inert singlino 
WDM is discussed in section [731 and the conclusions are summarised in section [7^ 



7.1 The EeZfSSM 

In the £;6Zf SSM, as well as being invariant under and either or Z^, 
summarised in tables [3^ and [7?H the superpotential of the E'gSSM ()3.6p is also 
invariant under an additional exact symmetry called • Under this symmetry 
only the two inert SM-singlet fields So. are odd. The couplings of the forms and 
therefore forbidden. This means that the fermionic superpartners of Sa — 
the inert singlinos a — are forbidden to have mass and do not mix with the other 
neutralinos. They interact only via their gauge couplings to the Z' boson which 
exist since they are charged under the U{1)n' gauge symmetry. 

All of the exact and approximate discrete symmetries relevant to the -Bg^fS^^ 
superpotential are summarised in table 17.11 

One may worry that the effects of the massless inert singlinos would have already 
been seen in precision measurements from LEP. The inert singlinos, although not 
mixing with the inert Higgsinos as they did in the -EgSSM inert neutralino dark 
matter scenario, still couple to the Zi mass eigenstate because of the non-zero Z-Z' 
mixing angle azz' defined in ()3.2ip . In the following we neglect the kinetic term 
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Table 7.1: The charges of the fields of the EeZ^SSM superpoten- 
tial under various exact and approximate symmetries that the 
superpotential may or may not obey. is already a symmetry 
due to gauge invariance. Either or is imposed in order to 
avoid rapid proton decay. is an approximate flavour symme- 
try. In the iSeZfSSM the extra symmetry Zf is imposed, forcing 
the inert singlinos to be massless and decoupled, i G {1,2,3} and 
a G {1,2}. 

mixing that is expected to be a small effect (see section 13. 2|) . For a given mz' ~ rn^^ 
in (|3.19p rn?zz'^ ^'^'-^ hence also azz'-, is maximised in the limit tan(/3) — )• 00. For 
m^2 ~ '^z' = 892 GeV the maximum value of ^ri^^, is 3270 GeV^ and the 
maximum value of azz' is then 4.15 x 10"'^. The Z\-a-G coupling relative to the 
Z-v-v gauge coupling R is equal to azz' • From ()6.28p the change in the effective 
number of neutrinos contributing to the invisible Z boson width at LEP due to the 
presence of massless inert singlinos is then AA^^^^ = 2i?^ = 2a^^, = 1.72 x 10~^ 
which is well below the experimental uncertainty ft^LEP = 8 x 10^3. When the Z2 

^^cff 

boson mass is large enough to avoid experimental detection limits the contributions 
of massless inert singlinos to the Z boson width and to other LEP precision 
measurements are expected to be within the experimental error. 

7.1.1 The neutralinos of the ^e^f SSM 

The chargino sector of the -Ee^f^SM is unchanged from that of the E'gSSM without 

. The chargino mass matrix is that given in (jS.lip . The same is true for the 
active Higgs scalar masses and mass matrices given in section 13. 3[ The situation in 
the neutralino sector, however, is quite different. 
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In the present study of the -EgZfSSM we define the term 'neutrahno' not to 
include the massless inert singhnos which do not appear in the superpotential and 
are decoupled. The neutralino mass matrix in the interaction basis 



B W-^ 



^3 B' 



-"da 



and again neglecting the small bino-bino' mixing, is then equal to 
/ 
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(7.2) 



where once again QduS U{1)n charges of down-type Higgsinos, up-type 

Higgsinos, and singlinos respectively, as given in table [3?TJ Typically g'^ ~ all the 
way down to the low energy scale. If the soft gaugino masses are unified at the 
GUT scale then we also have M[ ^ ^ M2/2 (see subsection 13. 2. ip . 

The Yukawa couplings in the off-diagonal blocks, marked out by lines, are 
suppressed under the approximate Z^. Given the smallness of these couplings, the 
inert neutralinos in the bottom-right block are pseudo-Dirac states with an 
approximately decoupled mass matrix 



s 

72 



A322 A321 
A312 A311 



\ 



in the basis 



pfO PfO ffO PfO 

-"d2 ^dl ^u2 ^ul 



A322 A312 
Y A321 A311 

They are approximately degenerate with the two inert chargino Dirac states. 
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The top-left block is the USSM neutrahno mass matrix (j5.3p and contains the 
states of the MSSM supplemented by the third generation singhno and the bino'. In 
the case where M^^ ~ M{ is smaU the hghtest neutralino mass state will be mostly 
bino. The bino' will mix with the third generation singlino giving two mixed states 
with masses around Qg g'lS. As Mi w M[ increases, the bino mass will increase 
relative to both the third generation Higgsino mass /u and the inert Higgsino masses 
given approximately by the biunitary diagonalisation of 

At the same time the state mostly containing the third generation singlino will have 
a decreasing mass as M{ increases relative to Qgg'iS. 

7.2 Dark Matter in the cEgZfSSM 

As discussed previously, in section 13.11 there is an automatically conserved i?-parity 
under which the charginos, neutralinos, inert singlinos a, and exotic D and D 
fermions, along with the squarks and sleptons, are all i?-parity odd, i.e. all of the 
fermions other than the quarks and leptons are i?-parity odd. We shall assume that 
the lightest neutralino A'^i is the lightest of all of the i?-parity odd states excluding 
the massless inert singlinos a. However, A^i cannot decay into a via neutralino 
mixing since the inert singlinos are decoupled from the neutralino mass matrix. 
Furthermore, the potential decay A^i — )• aa, allowed by the a-a-B' supersymmetric 
U{l)j\f gauge coupling, is forbidden if A'^i is lighter than the inert SM-singlet scalars 
a. In fact, in this case no kinematically viable final states exist that have the same 
quantum numbers as A'^i. Therefore A^i is absolutely stable and in the scenario 
presented is the DMP. The lightest inert SM-singlet scalar is not stable. There are 
no Yukawa couplings involving Sa, but the inert SM-singlet scalars are able to 
decay via the a-a-B' supersymmetric C/(l)jv gauge coupling. 

In the successful dark matter scenario presented in this section A^i is 
predominantly bino, with at least one of the two pairs of pseudo-Dirac inert 
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Higgsinos expected to be close in mass, but somewhat heavier, in order to achieve 
the correct relic density. This is due to a novel scenario in which the DMP is 
approximately the bino and inelastically scatters off of SM matter into heavier inert 
Higgsinos during the time of thermal freeze-out, keeping it in equilibrium long 
enough to give a successful relic density. In this section we discuss in detail how 
this novel scenario comes about in this model. 

7.2.1 The dark matter calculation 

In the considered model the DMP is not the lightest i?-parity odd state — an inert 
singlino — but the lightest neutralino A^i. We would like to use ()4.3p to describe 
the evolution of i?-parity odd states other than the inert singlinos — generically x- 
In this case we should also include in ()4.3p processes involving a and a particles 
that change the number of x particles by one. Since such processes necessarily 
involve inert SM-singlet scalars a, it is valid to neglect these processes in the case 
where these inert SM-singlet bosons have frozen out long before the freeze-out of 
dark matter. We will call this condition B, to go along with condition A, defined 
in subsection 14.11 and it should be satisfied given our assumption that the inert 
SM-singlet scalar mass eigenstates are heavier than the DMP, since they only 
interact via the heavy Z' boson. Assuming that both conditions A and B are 
satisfied we can use ()4.6p to describe the evolution of the number density n of 
-R-parity-odd states other than inert singlinos x- The value of n after the thermal 
freeze-out of A^i depends on annihilation cross-sections involving A^i and other 
i?-parity odd states close by in mass and n will eventually be equal to the number 
density of DMPs after other x states have decayed to A^i . 

7.2.2 The c^e^fSSM dark matter scenario 

In order to carry out the dark matter analysis in the constrained version of the 
model we have extended the RG code used for the study in ref. [5j |130) to include 
the Yukawa parameters and soft masses of the inert sector of the E'eZf SSM. The 
inputs are K^ij and A333 at the GUT scale, X^a^ at the EWSB scale, s, and tan(/3). 
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as well as the known low energy Yukawa couplings and gauge couplings. Given 
these inputs and the RGEs the algorithm attempts to find points with GUT scale 
unified soft masses mo, Mi 12-, and Aq. The low energy U{1)n gauge coupling g'^ is 
set by requiring it to be equal to the other gauge couplings at the GUT scale, which 
is calculated. 

For consistent points in the E^gSSM the lightest non-inert (USSM sector) 
supersymmetric particle is typically bino dominated. For the C-Eg^f^^^ ^"o.^ 
the same thing. The masses of the inert Higgsino states depend on s and on the 
Yukawa couplings X^ais and in the cSe^f^^^ lightest neutralino can be either 
the bino dominated state or a pseudo-Dirac inert Higgsino dominated state. In the 
latter case we find that the pseudo-Dirac inert Higgsino DMPs coannihilate with 
full-weak-strength interactions and lead to a too small dark matter relic density. In 
the former case the bino DMP normally annihilates too weakly and yields a too 
large dark matter relic density. If, however, there are inert Higgsino states close by 
in mass, they contribute significantly to {(Jv), allowing for the observed amount of 
dark matter. This relies on condition A being satisfied, i.e. the binos being 
up-scattered into inert Higgsinos at a large enough rate. 

Such points with an appropriate dark matter relic density can be found and 
three are presented in section [731 Condition B is satisfied since the inert 
SM-singlet scalars are so much heavier than the DMP and the Z2 boson mass is so 
large compared to the regular Z boson mass. Annihilation and scattering processes 
involving inert SM-singlets and singlinos must contain a virtual Z' boson, which is 
predominantly the Z2 mass eigenstate. 

To test condition A let us compare the rate for binos up-scattering into inert 
Higgsinos with the inert Higgsino coannihilation rate. We shall label the mostly 
bino state A'^i and the lightest pseudo-Dirac inert Higgsino states and N'^. The 
dominant up-scattering diagrams are of the form shown in figure [7T1 

As in section 16.11 we again define Rzij couplings such that the Z-Ni-Nj coupling 
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Figure 7.1: The form of diagrams for the up-scattering of the bino 
dominated DMP Ni off of SM particles X into the pseudo-Dirac 
inert Higgsino states N2 and -/V3. 

N2 



N3 

Figure 7.2: Full- weak-strength coannihilations of the pseudo- 
Dirac inert Higgsino states N2 and N3. 

is equal to Rzij times the Z-v-v coupling. In the -EgZfSSM we can write 

Rzi, = NfNf- ^Fn^, (7.3) 

D=3,7,9 [7=4,8,10 

where is the neutralino mixing matrix element corresponding to mass eigenstate 
i and interaction state a. D and U index the down- and up-type Higgsino 
interaction states respectively. For the pseudo-Dirac inert Higgsino states we have 
7713 ~ —m2 and Rz23 ~ 1, allowing for full- weak-strength coannihilations of the 
form shown in figure [7^21 

Using the notation from ()4.3p . the ratio of the rate for the mostly bino state 
up-scattering into the lightest mostly inert Higgsino state to the inert Higgsino 
coannihilation rate is given approximately by 

/ \ eq eq • v ' -^J 

To give an idea of the size of this ratio, if the SM particle X is relativistic and 
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mi ~ 1712 ~ ^3 then 



T 




RZ23J (|mi|r)3/2exp(-|mi|/T) 




(7.5) 



where again 



T 



X = 



This ratio is expected to be large because of the overwhelming abundance of the 
relativistic SM particle X, but it also depends on Rzi2- The value of Rzi2 depends 
on the Z|^-breaking couplings that mix the upper-left block of the neutralino mass 
matrix in ()7.2p — the USSM states including the bino — with the inert Higgsino 
states in the lower-right block. Since this symmetry is not exact we expect these 
couplings to be large enough such that we can still assume T ^ 1. Explicit 
examples of this parameter are included in table 17.31 in section 17.41 

With conditions A and B satisfied we use micrOMEGAs [91j to calculate the 
dark matter relic density for low energy spectra consistent with the GUT scale 
constrained scenario. The observed relic density of dark matter can arise in this 
model and examples are shown in table [73] in section [731 The most critical factor is 
the mass splitting between the bino and the lightest inert Higgsinos. Too large and 
there are not enough inert Higgsinos remaining at the time of the bino's thermal 
freeze-out to have a significant enough effect. Too small and {av) is dominated by 
inert Higgsino coannihilations, leading to a too small dark matter relic density. 

Since in this scenario the DMP is predominantly bino, the spin-independent 
DMP-nucleon cross-section cisi is not expected to be in the range that direct 
detection experiments are currently sensitive too. The spin-independent 
cross-section of a pure bino is suppressed by the squark masses, but is also sensitive 
to the squark mixing angles [131 1. For each flavour the cross-section vanishes for 
zero squark mixing. Since in practice the DMP will also have non-zero, but small. 
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active Higgsino components, there are also contributions to asi from t-channel 
active Higgs scalar exchange via the bino-Higgs-Higgsino supersymmetric gauge 
coupling. These contributions, though in fact dominant, are quite small, due to the 
overwhelming bino nature of the DMP. Estimates of fisi, using the same proton f^, 
fu , and fs parameters used in the study in ref . |132j , are included in table 17.31 in 
section 17.41 

7.3 The Inert Singlinos and their Contribution to the 
Effective Number of Neutrinos prior to BBN 

In the standard theory of BBN, which happens long after the thermal freeze-out of 
dark matter, the resultant primordial abundances of the light elements depend on 
two parameters — the effective number of neutrinos contributing to the expansion 
rate of the radiation dominated universe A'efr and the baryon-to-photon ratio t]. 

Whilst the primordial abundance of ^He is not the most sensitive measure of ry, 
it is much more sensitive to A'cfr than the other light element abundances. This is 
because prior to nucleosynthesis, when the equilibrium photon temperature is of 
order 0.1 MeV, the number of neutrons remaining, virtually all of which are 
subsequently incorporated into ^He nuclei, is sensitive to the expansion rate of the 
universe, which depends on N^s- The greater the expansion rate, the less time there 
is for charged current weak interactions to convert neutrons into protons. 

The analysis in ref. |133| . using the more recent neutron lifetime measurement 
from ref. [134j . gives Ncs = 3.801q'7q at 2-sigma, implying a more-than-2-sigma 
tension between the measured '^He abundance and the Standard Model prediction 
for of about 3. Although in ref. [135j it is suggested that these errors may be 
larger, similar results are also obtained for the effective number of neutrinos 
contributing to the expansion rate of the universe from fits to WMAP data |129j . 

In the £^6^2^^^ ^^e two massless inert singlinos would have decoupled from 
equilibrium at an earlier time than the light neutrinos, but nevertheless would have 
contributed to A^eff- Exactly when the inert singlinos would have decoupled from 
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equilibrium with the photon depends on the mass of the boson which determines 
the strength of an effective Fermi-hke 4-point interaction vertex that would have 
been responsible for keeping the inert singhnos in equilibrium. The various values 
for A'eff that can be achieved in this model ah fit the data better than the SM value. 

The implications of extra neutrino-like particles present in the early universe 
have long been studied and the methods used in following analysis rely on relatively 
simple physics (see for example ref. |136j ). 

The effective number of degrees of freedom contributing to the expansion rate of 
the universe during the run-up to nucleosynthesis is defined by 

ffeff = ff7 + V85.iVeff(Vll)"/' 

= 2 + V4iVeff(4/ll)4/3, (7.6) 

where g^^^ is the value of g^^ , as defined in ()4.9p , immediately prior to 
nucleosynthesis. Here = 2 is the number of degrees of freedom of the photon and 

= 2 is the number of degrees of freedom of a light neutrino. The three SM 
neutrinos are expected to decouple from equilibrium with the photon at a 
temperature above the electron mass whereas nucleosynthesis does not happen until 
the temperature is below the electron mass. When the photon/electron temperature 
is around the electron mass the electrons and positrons effectively disappear from 
the universqj. Their disappearance heats the photons to a higher temperature then 
they would otherwise have had, but the neutrinos, having already decoupled, would 
continue to cool at the full rate dictated by Hubble expansion. Because of the 
neutrinos' lower temperature at nucleosynthesis they would then contribute less to 
g^^Q per degree of freedom. In ()7.6p A'efr is defined such that in the SM A'^cfr = 3, for 
the three neutrinos decoupling above the electron mass, as we shall see. Extra 
particles, such as the -EgZ^SSM inert singlinos, decoupling above the muon mass 
would have had even lower temperatures at the time of nucleosynthesis and would 
therefore contribute to even less than light neutrinos per degree of freedom, 
much smaller number of electrons remains due to the small lepton number asymmetry. 
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7.3.1 The calculation of A^efr 



In the cEq'Z2SSM there is a typical scenario in which the massless inert singhnos a 
decouple at a temperature above the colour transition temperature (when the 
effective degrees of freedom are quarks and gluons rather than mesons) and above 
the strange quark mass, but below the charm quark mass. This has to do with the 
strength of the interactions that keep the inert singlinos in equilibrium which 
depend heavily on the mass of the Z2 boson. If the inert singlinos do decouple in 
this range then this leads to a definite prediction for Nes- We shall explain why the 
inert singlinos typically decouple in this temperature range in the following 
subsection. For now we derive the value of ATgg in this scenario as an example. 

We shall use the superscript to denote quantities at some temperature r° 
below the electron mass and the superscript e to denote quantities at some 
temperature T*^ above the electron mass and where all light neutrino species are 
still in equilibrium. We shall use the superscript s to denote quantities at some still 
higher temperature above the colour transition and the strange quark mass and 
where the inert singlinos are still in equilibrium. 

At the effective number of degrees of degrees of freedom contributing to the 
expansion rate is 

5eff = 9'y+9g + '^K9e+gfi + gu + 9d + gs + ^9iy + '^ga) 

= 2 + 16 + 7/8(4 + 4 + 12 + 12 + 12 + 6 + 4) = 651/4 (7.7) 
and at it becomes 

gts = 2 + V8(^6 + 4(^pyj (7.8) 

and at it becomes 

5eff = 2 + V8(6(^y + 4(gy^, (7.9) 
taking into account that the neutrinos and inert singlinos now have different 
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temperatures. With no subscript T always refers to the photon temperature, as is 
the notation throughout this thesis. 

From (|4.1U|) . the entropy within a given volume V due to a relativistic {boson, 
fermion} with number of degrees of freedom gi is given by 

S^ = {l,ys}9,^mfV. (7.10) 

Since we are assuming that the inert singlinos decouple before the strange quark 
threshold, in going from to we conserve the entropy in the comoving volume 
separately for the inert singlinos and for everything else. Specifically, for the inert 
singlinos 

{T'fV' = (r|)V^ (7.11) 

and for everything else 

[g-y + 9g + V<9e + 9^. + 9u + 9d + 9s + ^9u)] {T'fV = [g^ + ^Kge + 85.)] {T'fV 

613/4(r")V = W3/i{T^fV^. (7.12) 

This allows us to write 

{T'fV _ fT?Y _ 103/4 _ 43 



(7.13) 



In going from T*^ to we conserve the entropy separately for the neutrinos, for the 
inert singlinos again, and for everything else, giving 

[g, + r/sg,]{T'fV' = g,{Tyv\ (7.14) 
(T^V^ = (Tj^fV^, and (7.15) 



{Tiyv = {T^yv. (7.16) 
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Figure 7.3: Elcctrowcak interactions responsible for keeping the 
light neutrinos in equilibrium in the early universe. For all the light 
neutrinos there are the the processes on the left. For the electron 
neutrinos there is also the additional process on the right. 



This gives us 



" ^ - - and (7.17) 



(7.18) 



57 



T^V 43 43 4 



TV 247 5^ + 7/85e 24711' 



In this case the effective number of neutrinos contributing to the expansion rate 
prior to nucleosynthesis (at T^) is then 

/ 43 x4/3 

A^efT = 3 + 2 (^—j « 3.194. (7.19) 
7.3.2 The inert singlino decoupling temperature 

The light neutrinos are kept in equilibrium via their electroweak interactions. The 
relevant diagrams are shown in figure [731 

We express the cross-section for processes relevant for keeping muon and r 
neutrinos in equilibrium as 



{a.,.,u.v) = k,^lI^X\ (7.20) 



where k2, like ki from ()4.14p . is a constant defined for convenience and 

X' = (i(-^+si„^(%-)))%2^!^.o^03L (7,21) 
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Note that we are using the GUT normahsed U{1)y gauge couphng and so 

92 [5 gi 



cos{'dw) V 3 sm{'dw) 

The cross-section for electron neutrinos witli their extra diagram is then 



(7.22) 



= k^^I^Y', (7.23) 



where 



.2 V2 

We express the number densities of all Weyl fermions still in equilibrium with the 

photon as 

= «efl = n^L = "^iJ-R = "i^e = "i^^ = = hT^ (7.25) 
and the expansion rate is given by 



H = h^/g^T\ (7.26) 



The neutrino decoupling temperature T'^ can then be approximated by 

{ai,v)n^ = H 



{T^^n' = (7.27) 



(7.28) 



with K = ki/k2k3. A more detailed calculation finds that in the SM (with only 

neutrinos, electrons, and photons contributing to g^g) T^'"^'^ ~ 3.7 MeV and 
T'^'' « 2.4 MeV — the muon and r neutrinos decoupling earlier. 

At temperatures above the strange quark mass the processes relevant for keeping 
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^R' /^L' ^_R' ^pLt 

^ij' '^L' ^R' -^L' '^i? 



u1, U'^ji, d\, d^, S^, 




Figure 7.4: Interaction processes responsible for keeping the inert 
singlinos in equilibrium in the early universe. 



the inert singlinos in equilibrium are shown in figure \TM 

The part of the Lagrangian containing all of the fermion couplings in figure [731 
illustrating the relevant U{1)n charges, is 



L 



t c\ 



Qi 4. 4 ) '-^'^'^^'^i^ 



(l)e' 

and the total cross-section taking into account all of these diagrams is then 
(neglecting the small Z-Z' mixing) 



(7.29) 



, t2 , z^ 



«2 — 

m 



(40)2 ' 



where 



(7.30) 



Z^ = (5)2[2(2)2 + 2(1)2 + 3(1)2 + 3(1)2 + 6(1)^ + 6(2)2 + 3(2)2] 

= 1450, (7.31) 
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leading to an approximate singlino decoupling temperature of 





(7.32) 



( 



'(T 



) 



3 



V^e'ff \^z) sin4(^H^) 



(7.33) 



The only unknown variable here affecting the inert singlino decoupling temperature 
is then the Z2 boson mass mz2- Rearranging we find 



We now check which values of are consistent with our assumption that the inert 
singlinos decouple at a temperature between the strange and charm quark masses. 
For T"^ < rric we find that we require < 4700 GeV. For ~ 1000 GeV the 
situation is slightly more complicated. Firstly the temperature of the QCD phase 
transition is not accurately known and secondly the effective number of degrees of 
freedom is decreased by so much after the QCD phase transition that even if the 
inert singlinos were decoupled beforehand the universe may be expanding slowly 
enough afterwards that they could come back into equilibrium. After checking a 
range of scenarios we find that for 1300 GeV < < 4700 GeV our value of 
A'eff = 3.194 is valid. For mz2 ^ 950 GeV the inert singlinos decouple at a 
temperature above the muon mass, but below the pion mass, leading to a larger 
prediction of N^q = 4.373. At the time of the publication of paper II the 
experimental limit in the £^6^2 SSM was mz2 > 892 GeV, from ref. [68], which 
would allow a Z2 boson light enough for us to predict a value for A'cfj = 4.373. For 
Z2 masses in between these ranges the value of A'efT depends on the details of the 
QCD phase transition, but is somewhere between these predictions. For inert 
singlinos decoupling above the pion mass, but after the QCD phase transition, we 
have Neff = 4.065. All of these values are within the 2-sigma measured range 
A'egf = 3.80^Q7Q and closer to the central value than the SM result A'egf = 3. 




3/4 



(7.34) 



7.3.3 A^eff in the ^e^fSSM 
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Benchmark 


1 


2 


3 




tan(B) 


30 


10 


3 




s [TeV] 


5 


4.4 


5.5 


A333 


@ GUT scale 


-0.3 


-0.37 


-0.4 


A322 


^ EWSB scale 


0.1 


0.1 


0.1 


A311 (6 


^ EWSB scale 


0.0293 


0.0403 


0.0399 




@ GUT scale 


0.18 


0.18 


0.23 




Ml/2 [GeV] 


590 


725 


908 




mo [GeV] 


1533 


454 


1037 




^0 [GeV] 


1375 


1002 


413 



Table 7.2: The input parameters of the three cEqIi^SSM bench- 
mark points. 

Since the publication of paper III the limit on the Z2 mass in the E'gZf^^^ 
has increased to around 1350 GeV as discussed in subsection 13.2.21 This leads to a 
concrete prediction of A'cfr = 3.194 assuming that inert singlinos decouple at a 
temperature below the charm quark mass, i.e. < 4700 GeV. 



7.4 Benchmark Points 

In the tables [721 [731 and 17.41 we present three benchmark points in the cEe'^2^SM.. 

For all three points we fix A322 = 0.1 and A321 = A312 = 0.0001 at the EWSB scale. 

For the Z^-breaking couplings we also fix A332 = A323 = 0.012 and 

A331 = A313 = 0.005 at the EWSB scale. At the GUT scale we fix 

^333 = ^322 = H311 and K-^ij = for i 7^ j. The lightest (SM-like) Higgs mass is 

calculated to second loop order. 

We have chosen three points with quite different values of tan(/3) — 30, 10, and 
3. This illustrates the fact that tan(/3) can be quite low in this model since the 
SM-like Higgs mass is not constrained to be less than mz \ cos(2/3)| at tree level as it 
is in the A4SSM. 

The mass of the bino DMP A'^i is not directly constrained to be above above 
100 GeV. However, the lightest pseudo-Dirac inert Higgsino neutralinos A'^2 and 
are almost degenerate with the lightest inert Higgsino chargino Ci and therefore 
these are constrained to heavier than 100 GeV in order to be consistent with LEP 
constraints |93j . Furthermore, the thermal relic DM scenario outlined in section [7^2] 
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Benchmark 


1 


2 


3 




[GeV] 


-1086.7 


-1189.5 


-1405.5 


A322s/\/2 


[GeV] 


353.55 


331.13 


388.91 


A3IIS/V2 


[GeV] 


103.59 


125.38 


155.17 


Ni mass 


[GeV] 


94.07 


114.49 


143.50 


N2 mass 


GeV 

L J 


-105.12 


-126.45 


-156.57 


A^3 mass 


[GeV] 


105.14 


126.47 


156.62 


A^4 mass 


[GeV] 


167.05 


203.19 


255.47 


N5 mass 


[GeV] 


-353.77 


-311.29 


-389.12 


TVg mass 


[GeV] 


353.78 


311.30 


389.13 


Ni mass 


[GeV] 


-1092.5 


-1194.5 


1409.6 


A'g mass 


[GeV] 


1093.3 


1194.8 


-1411.2 


Nq mass 


[GeV] 


-1803.2 


-1572.3 


-1964.7 


iVio mass 


[GeV] 


1899.7 


1688.7 


2109.9 


Ci mass 


[GeV] 


105.04 


126.41 


156.52 


C2 mass 


[GcV] 


167.05 


203.19 


255.46 


Ca mass 


[GcV] 


353.78 


311.30 


389.13 


C4 mass 


[GcV] 


-1094.4 


-1196.1 


-1411.3 




[GeV] 


1850.4 


1628.4 


2035.4 




NeS 


3.194 


3.194 


3.194 






0.112 


0.107 


0.102 




T 


1.1 X 10*^ 


2.3 X 10^ 


2.3 X 10^ 




[cm-] 


1.9 X 10"^'^ 


2.5 X 10- ''' 


1.2 X 10"^*^ 



Table 7.3: The low energy ncutralino and chargino masses and 
associated parameters of the three benchmark points. The DMP is 
the lightest ncutralino A'^i which is predominantly bino in nature. 
There is a nearby pair of inert neutral Higgsinos N2 and N^i, and 
a chargino Ci into which A^i inelastically scatters diiring freeze- 
out, resulting in a relic density consistent with observation. The 
predicted values of and A^eff are also shown, as is the spin- 
independent iVi-nucleon direct detection cross-section asi. 
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Benchmark 


1 


2 


3 


h t TV 1 

/ llldoo 


1 VJC V 1 


122.2 


114.6 


115.3 




[GeV] 


1145 


987.1 


1522 


/13 mass 


[GeV] 


1890 


1664 


2080 


mass 


[GeV] 


2106 


1396 


1675 


mass 


[GeV] 


2103 


1393 


1673 




[GeV] 

L J 


1547 


518 


1084 




[GeV] 


1567 


611 


1156 


m „ , rriTj , 


GeV 

L J 


1561 


599 


1146 




[GeV] 


1483 


503 


1794 


'"■02^ Di 


[GeV] 


1443 


493 


1775 


771 - 


[GeV] 


2864 


2321 


3065 


Til ~ ID ~ 


fGeVl 


9840 


8 






[GeV] 


1 1 00 
ILZZ 


D/O.O 


1 1 1 n 




[GeV] 


1817 


1774 


1707 




[GeV] 


1470 


1069 


1546 




[GeV] 


1838 


1224 


1761 




[GeV] 


1434 


1009 


1512 




[GeV] 


1840 


1226 


1763 




[GeV] 


1748 


1265 


1818 




[GeV] 


1907 


1278 


1820 




[GeV] 


1500 


718.8 


1259 




[GeV] 


1655 


731.3 


1261 




[GeV] 


1708 


949.2 


1473 




[GeV] 


1775 


952.8 


1474 




[GeV] 


1705 


945.6 


1472 




[GeV] 


1774 


949.5 


1472 




[GeV] 


541.3 


626.9 


787.7 



Table 7.4: The remaining particle spectrum of the three bench- 
mark points. 
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requires and not to be too much more massive than A'^i. In practice the A^i 
is predominantly bino and its mass cannot be much less than 100 GeV. In 
benchmark 1, for example, it is 94 GeV. 

Requiring such values for the low energy bino mass Mi and requiring consistent 
EWSB in practice means that the SM-singlet VEV s cannot be too low. This in 
turn means that the Z2 mass is always more than about 1.5 TeV, automatically 
satisfying the most recent experimental lower limit. In these benchmarks from the 
constrained scenario the effective number of neutrinos contributing to the expansion 
rate of the universe prior to BBN N^q therefore takes on the lower value calculated 
in section [73] of around 3.2. This is more consistent with data than the SM 
prediction. 

In all benchmark points A'4 and C2 are predominantly wino. N^, Nq, and C3 are 
predominantly made up of the remaining inert Higgsinos states, with masses around 
^322s/V^, whereas N-j, iVg and C4, are predominantly made up of the active 
Higgsino states, with masses around /x. A'g and A^'io are mostly superpositions of 
the active singlino and bino'. 

The fact that T S> 1 indicates that the inert Higgsino components in the 
predominantly bino state A^i, though small, are large enough such that processes 
involving A'^i up-scattering off of a SM particle into N2 happen overwhelmingly 
more often than neutralino annihilation and coannihilation processes. In this way 
the ratios of the number densities of these particles are able to maintain their 
equilibrium values. 

The spin-independent DMP-nucleon cross-section cjsi, as estimated using the 
results in ref. |131j . is quite small for these benchmarks and is not currently 
detectable by direct detection experiments. This is due to the predominantly bino 
nature of the DMP as well as the large squark masses. 
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7.5 Warm Inert Singlino Dark Matter in the £"6881X1 



A model of dark matter is inconsistent if it predicts a DMP mass that is so light 
that the observed structure of the universe would have been erased. Such light dark 
matter is known as hot dark matter. Although such hot dark matter is inconsistent 
with observation, the dark matter also does not need to be cold, i.e. of a mass such 
that it was non-relativistic at freeze-out, to be consistent with current 
observations jl371 1138] . The intermediate scenario is known as warm dark matter. 
Limits on WDM from WMAP and Lyman-a forest data require the mass of a warm 
thermal relic particle responsible for all of the observed dark matter to be greater 
than 550 eV [137) . This is of the order of various other keV scale lower bounds on 
WDM particles [T38] . 

In the E'gSSM, if the symmetry was only approximate then the only stable 
supersymmetric particle or particles would be either the lightest of or both of the 
two light, predominantly inert singlino states. In this case the inert singlinos could 
form WDM. It is already known that if all of the observed dark matter is made up 
of gravitino WDM decoupling above GeV temperatures with g'^^ ~ 100 then the 
gravitino mass would have to be around 100 eV, contradicting the above limit [137J. 
For inert singlinos decoupling at a temperature T'^ between the strange and charm 
quark masses the singlinos would have undergone even less entropy dilution than 
such gravitinos and their masses would need to be even smaller in order for the 
observed dark matter relic density to be predicted. 

The number density of a single species of {boson, fermion} with temperature 
is proportional to oc T^^ 

n, = 5jl,3/4}Mi:3_ (735) 
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The number density of all three neutrino species today 



3^0 _ 3^Ml(rY--n° 



20 

Sj + 3sj, = -p^s^. (7.36) 



Conserving entropy between inert singlino freeze-out and today then gives 

ngyo = nlV (7.37) 

for inert singlinos and 



n°+3nO)y° = ^nOF° = (n^-2n|)y^ (7.38) 



for everything else which means that today 



n? 20 ni 



11 re^ — 2n| 



(7.39) 



if the inert singlinos were relativistic at freeze-out. 

For the case of two stable inert singlinos with masses mi and m2 the relic 
density today will be given by 



n^h' = nVi + n>2 — = ^ n ■ (7.40) 

^ ^ Pc Pc re^ 



If r^CT = r^DM and = rig = as derived above then this can be rearranged to give 

oPc 1 11 - 2n| 

mi +1712 = ^Dm/i To-nTTTT ^ • (7.41) 

h'^ rvif 20 n| 



Using pc = 1-05 X lO^/i^ eVcm-^ = 410.5 cm"'^ [10], and 

+ 3o..l 111 

(7.42) 



- 2n| _ + gg + 3/4(gg + + + grf + ffs + 3g^) _ 111 



n| 3/2 3 
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gives 

mi+m2 = 57 eV (7.43) 

in contradiction witli data. 

Therefore thermal WDM inert singhnos, hke thermal WDM gravitinos, cannot 
be responsible all of the observed dark matter. WDM inert singlinos with larger 
masses could only be responsible for all of the observed dark matter is there were a 
significant source of entropy dilution reheating the SM matter, but not reheating 
the inert singlinos, after the time of inert singlino freeze-out. Such entropy dilution 
would lower the number density of inert singlinos today relative to the known CMB 
photon number density. Thermal WDM decoupling at g^^ ~ 1000 could also lead to 
a successful WDM scenario, but such a situation is well beyond the framework of 
the i^eSSM, requiring, for inert singlino WDM, a much more massive U{1)n Z' 
boson and, more importantly, the existence of many new degrees of freedom, 
beyond those of the EqSSM, at some high temperature. 

However, the E'eSSM with an approximate symmetry does allow for another 
type of scenario, apart from having lightest inert neutralinos with masses of order 
half of the Z boson mass, in which the supersymmetric particles of the E^eSSM are 
responsible for less than the observed dark matter relic density. Such scenarios are 
consistent with, even if they do not explain, cosmological observations. If the 
symmetry was only approximate and WDM inert singlinos had masses significantly 
less than 57 eV then these inert singlinos would be the only stable supersymmetric 
particles and would contribute less than the observed amount of dark matter. 

7.6 Summary and Conclusions 

The difficulty in making the predominantly inert singlinos states predicted by the 
E'eSSM much heavier than 60 GeV makes them natural dark matter candidates, 
but has also led to a very tightly constrained scenario in which inert neutralino LSP 
dark matter is now severely challenged by the most recent XENONIOO analysis of 
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100.9 days of data. Furthermore we have not been able to show that such a scenario 
could be consistent with having universal (GUT scale constrained) soft mass 
parameters. 

In this work we discussed a new variant of the E'eSSM called the £^6^2 SSM that 
involves a novel scenario for dark matter in which the DMP is predominantly the 
bino with a mass close to or above 100 GeV which is fully consistent with 
XENONIOO data. A successful relic density is achieved via its inelastic up-scattering 
into nearby heavier inert Higgsinos during the time of thermal freeze-out. The 
model also predicts two massless inert singlinos which contribute to the effective 
number of neutrino species at the time of BBN, depending on the mass of the Z2 
boson which keeps them in equilibrium. For > 1300 GeV we find iVgff « 3.2. 

We presented a few benchmark points in the cE'eZf SSM to illustrate this new 
scenario. The benchmark points show that it is easy to find consistent points that 
satisfy the correct relic abundance as well as all other phenomenological constraints. 
The points also show that the typical Z2 mass is expected to be around 2 TeV, with 
the gluino having a mass around 500-800 GeV and squarks and sleptons typically 
having masses around 1-2 TeV. The DMP-nucleon spin-independent direct 
detection cross-sections are well below current sensitivities. 

Although very light inert singlinos in the E^eSSM provide a candidate for WDM, 
in order to account for all of the observed dark matter thermal WDM inert 
singlinos would need to be too light — lighter than would be consistent with other 
cosmological observations. Inert singlino WDM contributing less than the observed 
dark matter relic density would, however, provide another scenario in which the 
^eSSM predicts less than the observed amount of dark matter and is consistent 
with all observations. 
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Chapter 8 

Summary and Conclusions 



In chapter [S] the first study of the inert neutrahno sector of the ii^gSSM is presented. 
It was found that in the ii^gSSM the dark matter naturally arises from this 
approximately decoupled sector. The inert neutralino dark matter scenario was 
studied both analytically and numerically. It was found that in order for the inert 
neutralino LSP not to be too light and singlino dominated, leading to too large a 
dark matter relic density, certain trilinear Higgs Yukawa couplings relevant to the 
inert sector should be large and the ratio of Higgs VEVs tan(/3) should be relatively 
close to unity. If the LSP mass is allowed to increase to around half of the Z boson 
mass then the LSP also contains larger inert Higgsino components and can 
annihilate more efficiently in the early universe, leading to a reduced dark matter 
relic density. Imposing that the LSP has a mass greater that half of the Z boson 
mass, to avoid potential confiict with LEP data, and accounts for all of the 
observed dark matter implies that tan(/3) should be less than about 2, depending on 
the sizes of various Yukawa couplings that one is willing to allow. The inert 
neutralino dark matter scenario relies mostly on parameters that only affect the 
inert sector physics. As a result the parameter space of the MSSM-like sector of the 
-EgSSM is less constrained compared that of the MSSM since in the E'eSSM these 
parameters are not constrained from dark matter considerations. The exception is 
tan(/3) which strongly affects the LSP mass, with the LSP mass being 
approximately proportional to sin(2/3). 
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In chapter [H] a more in-depth study of the inert neutrahno and chargino sectors 
of the E'gSSM, with a particular focus on physics relating to the Higgs boson, is 
presented. The condition that Yukawa couplings remain perturbative up to the 
GUT scale is imposed and the LSP and NLSP masses cannot then be made greater 
than about 60 GeV. Scenarios where the LSP and NLSP masses are around half of 
the Z boson mass are found that produce less than or equal to the observed amount 
of dark matter. It is found that inert neutralino masses below half of the Z boson 
mass can be consistent with LEP data provided that tan(/3) is not too large. In 
plausible scenarios consistent with observations from both cosmology and LEP it is 
found that the couplings of the lightest inert neutralinos to the SM-like Higgs boson 
are always rather large. This means that the SM-like Higgs boson has a large 
branching ratio into invisible final states and this has major implications for Higgs 
boson collider phenomenology. The branching ratio into SM particles is reduced to 
around 2-4%. It also leads to large spin-independent LSP-nucleon cross-sections 
and because of this scenarios in which i^gSSM inert neutralino LSPs account for all 
of the observed dark matter are now severely challenged by recent dark matter 
direct detection experiment analyses. 

In chapter [7] a new variant of the E'gSSM called the Sg^fSSM is presented in 
which the dark matter scenario is very different to the inert neutralino CDM 
scenario and in which the presence of supersymmetric massless states in the early 
universe modifies the expansion rate of the universe prior to BBN. In the dark 
matter scenario the DMP is the bino and a successful relic density is achieved via 
its inelastic up-scattering into nearby heavier inert Higgsinos during the time of 
thermal freeze-out. The nearby pair of inert Higgsino neutralinos form a 
pseudo-Dirac pair with masses approximately equal the corresponding inert charged 
Higgsino Dirac mass and cannot have masses below about 100 GeV. In the 
£;6Zf SSM the two inert singlino states are exactly massless and contribute to the 
effective number of neutrino species at the time of BBN, depending on the mass of 
the Z2 boson which keeps them in equilibrium. For mz^ > 1300 GeV we find 
A'cff ~ 3.2. The dark matter scenario is consistent with having universal (GUT 
scale constrained) soft mass parameters and the DMP-nucleon spin-independent 
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direct detection cross-sections are well below current sensitivities. 

In the ii^eSSM light inert singlinos contribute too much CDM if they are 
non-relativistic at freeze-out — more than the observed dark matter relic density. 
However, in section 1731 we showed if the inert singlinos have masses less than 
around 50 eV then they will contribute WDM less than the observed dark matter 
relic density. 

In the future it would be interesting to study more theoretical aspects of the 
E'gSSM and Sg^f such as how much fine-tuning these models involve and what 
the effects are of potential non-renormalisable terms in the superpotential. At the 
same time, now that the LHC is taking data it is important to study the collider 
phenomenological predictions of these models. Gluino cascade decays in which the 
gluino sequentially decays into the DMP, giving off pairs of fermions at each stage, 
is the subject of a further paper currently in preparation [139] . In this paper we try 
to identify how the ii^gSSM could be distinguished from the MSSM at the LHC. 
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Appendix A 



Weyl, Majorana, and Dirac 
Spinors in 3 + 1 Dimensions 



We represent the Dirac gamma matrices in the Weyl basis 



M _ 



^ 



CT^ 



where cr° = cr° = 1 and a'^ = —a'^ and write a general Dirac spinor 



* = 



with tpL a LH Weyl spinor and a RH Weyl spinor. We write a general 
infinitesimal Lorentz transformation on a Dirac spinor ^ as 



Ai/2(cM,M) 



\ 



\ 



Using the mathematical identity 



2 * 2 
a a = —aa 
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we can see that cr^J transforms as RH spinor and cr^^ transforms as a LH spinor 



= (l + V2id^.£ - V2d^-£) cr'^i^l and 



= (l + V2id2?.a + y2d§_.aj a^il^R. (A.5) 

We therefore define the charge conjugation operation acting on a Dirac spinor to 
be 



I 



(A.6) 



Since — cr^cr^* = 1, applying the charge conjugation operation twice yields = 
as long as wc define w^w^ ~ ^L^*R ~ define ur = —ul = — ^ implying that 

= 1. We define the CP conjugation operation acting on a LH Weyl spinor so 
that the RH spinor 



(A.7) 



and on a RH Weyl spinor so that the LH spinor 



(A.8) 



The gauge and Lorentz invariant part of the Lagrangian for a massive Dirac 
spinor with mass m 

= ^[ia''V^^^ + ^ya''Vi,'il^j^-m'il^y^-m^yj^. (A.IO) 



146 



This may be rewritten in terms of the two LH Weyl spinors ipL ^'^'^ ^ 



C 



D 



ij^ia'^V^il;^ + VR^ia^P^Vp - fmV^R + c.c.) (A. 11) 



up to a total derivative, revealing that a Dirac spinor with mass m is formed from 
two Weyl spinors of the same handedness with a mass matrix 



m 
m 



(A.12) 



The covariant derivative acting on in (jA.lip is the complex conjugate of 
covariant derivative acting on tp^^ so that if 'ijjj^ is in some representation r of some 
gauge group such that 



(A.13) 



then 



(A.14) 



and tjj'^ is in the conjugate representation r. 

The Lagrangian for a single LH Weyl spinor ■0 with mass m 



(A.15) 



may be written in terms of a Major ana spinor 



M 



(A.16) 
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as 



■M 




(A.17) 



The Majorana spinor is nothing but a Dirac spinor that is self-charge-conjugate. 
The Weyl or Majorana mass matrix for a Dirac particle ()A.12p is diagonalised to 



in the Weyl or Majorana mass eigenstate basis. Conversely, if two mass eigenstate 
Majorana spinors have equal and opposite masses and all other quantum numbers 
equal then together they form a Dirac spinor. If the masses of two such Majorana 
spinors are opposite, but not quite equal then together they are said to form a 
pseudo-Dirac state. 




(A.18) 
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Appendix B 

The Pseudoreality of the Spinor 
Representation of SU{2) 

Let the field ip form the representation (2,r) under the gauge group SU{2) G so 
that an infinitesimal gauge transformation acting on (p can be written 

H(da,d/3)(/7 = [l-ida^T" -idl3^T^^^, (B.l) 

where r" = cr"/2 are the generators of 2 and are the generators of r. We also 
define a field (p that is in the representation (2,f) so that 

E{da,d/3)ip = (l + ida'^T''* + id/3^T^*^ip. (B.2) 

The field may be redefined as the equivalent field 2a;r^(^. Using ()A.4p again, 
this time in the form r^r"* = — r'^r^, we see that this field transforms with 

S(da,d/3)2a;r2(^ = 2ujt^ (l + ida^r"* + id^''T^*^<fi 

= (l - ida^r" + zd/3^r,^*) 2tJT2(^, (B.3) 

meaning that that it is in the representation (2,f). We have somehow managed to 
redefine the field so that it transforms in the 2 rather than 2 representation of 
SU{2). Therefore, even though the spinor representation is not real, the 2 and 2 
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representations are somehow equivalent. This representation is sometimes called 
pseudoreal. 

Antidoublet representations of SU (2) can always be redefined to be doublet 
representations. If two fields are in the doublet representation of SU{2), a gauge 
invariant bilinear may be formed by transforming one of the two fields such that it 
is in the antidoublet representation. We thus define the gauge invariant product of 
two doublet representations of SU{2) 




= Ilt2 — tll2 ■ 



(B.4) 
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